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PREFACE TO THE FIRST GERMAN 
EDITION 


The scheme of arrangement of the present volume runs 
along precisely the same lines as those which guided. me 
in preparing the three first volumes of the present series, I 
which aims at giving a thorough introduction to theoretical 
physics. In view of the great area over which theoretical 
optics now extends, we have in this case, too, been able 
to deal with only a meagre selection from the very 
abundant material that is available. The choice made 
was governed first and foremost by the desire to restrict 
the discussion to the framework of the classical theory as 
applied to bodies of continuous space-distribution. This 
enabled me to lay greater stress on the system used in 
arranging and developing the theorems, and on their 
inks with the other branches of theoretical physics. For 
his reason there are numerous references to the preceding 
volumes of this series; the Roman figure I refers to the 
volume on General Mechanics, II to that on the Mechanics 
f Deformable Bodies, III to that on the Theory of 
electricity and Magnetism. 

Although the assumption of matter having absolutely 
mtinuous properties could be maintained in all the 
.needing volumes, it was found necessary here to overstep 
his assumption in treating the problem of dispersion. 
\nd since it is impossible to leave out dispersion in giving 
an account of theoretical optics, I have included a first 
nt induction to the atomic point of view in the last part 
of this volume and have seized this opportunity to make 
an attempt to describe the way in which it links up 
naturally with quantum mechanics. For the circum- 
stance that access to quantum mechanics, as well as to 

v 
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A 

the theory of relativity, can be gained most readily from 
the side of classical theory by making an appropriate 
generalization seems to be indicated not only from the 
didactic point of view but also from considerations of 
concreteness of expression. 

An index to the definitions and theorems is appended. 

Max Planck. 

Berlin - Grime wald , 

July 1027 . 


PREFACE TO THE SECOND EDITION 

The new edition differs from the old only in a few 
minor points and in having a few additions. 

Max Planck. 

Berlin -Grime wald, 

December 1030 . 
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INTRODUCTION 


§ 1. Physical optics is a special department of electro- 
dynamics — namely, that which comprises the laws of 
rapidly varying fields. Its particular significance consists 
in the fact that it represents the branch of physics in 
which the most refined measurements are possible, and 
which consequently enables us to penetrate furthest into 
the details of physical phenomena. At the same time, 
optics presents a clearer illustration than any other 
branch of physics of the peculiar tendency of progressive 
scientific research to leave the original point of departure — 
namely, the specific sense-impressions — and to place 
physical concepts on more objective foundations. For, 
whereas the most important optical concepts, those of 
light and colour, were originally derived from the im- 
pressions on our eyes, these concepts have nothing at all 
to do with the immediate sensation of sight in present-day 
physics, but relate rather to electromagnetic waves and 
vibration periods — a trend of development which has 
justified itself in the abundant fruit which it has borne. 

§ 2. We can progress most easily by linking up with the 
general system of Maxwell’s equations for the electro- 
magnetic field in stationary bodies, particularly if we 
use the special form which they assume for transparent 
and non-magnetic bodies. Since the transparency of a 
body is associated with the condition that no transforma- 
tion of electromagnetic energy into heat occurs in it, all 
transparent bodies are electrical insulators in which the 
vector J of the electric flux vanishes everywhere and at all 
times. Besides excluding conductors, this also excludes 
strongly magnetizable bodies ; for other bodies we may, 
without introducing an appreciable error, identify the 
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magnetic induction B with the magnetic intensity of field 
H. Then, by III (31 ) the field equations assume the simple 
form : 

D = c curl II, il =- c curl E . . . (1) 

together with the supplementary equations III (49) and 
(51): 

div D = 0, div II = 0 . . . . (2) 

Here E denotes the electric intensity of field, Jl the 
magnetic intensity of field, D the electric induction, c 
the critical velocity, all quantities being measured in the 
so-called Gaussian system of units (III, § 7). 

The above system of equations embraces the optics of 
all transparent substances. But the variables that occur 
in them play the part only of auxiliary quantities, since 
they are not directly measured. There is one quantity, 
to determine which is the goal of all optical measurements 
and to calculate which is therefore the proper task of 
every optical theory. This quantity is the vector of the 
electromagnetic flux of energy : 

s = g[E,u] (:i) 

which gives the intensity and direction of the intensity of 
radiation [see III (26)]. 

For the subsequent treatment of these equations we 
have to take into consideration the particular relation 
which connects the vector of the electric intensity of field 
E with the vector of electric induction 1) and which 
endows a substance with its characteristic optical behaviour. 
Accordingly we find it appropriate to divide the material 
into three parts, so that we successively discuss the optics 
of isotropic homogeneous bodies, the optics of crystals and 
the optics of non-homogeneous bodies in which tlxo 
phenomena of dispersion and absorption are included. 



PART ONE 

OPTICS OF ISOTROPIC AND HOMOGENEOUS 
BODIES 



CHAPTER I 


REFLECTION AND REFRACTION 

§ 3. In the case of an isotropic and homogeneous 
substance the relation between electric induction and 
electric intensity of field is expressed by the equation 
III (28) : 

D = e . E (4) 

where e denotes the dielectric constant. The field- 
equations (1) then become : 

eE = c curl H, H = - c curl E (5) 

We shall consider as the simplest particular solution of 
these differential equations the case of a plane wave which 
propagates itself in the direction of one of the co-ordinates, 
say in that of the positive ^-direction. Then all the 
field-components are independent of y and z and we get 
from (5) and (2), since static fields do not come into 
question for optics : 

Ex = 0, IE = 0 

whereas the following differential equations hold for the 
other components : 

dE y _ dlh dEz_ dHy 
6 dt ~ C dx 5 e dt C 8z ’ 

dHy dEx dHx dEy 

dt ~ ° dx’ dt C dx 

Thus there are two pairs of connected quantities among 
these four field-components; namely E y is connected 

5 
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with Hz and Ez with //„, and the same differential equation 
holds for each individual component, namely : 

d^Ej, c 3 ?)~Ej/_ 

~W ~ e ' W ‘ 1 

So if we set : 

7 = ( 7 ) 

it follows from the general integral already derived in II, 
§ 35, for the differential equation ((>) that the most general 
expression for a plane wave which propagates itself in a 
homogeneous isotropic medium in the direction of the 
positive cr-axis is : 

Ex — 0 II x == 0 



where / and g represent arbitrary functions of a single 
argument. 

As we see, both field-strengths are perpendicular to the 
direction of propagation; hence the wave is called 
£ 4 transversal.” It resolves into two components which 
are in general independent of one another, and which 
lie in the direction of the co-ordinate axes. 'In the case 
of each component the electric and the magnetic field- 
strengths are proportional to one another. Their signs 
are determined by the theorem that the directions of the 
electric field-strength, of the magnetic; field-strength and 
of propagation form a right-handed system. 

§ 4. If we now propose to ourselves the question as to 
what is to be measured in this electromagnetic wave and 
which of its properties can hence be ascertained 
objectively, we find the answer in the vector of energy - 
radiation (3) which in the present case reduces to its 
^-component : 

*=-~(EyH, ■ EzHy) — ^ (,/" i- f, 2 ) 


Si 
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Thus in an isotropic body the direction of the energy- 
radiation coincides with that of the wave-normal x, and 
the amount of energy radiated in the time dt through a 
surface F which lies in a wave-plane is : 

S x . F . dt = + g 2 ) . Fdt . (9) 

Since, however, appreciable effects of radiation always 
require a finite time, we never measure the radiation 
vector Sz itself, but rather only its time-integral or its 
mean value in time taken over a sufficiently great interval 
of time T. Hence if we use the following abbreviations 
for the mean values : 

= F l gHt = ^ • • • (10) 

then the amount of energy radiated through the surface 
F in unit time is : 


£(p + g*),F (11) 

which can be recorded by any instrument that takes up 
the radiant energy completely, and provided that it is 
sufficiently sensitive (bolometer, radiometer, thermopile). 

After the total radiation of the wave has been measured, 
its further analysis presents a two-fold problem ; firstly, 
we must separate the two summands / 2 and g 2 from each 
other ; secondly, we must pass from the mean time values 
to the functions themselves ; that is, we must investigate 
the exact form of 'the wave-functions / and g. For this 
purpose we require special optical contrivances the theory 
and action of which we must derive in the sequel. At 
this stage nothing at all can be stated about them. In 
particular there is no reason for assigning any sort of 
periodicity to the functions f and g . Actually there are 
in optics no waves which have a sharply definite period in 
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the mathematical sense, such as we have, say, in acoustics. 
We therefore do best by leaving the question of the form 
of the waves completely aside for the present, taking it 
into consideration only when it becomes really necessary. 
There is only one assumption which we may make from 
the very outset, namely, that the mean time values of 
/ and g vanish, that is : 

/ = 0 and g = 0 , (12) 

For if a wave-function has a mean value different from 
zero, we can imagine the wave in question to be replaced 
by another wave for which the conditions (12) are fulfilled, 
with a statical field superposed on it, the field being 
characterized by the mean value, which is not equal 
to zero. The presence of this field can be made mani- 
fest by its ponderomotive action on a charged test-body 
(electron) and can so be separated from the true optical, 
wave. 

§ 5. A plane-wave of unlimited cross-section cannot, 
of course, be realized in nature. Nevertheless wo can 
produce waves which approximate appreciably to the 
character of plane-waves. For let us imagine a point- 
like source of light which begins to emit light at a definite 
moment of time — say when t = 0. Then, since the 
surrounding medium has been assumed to be homogeneous 
and isotropic, the light will propagate itself uniformly in 
all directions. The bounding surface which has been 
reached by the light after a definite interval of time is 
called a wave-front. There is thus a wave-front corre- 
sponding to every moment of time, and the whole of the 
surrounding space is hence filled by the system of successive 
wave-fronts which enclose one another. In the present 
case these wave-fronts are obviously spherical surfaces 
which surround the source of light concentrically, and so 
a small portion of a sufficiently great spherical surface can 
be regarded to a sufficient degree of approximation 
as a plane wave-front or a wave-plane. Its normal is the 
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First body 


Second body 


corresponding radius of the sphere, and the radiation 
vector S points in the same direction. 

§ 6. Let us now investigate the phenomena that occur 
when the plane-wave (8) falls on the plane bounding 
surface of a second isotropic body. We shall take the 
normal of this bounding plane, the so-called incident 
normal, as the £-axis of a new co-ordinate system directed 
towards the interior of the second body, whereas the origin 
0 of the ^-system is to be coincident with the origin of 
the ^-system. Without reducing the generality of the 
case we can then make the ?/-axis and also the rj- axis lie 
in the plane defined by x and £, the so-called incident 
plane, and take this as the plane 
of Pig. 1. Here all points for 
which | < 0 denote the first body 
(on the left), from which the 
wave (8) comes, and all points 
for which £ > 0 denote the second 
body (on the right) ; the points 
for which £ = 0 (the ^-axis) 
constitute the boundary plane. 

The a>axis is the direction of the 
ray which comes from the first 

body— -that is, from the left-hand side; it makes the 
angle 6 with the incident normal £. The y - axis denotes 
the wave-plane of the incident ray; this wave-plane is 
perpendicular to the plane of the figure and also makes 
the angle 6 with the boundary plane. It has been 
omitted in the figure so as not to multiply the directions 
to be shown unnecessarily. The 2 -axis coincides with 
the £-axis and points from the plane of the figure towards 
the observer. 

We base the solution of the problem before us on the 
reflection that every system of waves which satisfies the 
differential equations in the interior of the two bodies and 
also the boundary conditions, represents a process which 
is possible in nature. 

In order to have the differential 
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the second body we imagine a plane-wave in it also, after 
the model of equations (8), which has the rav-direction x x 
(see Fig. 1), inclined at an angle d 1 to the f-axis, and the 
wave-plane y^, where we shall again suppose z x to co- 
incide with z and £. Then the equations (8) hold for the 
six field-components E X1 , E n , E n , n xi , II n , Jt n , except 
that the co-ordinate x x now occurs in place of x on the 
right-hand side of these equations, the functions j\ and g x 
replace the wave functions / and g, while the constants 
e and q are supplanted by the dielectric constant e x and, 
by (7), the velocity of propagation : 





( 13 ) 


in the second body. 

But this assumption does not suffice. For by III, § 0. 
the boundary conditions require that for £ = 0 the values 
of the tangential field-components — that is, the quantities 
Erj, E%, H m H ( : — are coincident in both bodies. This gives 
four equations connecting the wave-functions ; to satisfy 
them, however, we have only the two functions j\ and g x 
available, since the functions / and g are initially given. 
To generalize our initial assumption still further, therefore, 
we assume a second wave in the first body; this wave is, 
of course, also represented by the equations (8), except 
that it has a different ray-direction x' 9 which we shall 
assume to make an angle d' with the f-axis (see hig. I), 
and has the wave-plane y'z', where again z r = z. The 
six field-components E X ’> E, r , E z >, ./fr, H, r , H z > are given 
by the equations (8), if we substitute in them the wave- 
functions f and g' and the co-ordinate x f , the constants 
e and q remaining the same. 

We have now approximately generalized our assumption 
for the interior of the two bodies, and can proceed to set 
up the boundary conditions. In the first body there is 
an electromagnetic field which results from the super- 
position of the two plane waves that we have assumed. 
Hence, remembering that the field-components E x , 
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Ex', M x > vanish, we get for the field-components which 
interest us in the first body : 

JB, = Ey cos 9 + Ey. cos 6 ' = . f + .f 

V e ye 

jEj = JS* + JB*. = A • £ + A • 0' 

V<= V e 

Hr? = H y cos 0 4- H y . cos 9' — — cos 9 .g — cos 0' . g r 

h c = a + »>=/ + /' 

On the other hand, we obtain for the second body 
(i > 0), remembering that the field-components Ex, and H x , 
vanish : 

Er, = Ey, COS 9, = ^-7l 

E l = E n= ^Si 

H v = H yi . cos 0 X = — cos 0 X . </ x 

H ? = H n =/i. 

Hence if we use the abbreviation : 



we must have for the boundary plane | = 0 ; 

cos 9 ./ + cos 0'. /' = C -^l -/i 


(14) 


cos 9 . g + cos 0' . cf = cos 0 X . </ x 
/+/'=/! 

These four equations comprehend all the details of the 
theory of reflection and refraction. As we see, they fall 
into two groups, one of which contains only the /-waves, 
and the other only the g- waves. Thus these two kinds of 
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waves behave quite independently of one another ; each 
obeys its own laws. 

§ 7 . By means of the last four equations wc first cal- 
culate the unknown wave-functions /', f x , g', g x from the 
given wave -functions / and g. We get : 


r 


f 1 = 


n cos 6 — cos 9 X 
cos 8 1 — n cos ( " 

n [ cos 6 — cos 8') 
nns 0, — n cos O' 


-f-r-f 


(15) 


(10) 


cos 9 — n cos 9 X 
n cos 6 1 — cos 9' ® 


a.g 


9i = 


n(cos 9 — cos O ') __ 

n cos 9 1 — cos ~ a i-. 


(17) 


(IS) 


As for the arguments of these functions, we have : 


t — ~ in /and <7 
t -*Hnf 1 a,ndg 1 


t — — in/' and g' 


And £ = 0 everywhere, so that in transforming to the 
co-ordinates f, 77 , £ we have : 

x — r] sin 9, x x = rj sin 0 X , a;' — rj sin 0' 

which makes the arguments assume the values : 


2 r) sin 0 ^ 77 sin 0 X 77 sin 0 ' 

<T~’ ?i ’ ~~ <1 

Since the functional equations (15) to (18) must be 
satisfied for all times t, and for all points 77 of the boundary 
surface, it follows that these three arguments must bo 
equal to one another — as can also be seen directly if wc 
differentiate one of the functional equations by parts first 
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with respect to t and then with respect to tj, and divide 
the resulting equations by one another. We get : 

sin 6 _ sin 6 X _ sin 9' ,, Q * 

<7 ~~lh~ ~~q~ 


and hence arrive at the law of refraction 


sing _ q 
sTn0i _ ft 


A 




( 20 ) 


and the law of reflection : 

8' = it — 6 ( 21 ) 

If we call the angle which the reflected ray makes with 
the reversed incident normal the angle of reflection, then 
the angle of reflection is equal to the angle of incidence. 

§ 8. Snell’s law of refraction (20), which states that the 
ratio of the sine of the angle of incidence 6 to the sine 
of the angle of refraction 0 1 is equal to the refractive 
index n of the second body with respect to the first or 
to the ratio of the velocities of propagation q and q x> has 
been accurately confirmed by innumerable measurements. 
The refractive index of a substance is usually referred 
to air as the first substance. Thus the refractive index 
of water is equal to 1*3, that of glass to 1*5. We then 
obtain the refractive index of a substance with respect 
to any other substance by writing down the ratio of their 
refractive indices with respect to air. If we exchange 
the substances, the refractive index assumes the reciprocal 
of its previous value. Accordingly, the refractive index 
with respect to a vacuum — the so-called “ absolute ” 
refractive index — is the product of the refractive index 
with respect to air and of the absolute refractive index 
of air — namely 1*0003; as we see, its value differs in 
most cases only inappreciably from the ordinary refrac- 
tive index. 

If we allow the angle of incidence 8 to vary from 0 
(normal incidence) to ^ (grazing incidence), the angle 
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of refraction 0, increases from 0 to sin 1 - (limiting 

angle). But there is a point of fundamental importance 
which must not be overlooked. It is only when n>l, 
or if, as we say, the second substance is optically denser 
than the first, that the limiting angle is real. Then the 
angle of refraction d 1 is always smaller than the angle of 
incidence 0 — that is, the ray is bent towards the incident 
normal by the refraction, and the limiting angle denotes 
the greatest value which the refractive index can assume 
at all. But if n< 1 — that is, if we exchange the two 
substances with each other — the angles of incidence and 
refraction also exchange their roles, and the angle of 
refraction becomes greater than the angle of incidence ; 

it attains the value \ only when the angle of incidence 

has reached the value of the limiting angle. If the angle 
of refraction is allowed to go beyond the limiting angle, 
then (20) leads to an imaginary value for the angle of 
refraction, and the solution which we have found for the 
problem of refraction becomes meaningless. As there 
is nothing to prevent our giving the angle of incidence 

any arbitrary value between 0 and a special question 

arises here, which we shall, however, deal with on a biter 
occasion (§ 12) ; for the present we shall restrict ourselves 
to considering those cases for which the law of refraction 
yields a real value for the angle of refraction 6 1 . 

§ 9. But the electromagnetic theory of the refraction 
of light states more than that the refractive index is 
independent of the value of the angle of incidence ; it also 
tells us the value of the refractive index. For by (20) 
this is equal to the square root of the ratio of the dielectric 
constants, or, if we take as our basis the absolute refractive 
index : 

n = Ve x ( 22 ) 

If we compare this relationship with observed facts, we 
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find, in general, crass disagreement. For example, for 
water n = 1*3, while = 80. But even apart from 
this the fact that (22) cannot be generally valid follows 
from the fact that by definition the dielectric constant c 
is independent of the form of the wave-functions / and g , 
whereas the refractive index n , in the case of all substances, 
depends more or less markedly on the form of the light- 
waves, that is, on the colour of the light. This phe- 
nomenon, dispersion, long constituted a serious obstacle 
to the acceptance of Maxwell’s theory. If we wish to 
take adequate account of it in the theory here described, 
nothing remains but to conclude that the fundamental 
assumption which was introduced at the beginning of 
this chapter into the field-equations for the optics of 
homogeneous and isotropic bodies — namely, the relation 
(4), which states that the electric induction is proportional 
to the electric intensity of field — does not in general 
correspond with reality in the case of rapid optical 
vibrations. To obtain a satisfactory theory of dispersion 
we shall therefore have to replace this relationship by 
one that is more general. This will be done in the third 
part of the present volume, where it will be found that 
this generalization will have to be based on the circum- 
stance that in the case of refined optical phenomena in 
nature the assumption that matter is absolutely con- 
tinuous and homogeneous is no longer justified, but must 
be modified by the introduction of characteristic structural 
properties to a certain extent. 

If this view is correct, an important significance will 
still have to be attached to the relation (22) — namely, 
that of a limiting law which is the better fulfilled the 
less the dispersion makes itself observed. If we carry 
out a test in this direction, the relationship in question 
is found to be definitely confirmed. For the substances 
which disperse least are gases, and the earliest measure- 
ments, by L. Boltzmann, have accurately confirmed the 
formula (22) in their case. A particularly noteworthy 
feature is the exact quantitative parallelism between the 
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dependence of the refractive index and of the dielectric 
constant on the pressure in the case of gases, and this occurs 
in the sense of equation (22). Hence, with reference to 
this equation, we are right in speaking of a far-reaching 
confirmation of the electromagnetic theory within the 
admissible range of application. 

§ 9a. But besides giving the directions of the reflected 
and the refracted rays, the theory also gives the form 
of the reflected and the refracted rays, by demanding 
that the wave-functions in question shall be proportional 
to the corresponding wave-functions of the incident wave. 
If in the formulae ( 15 ) to ( 18 ) we replace the refractive 
indices n according to (20) by the angles 0 and 0 V and 
the angle d' by tt—6, the constants of proportionality 
assume the following values : 

for the reflected wave (/', g') : 

tan (9 — 9 X ) sin (0 — 0,) 

'‘-fi5rr+8j' ,7 = »i"(« MM ' ■ < - i) 

for the refracted wave (/, , </,) : 

sin 20 _ sin 2d 

/J ' 1 ~ sin (0 -|- 0j) cos ( 0 — d t )’ ‘ sin ( 0 | 0 t ) ' 

According to these formulas (known as Fresnel's 
formulas) there is a fundamental difference between the 
two wave-functions / and g, which corresponds to the 
physical circumstance that, according to (8), in the ease 
of the /-wave the electric intensity of field lies in the 
plane of incidence, whereas in the ease of the </-w ave the 
electric intensity of field is in a direction perpendicular 
to the plane of incidence. The coefficients p correspond 
to the former, the coefficients a to the latter. 

To test the theory we have to measure the radiant 
energy. Let us first consider the reflected ‘ware. From 
equation ( 11 ), using ( 15 ) and ( 17 ), we get for the ratio 
of the intensity of radiation of the reflected wave to 
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the intensity of radiation of the incident wave, that is, 
the “ reflection coefficient 55 : 

fl + F _ 

f + f P P + ~f 

As we see, the reflection coefficient p lies between the 
values p 2 and a 2 which correspond with the limiting cases 
where one of the two waves / and g is vanishingly small. 
In general, by measuring p we find for the ratio of the 
intensities of radiation of the waves f and g : 

/ 7; ^ 2 = 7^772 ( 2S ) 

This relation may be tested experimentally by measur- 
ing the reflection coefficient p of a definite incident 
wave (/, g) for different angles of incidence 9 and, after 
introducing the values of p, and or calculated from (23), 
investigating whether the same value comes out for the 
intensity ratio (25) for all angles of incidence. This has 
been confirmed experimentally in all cases where the 
surface of the body used has been sufficiently smooth 
and where sufficient precautions have been taken to 
avoid the chemical impurities that often contaminate 
surfaces, including those due to the use of polishing 
agents. 

Light derived from a body which has been made, to 
glow by having its temperature increased is called e£ natural 
light. 55 For this, measurement gives, as we should 
expect, p = g 2 , and hence : 



The dependence of the reflection coefficient p on the 
angle of incidence 9 is obtained from (23). For normal 
incidence we have, since 9 and 6 1 are infinitely small : 


9-9. 

p= we; 


'n — 1 
n + 1 


(27) 
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The reflection is the greater, the more the refractive 
index n differs from the value unity, that is, the more the 
two bodies differ optically. For grazing incidence, 

0=5, or, for the limiting angle, B x = ~, the reflection 

coefficient attains its maximum value 1 ; in that case no 
radiation penetrates into the second body at all. 

But the inferences to be drawn from the theory and its 
possible applications go much further. For if we consider 
the constitution of the reflected radiation we find that if 
the incident light consisted of natural light, the reflected 
light does not, but that rather, by (15) and (17) : 


p :(f , 


g, 2 __ cos 2 (0 + 9 X ) 
a 2 ~~ cos 2 (6 — 6 X ) 


(28) 


It is only for normal and grazing incidence (or in the ease 
of the limiting angle, respectively) that this ratio is equal 
to 1 and that the reflected radiation consists of natural 
light — which is obvious at once, since in the f irst ease, the 
incident plane is indeterminate, and hence there is no 
physical difference between the waves f and (f, whereas 
in the latter case all the light is reflected. But in general 
the two components have different intensities that is, 
the reflected radiation is “polarized,” the intensify of\the 
f- wave, whose electric intensity of field lies in the incident 
plane, being smaller than that of the fy-wave. I n fact, for 
the case : 


-f 0 X = or tan 0 n (29) 


that is, when the reflected ray is perpendicular to the re- 
fracted ray, the quotient (28) is equal to zero, so that 
the /-wave is absent altogether in the reflected light, and 
the gr-wave alone remains. The radiation is then said to 
be “ completely ” polarized “ linearly,” because in the 
present case the electric and the magnetic intensity of 
field of the wave have perfectly definite directions. The 
angle of incidence defined by (29) is called the “ angle of 
polarization,” the plane of incidence the u plane of polariza- 
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tion ” and the reflecting body (mirror) is called the 
“ polarizer/ 5 If we call the plane defined by the direction 
of the ray and the electric intensity of field the “ plane 
of vibration 55 of the ray, then the plane of vibration in the 
light reflected at the angle of polarization (29) is constant 
and perpendicular to the plane of incidence (Brewster’s 
Law). This law can be directly tested experimentally by 
allowing the linearly polarized ray to be reflected from a 
second body of the same substance — namely, at the 
polarizing angle — but in such a way that at the second 
reflection the incident ray assumes the role of the /-wave — 
that is, so that the plane of vibration coincides with the 
new plane of incidence. This occurs when the plane of 
incidence at the second reflection is perpendicular to the 
plane of incidence at the first reflection. In this case 
no light is reflected at all (Malus’s mirror experiment). 
The second reflecting body then acts as an <£ analyser.” 

Thus the reflection at the angle of polarization forms, 
quite generally, a means of measuring the two wave- 
components / and g present in any given radiation. For 
by measuring the intensity of the whole reflected light we 
obtain, since f vanishes : 


U 1 = tfV 


n A • 




n 2 + 1 


r ■ 


(30) 


and from this the intensity of the g - wave ; and in the same 
way, by reflection at an incident plane turned through 

the intensity of the /-wave is obtained. 

When the reflection takes place at any arbitrary angle 
of incidence 9 the reflected ray is only partially polarized. 
The intensities of the two wave-components f 2 and g' 2 
are then obtained by multiplying the corresponding 
components f 2 and g 2 in the incident wave by the corre- 
sponding reflection coefficients. For the light whose 
vibrations are in the plane of incidence (/-wave) we have : 

0 tan 2 (6 — 0 X ) 

Pi\ = F* = ta^Jf+d^j 


( 31 ) 
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For the light whose vibrations are perpendicular to the 
incident plane (g-wave) we have : 


sin 2 (9 - e x ) 
sin 2 ( 9 + 0J 


(32) 


§ 10. Exactly similar considerations apply to the 
radiation transmitted by the boundary surface. They can 
be reduced to those of the preceding paragrapli by re- 
collecting that according to the principle of the conserva- 
tion of energy the transmitted radiation is equal to the 
excess of incident radiation over reflected radiation. 
For the boundary conditions (cf . Ill, § 6) provide that no 
energy of radiation still becomes lost. If we take as the 
definition of the “ coefficient of transmissibility 55 the 
fraction of the incident radiation which passes through 
the boundary surface into the second body, the trans- 
mission coefficient for the light which vibrates in the plane 
of incidence (/-wave) is : 


2 sin 29 . sin 2 0 X 

^ ~~ sin 2 (9 + 0/) cos 2 (9 - 


(33) 


whereas for the light which vibrates perpendicularly to 
the incident plane (g-wave) : 


1 " P± = 1 


sin 2 9 . sin 2 9 X 
sin 2 (0'-|- 0 X ) 


(34) 


Of course, we can also calculate the transmission co- 
efficient directly from the ratio of the intensity of the wave 
fi or 9i which penetrates into the second body to that 
of the incident wave / or g. Now we may not, by referring 
say to (16) and (18), set the transmission coefficient in 
question equal to /x x 2 or oq 2 , but must revert to the ex- 
pression (11) which gives the energy radiated per unit of 
time through a surface F which lies in a wave-plane. This 
is clear if we consider that in the passage into the second 
body both the velocity of propagation q and also the 
extent of the surface F , which represents the normal 
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cross-section of a selected cylindrical beam of rays 
change. A simple geometrical reflection shows that : 

(3S) 

cos e ~ cos B t ' ' ' 1 ‘ 

For this is the size of the surface which is cut out of the 
limiting plane (£ = 0) on the one hand by the cylinder of 
the incident rays, on the other hand by the cylinder of 
transmitted rays which results from it. 

Hence we get as the ratio of the transmitted energy of 
radiation to the incident energy for the /-wave by (11), 
(20), (35) and (16) : 

Vif/Fi „ sin 0 X 2 cos 9 1 
qp]p ~~ sin 0 cos 9 

and correspondingly for the tj - wave 

( h ( Ji^\ = sin 9 1 ' 2i cos 6 X 
qg*F sin 6 Gl cos 6 

and these expressions, by (24), agree with (33) and (34). 
The ratio of the two components of radiation in the trans- 
mitted radiation is : 


,/i 2 MF i 

gf ^l 2 cos 2 (<9 - 0^ 


(36) 


Since this expression can become neither zero nor 
infinite, the transmitted radiation is never completely, 
but only partially polarized in the case of incident natural 
light ; and the component whose vibrations lie in the plane 
of incidence is the more intense ; this is the reverse of what 
happens in the case of reflected radiation. 

Finally we must refer to a relationship which expresses 
a general law of optics. Both the equation (20) of the 
law of refraction and the values of the reflection and the 
transmission coefficients remain unaltered if we exchange 
the two bodies and simultaneously the angles 6 and 6 X — 
that is, if we allow the wave to fall in the reverse direction 
from the second body on to the boundary surface of the 
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first body. In other words, light does not only pursue the 
reversed path, but is also divided in the same ratio into 
reflected and transmitted radiation. A boundary surface 
of two bodies reflects equally well in both directions and 
also transmits equally well. This law is an application 
of a law of reciprocity first enunciated by Helmholtz, 
which states that the weakening of intensity which a 
definite ray of light undergoes in its passage through 
arbitrary different bodies in consequence of reflection, 
refraction, diffusion and absorption is independent of the 
direction in which the ray travels. 

§ 11. Now that we have seen that the components / 
and g of a given plane wave (8) which vibrato per- 
pendicularly to one another can be separated from one 
another by measurement and have seen how this can bo 
done, the further question arises whether the values of the 
two intensities/ 2 and g 2 also determine the intensity of the 
component of radiation which vibrates in another direction 
perpendicular to the direction of propagation, analogously 
to the way in which two mutually perpendicular com- 
ponents of a force also determine the component in any 
other direction. We shall find that this is not the case 
and shall see wherein the characteristic difference lies. 

To decide this question, we again consider the most 
general case (8) of a plane wave, but for convenience of 
notation, we shall take the z - axis as the direction of pro- 
pagation instead of the rr-axis, so that the ay-plane be- 
comes the wave-plane of the radiation. Let the /-wave 
vibrate in the ^-direction, and the g - wave in the y-diroction. 
If we now replace the axes x and y by two other mutually 
perpendicular axes x' and y' that lie in the same plane and 
that make an angle <f> with the former axes, then since 
the wave-functions /' and g ' by (8) transform like the 
components of the field-intensity, they become, when 
referred to these new axes : 

f = /cos<£ + gshi<f> . 

g’ = —f sin <f> + g cos 


( 37 ) 

( 38 ) 
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and hence the corresponding intensities : 

f 2 = p cos 2 cf> + 2 \fg sin <f> cos <£ + g 2 sin 2 <f> 
g' 2 = /' 2 sin 2 — 2fg sin <f> cos cf> + g 2 cos 2 <£. 

Adding these expressions we get : 

/* + **-/* + ? 

that is, we always get the total intensity of radiation of the 
whole wave if we add together the intensities of two 
components which vibrate in perpendicular directions; 
this is, indeed, obvious from the fact that the original 
^-system can be chosen quite arbitrarily. But the 
intensity of any single component, say that whose direction 
of vibration makes the angle </> with the x-axis, does not 
alone depend on the intensities/ 2 and g 2 , but also on a third 
quantity, fg. So if we use the abbreviations : 

P^A, g% = B, fg =* C (39) 

then the component of the intensity of radiation vibrating 
at an azimuth </> assumes the value : 

/' 2 = J</, = A cos 2 <f> + B sin 2 (j> + 2(7 sin <f> cos (f> . (40) 

Accordingly, the intensity of radiation is not represented 
by a vector, but by a tensor of the second degree (II, § 13, 
§ 20) — namely, by a plane tensor which has only three 
components and which is characterized by not having the 
quantity (40) dependent on the choice of the co-ordinate 
system. The component C can have any value between 
— \/AB and + \/AB, but it cannot overstep this range, 
since J is positive for all values of J. 

For one of the limiting cases, G = ± V AB, J by (40) 
becomes a perfect square : 

J = (VA . cos </> ± VB . sin^) 2 . . (41) 


So J = f ' 2 vanishes for the azimuth tan = T 



and the whole radiation reduces itself to a single component 
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g“i which vibrates in the perpendicular azimuth. In this 
case the light is linearly polarized. From /'- = 0 it 
follows that/' = 0 and by (37) : 

-^ = - tan / = ± (42) 

that is, the two wave-functions / and g bear a fixed ratio 
to each other which is independent of their argument. 
In such a case, where the one wave-function is completely 
determined by the other, the two corresponding com- 
ponents of radiation are said to be cc completely coherent.” 
Thus linear polarization always implies complete coherence 
of the two components into which the polarized ray can 
be resolved. It follows directly from the relationships 
(15) to (18) that linearly polarized light remains linearly 
polarized after refraction and reflection, since then the 
ratios f x : g x and f : g' are also constant. 

In general the light is only partially polarized, so that 
the two components are partly non-coherent. Natural 
light forms the opposite limiting case, for which the 
intensity of radiation J (j) is quite independent of the 
azimuth of vibration <j>. By (40) this not only 
leads to A = J5, but also to C = 0. As may easily be 
seen, the latter condition is always fulfilled in actual fact, 
if the two vibration components / and g are completely 
non-coherent — that is, fully independent of each other. 
For this means that to a given value of / there belongs not 
one definite value of g , but an enormous number of values 
of g, and that the mean of all these values, g, is inde- 
pendent of the value assumed by /. From this there 
follows as the mathematical expression for the complete 
non-coherence of /and g : 

N — f -9 (43) 

and by (39) and (12) : 

0 = fy = 0 (44) 

Since the functions of / and g depend uniquely on their 
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common argument, a necessary preliminary condition for 
the non-coherence of / and g is, that corresponding to a 
definite value of / or g there must be an enormous number 
of values of the argument, or that the value in question 
of the wave-function must recur enormously often in the 
course of time. 

Although the relation (44) has been shown to be a 
necessary consequence of the complete non-coherence of 
the two wave-functions / and g , this relation is by no means 
sufficient to ensure non-coherence. Indeed, we shall 
soon see (§ 20) that the components of /and g can actually 
be completely coherent even when the equation (44) is 
simultaneously fulfilled. Of course / and g are not pro- 
portional to one another, as in (42) — that is, the light is not 
linearly polarized — but nevertheless it is completely 
polarized, since the one component is completely deter- 
mined at every moment by the other component. 

If we now again turn to consider the general case, we 
easily see that the intensity of radiation J ({) has by (40) a 
maximum and a minimum, both of which result from the 
equation : 


and occur for the angle which is obtained from : 

tan 20 = .... (45) 

These two directions, which are mutually perpendicular, 
are called the principal planes of vibration of the ray ; 
the corresponding intensities — that is, the maximum and 
minimum of J : 

J = A+ B ± ~V(A - B'f + 4G' 2 (46) 

are called the principal intensities of radiation. 

For linearly polarized light we have, as we saw, 
0 = ± VAB. Their by (46) the maximum intensity is : 

J = A + B 
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and the minimum intensity is : 

J = 0 

For natural light A = B and (7 = 0; hence, by (46), 
J = A = B, whereas the principal planes of vibration 
become indeterminate, by (45). 

If the co-ordinate axes x and y are made to lie in the 
directions of the principal planes of vibration, the equation 

(45) is satisfied by <f> = 0 and <£ = Consequently 

C = 0 in that case, and the component of the intensity 
of radiation in any arbitrary direction <f> becomes, by (40) : 

J $ ~ A cos 2 <f> 4- B sin 2 <j> (47) 

So A and B are the principal intensities, in agreement with 

(46) . 

Since in general the intensity of radiation. J t{> depends, 
according to (40), on three constants, it follows that to 
characterize the state of polarization of the radiation it is 
not sufficient to measure the two components J 0 = •= A and 

Jn — B . Actually, for example, when A -.= B, the light 
2 

might be natural light or light linearly polarized at an 
azimuth of 45°, or it might be composed of some mixture 
of natural and linearly polarized light. To determine the 
third tensor-component fg = G , it is thus necessary to 
make a third measurement, say of the component which 

vibrates at the angle ~ : 

° 4 

J n _ 4+J? + c . . (48) 

4 

Then A, B and G are fully determined. 

§ 12. Total Reflection. Having disposed of the questions 
relating to the phenomena of reflection and refraction, 
we have still to deal with the ease which wo deferred 
provisionally at the conclusion of § 8 because the solution 
we found was inapplicable to it. This is the case where 
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sin0>^ (n< 1) because the angle of refraction 6 X then 
becomes imaginary, by (20). 

It is obvious that it is not sufficient to say here that 
since the angle of refraction is imaginary, no refraction 
can occur, and consequently all the light is reflected. 
For the equations of the problem are by no means satisfied 
by setting the wave-functions of the refracted wave equal 
to zero and those of the reflected wave equal to those of 
the incident wave. It is impossible to evade the task of 
finding expressions for the wave-functions which fulfil 
all the conditions demanded by the theory. On the other 
hand, we must bear in mind that a complex solution of 
an equation is also a solution, and that even if it has no 
real significance, it may easily serve as a guide in finding 
a real solution. 

In fact, every equation between complex quantities 
resolves into two equations between only real quantities. 
If we remember also that all the equations between the 
field-strengths of the incident, reflected and refracted 
waves — both those which hold in the interior of the body 
and also those which express the boundary conditions, 
are linear and homogeneous with respect to the field- 
strengths and have real coefficients — it follows that if all 
the equations are satisfied by certain complex values of 
the field-strengths, the real parts of these complex values 
also satisfy these equations and hence represent a real 
physical solution of the problem. Hence we immediately 
derive a real solution from every arbitrary complex 
solution merely by omitting the purely imaginary parts 
from the complex values of the field-intensities and 
retaining only the real parts. In putting this idea into 
practice we reflect that since the angle of refraction 9 X is 
complex, the coefficients n and a and the wave-functions 
of the refracted and reflected waves will also assume 
complex values. For the sake of complete generality we 
shall also assume the same for the incident wave ; that is, 
we now take / and g to stand for two arbitrary complex 
functions of a single complex argument. If we denote 
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the real part of any complex quantity by prefixing an 
R, then we get for the real field-components of the incident 
wave, from (8), the expressions : 


Ex 

Ey 


-7=b/ 

V e 


Hz 

H, 


0 


llcj 


( 49 ) 


Ez = —7=Rg 
VC 


Hz = Rf 


where the quantity t - ^ is again to l>c put as the 

argument of / and g. As hitherto, * is the wave-normal 
whereas the y-axis again lies in the incident plane and 
the z-axis is perpendicular to it. 

As for the refracted wave its wave-normal aq is complex. 
Hence we shall refer its real field-components right from 
the outset to the real co-ordinate system £, rj, £, winch is 
determined by the incident normal £ and the plane of 
incidence £ v (Fig. 1). For this purpose we shall first 
form the complex field-strengths according to the 
equations : 

E s = E Xl . cos 6i - E n . sin d 1 
E„ = E Xl . sin 6 1 + E H . cos 

and so forth, and replace the components Kr„ E Vl , 
Ez ... in it by the wave-functions j\, g x after the model 
of (8), and then replace these in turn, by (U>) and (18), by 
the wave-functions / and g of the incident wave. We then 
get the following expressions for the real field-components 
of the refracted wave : 

Ee = -Kgx/sin d x 

V«i 

JS = -4= • B/xi/cos 0 X 
v € 1 

E t = y== • BojJ 


H s = Rcrtf sill 0 X 

Hr, = - Rcr x g cos 0 X • ( 50 ) 


i/c = iW 
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Here we have to insert t — ^ everflWSt&re as the 
argument of / and g, in that : 

x x = £ cos ^ sin . 


Finally, for the reflected wave, whose normal x' is real, 
we obtain, by (8), (15) and (17), the real field-components : 

o 

II 

bS 

o 

1! 

J? 

$ 

Eyr = -L Rjl f 

Hr=-Rog (52) 

E Z ’ = ~^r=Rag 

V e 

cs 

II 

£ 


with the argument t — —• of / and g. 

§ 13. Although the preceding equations solve the 
problem of total reflection in principle, the theory 
cannot be directly applied to the general case of arbitrary 
wave-functions because it is not possible to specify the 
real part of a complex expression so long as its form is 
not known. For this reason we are compelled, if we 
wish to arrive at real expressions, first to consider par- 
ticular solutions of the problem, and then to dispose of 
the general case by appropriately combining these par- 
ticular solutions. 

In choosing the particular solutions that can be used 
in optics, we bear in mind that, to represent natural 
light, we must consider, as we saw in § 11, only such 
wave-functions as can assume each one of its values 
enormously often. Hence we choose as the simplest 
particular solution a simple periodic function by setting : 

f(t) = e ia)t = cos cot i sin cot . . . (53) 

and we treat the two waves / and g separately by first 
taking : 


g(t) = 0 . . 


• (54) 
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Then, by (49), we get for the incident wave : 

Ex = 0 Hx = 0 

= cobc(<-|) H,= 0 (55) 

E. = 0 H 2 = cos co (if — 

This is a special sine wave of “radian frequency” w 

and “ revolution ” frequency or vibration number v — “L 

Zrr 

and wave-length A = (II, § 40). 

To find the real field-components for the refracted wave 
from (50) we first calculate the complex quantities cos 9 t 
and ji v For the former we get from (20) : 

cos 8 X = ± id' 

where we set the positive real quantity : 


'sin 2 6 
n* 


(56) 


(which is not to be confused with the angle of reflection 
in (21)). The coefficient /x x is given by (24). Wo resolve 
it into its real and its imaginary parts by writing : 

_ 2 sin 8 cos 0 

^ sin 9 cos 9 -h sin 8 X cos 9 X 

— 2 si** 2 0 cos 2 # T i • 2 sin 9 cos 0 . O' . sin 0 X 
sin 2 9 cos 2 0 -b O'- sin 2 ^ x 

or, if we use the abbreviation : 


6' sin 8 1 
sin 8 cos 9 


tan (0 < 5 


7T) . 


Ml 



(57) 


The choice of sign is determined by considering the 
expressions (50). 
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These contain the quantity f{t — ‘- 1 - ) where, by (51) : 

A 

%i = ± + rj sin 0,. 

If we therefore substitute the expression (53) for/, the 

exponent contains the real member ± ^ , and, if the 

solution is to be of use, this must not become co for 
£ = oo — that is, at an infinite distance from the boundary 
surface. Consequently the loiver sign is everywhere to 
be taken and we have : 

cos 6 X = — id' . (58) 

and correspondingly : 

/x 1 = 2cos|e^ (59) 

A 


Accordingly, by (50), the field-components of the 
refracted wave become : 


3 


2 sin 0, 8 - — f 

E s = cos ^ • e cos 1 oj [t — 4 1 + 

Vei z l ?i 1 




20' 8 - m p . f ( v s in0! 

E„ = - cos - ■ e sin • a> ( i 

V<u 2 l V 01 

/J. = 0, /A = o, i/„ = o 

> ^ ^ w3 


+ 


si 


H,= 2cos| e -^cos{cu( < -^) + |} 


(60) 


(61) 


As we see, even in the case of total reflection appreci- 
able vibrations occur in the second, optically less dense, 
medium. But these vibrations have the peculiarity that 
their phase does not depend at all on the co-ordinate £ 
which is normal to the boundary plane. This causes 
the wave to advance in the second body in a direction 
parallel to the boundary plane, in the direction rj, and 
consequently does not penetrate into the interior of 
the body at all. Rather, it remains confined within a 
boundary layer whose thickness is of the order of mag- 
nitude of a single wave-length A x . The greater we 
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assume the angle of incidence to be, the less the wave 
penetrates into the second body, and the more slowly it 
advances along the boundary plane. 

For the reflected wave the expressions become much 
simpler, in that here the wave-normal : 

x' = — £ cos 6 + i] sin 9 ... . (G2) 


which enters into the argument t of the wave- 

function, is real. On the other hand, by using (58) and 
(57) we get for the coefficient /x from (23) the complex 
value : 

/x = e ? ' 5 (63) 


Hence from (52) we get for the field-components of the 
reflected wave : 


Ex* — 0, Ey' — -p= cos i. co (t -f- si- 5 Ez* ~- 

V € [ V J 

IIx = 0, Hy = 0, Hz = cos jtuf^ - ^) + sj 


0 


((’,4) 


We can, of course, assure ourselves subsequently, by 
means of a simple calculation, that the expressions (55), 
(60), (61), (64) for the incident, refracted and reflected 
waves in actual fact satisfy all the conditions in the 
interior (£ g 0) and at the boundary (£ = ■= 0) of the two 
bodies. 

It is a characteristic of the reflected wave that it has 
the same amplitude as the incident wave, but in contrast 
with ordinary reflection, a phase which is displaced by 
an angle 8 with respect to the incident wave. 

This abrupt change of phase is equal to zero, by (57), 
for the limiting angle (sin 6 = n, 9 0), where the 

total reflection coincides with ordinary reflection, and 
increases in the same sense as the angle of incidence 6 
increases until it assumes the value rt for grazing incidence 



The above results hold for the /-wave which vibrates 
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in the plane of incidence. It is obvious that the g-wave 
which vibrates perpendicularly to the incident plane obeys 
fully analogous laws, which nevertheless are distinguished 
in a characteristic way from those which have hitherto 
been obtained. Since the process of calculation is 
exactly the same it will be sufficient to summarize the 
results brieflv : 


Ex = 0 
E,j = 0 

1 , , 

Et = — .= COS co { t 

Ve 

and if we set : 


H x = 0 

/L = — cos co ( t 

IE = 0 


8' sin 9 


cos d sin 9 1 
we get from (24) and (23) : 


= tan ^ (0 < r < it) 


(65) 


( 66 ) 


o'! = 2 n cos - * e 2 , o = e iz 


(67) 


and hence, by (50), the field components of the refracted 
wave become : 


E- = - - cos 

V « 


E ; = 0 , E,, = 0 

r ~ f ft V sin 0 A t \ 

s - e cos 1 co ( t — + - 1 

2 V (7, / 2 




II S = 2 sin 9 cos ^ e »* cos j <o (t — ^ 1 

2 (V / 2j 

rr o O' T ~~a~ ■ f ft V sin dt T 

H„ = — 2 n6 cos - e am ■ w i — + - 

2 t V ?i 2 


( 68 ) 


H t = 0 

and by (52) the field-components of the reflected wave 
become : 


Ex' = 0, Ey = 0, Ex' = cos Ico (f — — ) + t! 

I V ?' J (69) 

JL' = 0, Hy — — cos | co (i — ) 4- t II 0 1 
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The wave g which vibrates perpendicularly to the 
incident plane thus experiences a different phase change ' 
r at total reflection from the wave / which vibrates in 
the incident plane; by (57) and (66) we have in general 
that § > r. It is only in the extreme cases of the limiting 

angle ($ x = and of grazing incidence (d =---- that 

§ and t coincide. 

§ 14. To solve the problem of total reflection entirely 
we have yet to generalize the particular solution that has 
been found for the case of any arbitrarily given waves 
/(f) and g(t). This can be done simply by first multiply- 
ing the expression (53) by means of any arbitrary complex 
constant, by which the vibration acquires an arbitrary 
amplitude and an arbitrary phase-constant, and then * 
performing a summation over the whole multiples" of a 
definite radian frequency w. This causes the real part 
of /(f) or g(t) respectively to assume the form of a Fourier 
series (II, § 38) : 

SC, cos (^-b ft, . (70) 

and it is known that every function which is periodic with 
respect to the time T can be represented by such a series, 
By (12) its initial term O 0 is set equal to zero. In assum- 
ing the wave-functions / and g to be periodic we are not 
actually introducing a limitation, for we can make the 
period T so great that for the measurement only those 
times t come into question which lie within a period— 
that is, between 0 and T. This occurs, for example, 
if we denote the moment of timo at which the source 
begins to emit light by t = 0, and the moment of time 
at which the last measurement has been performed by 
t = T. 
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SPECTRAL RESOLUTION. INTERFERENCE. 
POLARIZATION 

§ 15. In deriving the laws of reflection and refraction in 
the preceding chapter no assumptions were made about 
the form of the wave-functions / and g which were used. 
This also applies to the laws of total reflection. For the 
expansion of a wave-function into a Fourier series is only 
a particular form of mathematical representation. In 
particular we must guard ourselves against thinking that 
the wave-functions must have any properties of periodicity 
within the large fundamental period T. Rather they can 
even have the value zero at times, and can behave quite 
arbitrarily before and after these zero values. Hence it 
also follows that no objective meaning for the form of the 
wave-function can be ascribed to the amplitudes On of 
the Fourier series (70), nor to the fundamental period T 
on which they are dependent. 

To get information about the form of a wave it is there- 
fore necessary to make special measurements. The 
method which suggests itself most readily in principle is 
that which was used by II . Hertz to find the form of the 
electromagnetic waves which he had discovered — namely, 
to produce stationary vibrations by means of reflection 
from a perfectly conducting and hence perfectly reflecting 
surface (III, § 92). For it is in this process that the 
periodic character of a wave manifests itself, owing to the 
appearance of equidistant nodes and anti-nodes, which 
allow the wave-length, and hence the frequency, to be 
measured. The corresponding experiment with optical 
waves was first successfully performed by 0. Wiener 
(1890), who allowed the violet rays of an arc lamp to fall 

35 
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normally on a silver mirror. To make the resulting nodes 
and anti-nodes of the electric field-strength visible, a very 
thin membrane-like layer of collodion containing silver 
chloride, which is sensitive to light and does not appreciably 
disturb the vibrations of light, was placed obliquely at a 
small angle a over the mirror (Fig. 2). No photographic 
action was then found to occur at the places where the 
membrane encountered a nodal plane of the electric field- 
strength, for example at the mirror itself, whereas a 
maximum darkening of the membrane occurred at the 
places midway between the nodes, where the planes of 
the anti-nodes (shown by the broken lines in the figure) 
are intersected by the membrane. The wave-length can 
then be directly calculated from the distance between 
the dark stripes and the angle of inclination of the mom- f 
brane to the mirror. 



Fio. 2. 


On account of the smallness of the wave-lengths of 
ordinary light, it is customary to use as the unit of length 
in optics the micron, I/x = 10“ 4 cm., or I/x/x (millimicron) = 
10“ 7 cm., or, more often, the Angstrom unit, lA 10 8 
cm. The visible spectrum stretches from about 0-4 /a 
( violet) to about 0*8/x (red). 

§ 16. Although this experiment proves the periodic 
character of a light-wave, it by no means follows from this 
that the Fourier series (70) reduces for this wave to a 
single term, a single order number n. For the measure- 
ment, even in the case of the sharpest spectral lines, 
always admits of a very great range of order numbers, 
whose frequencies vary but imperceptibly. For since : 


OJ 


27 rn 

2 r 


2ttv 


(71) 


and since T is enormously great compared with the time 
occupied by a single light vibration, the order numbers n 
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that come into question are enormously great and a 
change in the order number n by one or more units effects 
no appreciable change in co. 

Thus we must conceive of light of a definite colour, so- 
called monochromatic or “ homogeneous 55 light, not as a 
wave of a single period, like the sound-wave of a definite 
tone in acoustics, but as composed of numerous waves of 
almost equal periods. There is no monochromatic light 
in the absolute sense, but only in a more or less approxim- 
ate sense. Expressed mathematically, the condition for 

homogeneous light of frequency oj 0 = is that in the 

Fourier series (70) only those amplitudes C n differ from 
zero for which : 


n 


o < < i 


(72) 


The degree of homogeneity is determined by the greatest 
value which the ratio (72) can assume without C», vanishing. 

A particularly clear picture of the behaviour of the 
wave-function of homogeneous light is obtained by 
substituting in the Fourier series (70) : 

n = n Q + (n ~ n 0 ) 

and expanding the cosine correspondingly. The series 
(70) can then be written as the single term : 

C'o-cos (^~°- + e 0 ) .... (73) 

where : 

C'o COS e ti = ZCn COS (— ': >l - ~^ 1 - +8„ 

C\ sin 9 0 = SCn Sin 


and can be regarded as a single vibration of frequency ct> Q , 
amplitude C 0 and phase-constant 6 0 . It is true that C 0 
and 0 O are not strictly constant, but on account of (72) 
they only change relatively slowly, and in general 
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irregularly, with the time t ; in fact, the more slowly the 
more homogeneous the light. 

This shows the close relationship which exists in the case 
of homogeneous light between the fluctuations of the 
amplitude and the phase on the one hand, and the degree 
of homogeneity on the other. An absolutely constant 
amplitude and an exactly regular phase could be possible 
only in the case of absolutely homogeneous light ; every 
kind of fluctuation denotes a lack of homogeneity. 

§ 17. A far sharper resolution than that given by 

is obtained by the re- 
flection or transmission 
of light through a trans- 
parent plane - 'parallel 
plate. This problem, 
too, may in general bo 
solved by means of the 
methods here developed. 

If a. plane wave falls 
on a transparent plane 
plate of thickness I) at 
an angle of incidence 6 
(Fig. 3), it is partly re- 
flected, it partly pene- 
trates into the plate ; the 
penetrated portion of 
the wave partly travels to and fro within the plate 
between its boundary faces and partly passes out through 
the rear or the front face of the plate. 

The complete solution of the problem is contained in the 
following assumption. In the first medium (air) there are 
two waves, one directed towards the plate at the angle of 
incidence 9, the other directed away from the plate at 
the angle of reflection 9. In the second medium (glass) 
there are also two waves, both making an angle 0 X with the 
normal to the plate, the one directed from the front to 
the rear boundary surface, the other moving in the 
contrary direction. In the third medium (air) there is a 


reflection at an opaque mirror 
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single wave which advances in a direction parallel to the 
original incident wave. 

These five waves are connected by four boundary 
conditions — namely, one each for the electric field- 
strength and the magnetic field-strength at each of the two 
parallel boundary surfaces. From these it is possible 
uniquely to determine the other four waves, when the 
incident wave is given. 

We achieve the same object, not so directly but more 
concretely than by the preceding method, if we take as 
our starting-point the course of the process in time ; this 
has the fundamental advantage that it can then also be 
applied when, as always happens in reality, the cross- 
section of the incident wave and the surface of the plate 
does not exceed all limits. 


Let the incident 


wave 




be 


represented 


graphically in Fig. 3 by any arbitrarily selected ray, which 
we assume strikes the first surface of the plate at A v so 
that x represents the wave-normal at the point A v The 
reflected ray (that is, the reflected wave) then passes from 
A x out into the air again with the coefficient of absorption 
or “ weakening ” ( Schwachungskoeffizient ) which is given 
by (23) ; whereas the refracted ray, which has the co- 
efficient of absorption i u 1 , gets as far as B x in the plate. 
Here it splits up again, partly to enter the air in the 
forward direction and partly to turn back in the plate 
towards A 2 . And so the same process continues until 
finally the energy of the ray is essentially used up and the 
remainder can be neglected. According to this method 
of reasoning, the four waves which were considered above 
and which are to be calculated from the incident wave 
present themselves as four sums having an infinite number 
of terms each of which is known and can be specified. 

Let us first consider the transmitted wave by choosing 
any wave-plane B' X) B' z . . . sufficiently far from 
the plate. This wave arises from the superposition of 
those waves which are represented graphically by the rays 
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AiB-yB'i, A 1 B 1 A 2 B 2 B' 2: .... If we denote the length 
A X B X = B X A 2 = A 2 B 2 = ... by S, and the lengths 
B-iB'^ B 2 B' 2 , ... by p x , p 2 , . . • then the expressions 
for the ware-functions which meet in the selected wave- 
plane and which correspond to the individual rays result 
if we substitute the light-paths quoted in the argument of 
the function and add, besides, the absorption coefficients 
fx,, fi x , jj!, which present themselves at every reflection 
and refraction. Here p' and p' x denote the coefficients 
which are valid for a wave which passes from glass into 
air. They are derived from the expressions for /a and pq 
by exchanging the angles 9 and 8 X with each other. 

Accordingly the wave-function of the first ray, which 
has undergone two refractions, at A x and at B x , runs, at 
the point B\ : 


yiM 1 


•A* 


x 

Q 


3 M 

(Ji Q J 


and the wave-function of the second ray, which has under- 
gone two refractions at A x and at B 2 and two reflections at 
j and A 2i is, at the point B' 2 : 

yiy'V'i •/(* - " - ^ - ”)» an( l so forth. 


Hence we get for the wave-function of the whole of the 
transmitted light : 

Y rf) Vi/(* ^-±- 1)S - (74) 

P = 0 x H ( ll ( I 

Here the coefficients \x! and fi\ can bo reduced in a 
simple manner, in view of (23) and (24), to terms of the 
single coefficient [jl by means of the equations : 

P ~ Mi = 1 + /x, = 1 — [l . . (75) 

Concerning the argument of f, we can, by (14), express 
the traversed lengths in terms of a common denominator : 

— x ~f~ (2 p -f- 1)5??/ pp -i- x 

i 


( 76 ) 
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and so reduce the total path of the light to the equivalent 
path in air (“ optical length of path ”). We can then 
reduce the lengths p P + 1 to terms of the first of these, p 1; 
by means of the following relation which is immediately 
evident from Fig. 3 : 

Pi - Pa = BiC ==«' 

Pi - P3 = 2 ft • • • ( 77 ) 

Pi - Pi- n = 

Finally we have the following expressions, which can be 
easily obtained for the lengths § and e which have been 
introduced : 

S = (78) 

cos 9 X 

c = B X B >> . sin 0 = 28 sin 9 X sin 9 \ . 

e = 2D tan 9 X Bin 9 j 

If wc introduce the following abbreviations for the terms 
of the argument /not involving p : 

+ = a .... (go) 

and for those multiplied by p : 

28 n — e .. n __ 


then the wave-function (74) now runs : 

. . . (82) 

§ 18. In order to be able to perform this summation 
we require to make some assumption about the form of 
/. We therefore suppose / to be expanded into a Fourier 
series and first consider a single term of this series of 
frequency o>. Since it is often more convenient to cal- 
culate with exponential functions than with trigono- 
metrical functions, we regard the real vibration as the real 
part of a complex vibration and so write : 

f(t) = e (0t . 


(83) 



42 


THEORY OF LIGHT 


CHAP. 


where, for the sake of brevity, we follow the usual custom 
and suppress the symbol R. 

The expression (82) then becomes : 

2(1 - 

p = 0 

and this sum is equal to : 


(1 - /x 2 ) e icoa 
1 


(84) 


The real part of this complex expression represents the 
required wave of the transmitted light. Since the 
variables t and p x are contained only in a, it is simply 
periodic and has the radian frequency to. By (10) its 
intensity is proportional to the square of the amplitude. 
This square can be most simply calculated by observing 
that it is at the same time the square of the absolute value 
of the complex quantity (84) ; that is, it is the product of 
(84) and its conjugate imaginary. Since we are concerned 
only with the ratio of the intensity of the transmitted 
light to the intensity of the incident light, and since the 
absolute value of (83) is equal to 1, we obtain for the ratio 
of the intensity of the transmitted light to that of the 
incident light : 


T _ (1 — n 2 ) . e i<oa (I — pfi) . e~ ia)a 

j (i - w 

1 — 2/x 2 cos col 3 + yfi 


(85) 


By (81), (78), (79) and (20) we here have : 


$ = 


2nD cos & 2D 


Thus : 


o/3 


9 

4:77 D 

X 


ffl 


cos 6 X 


■cos i 


• ( 36 ) 


where X x denotes the wave-length in the substance of the 
plate, and ; 
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Maximum intensity is attained when £ = 2rr . p, that is, 
when : 

2D 

cos ff x = p (integer) .... (88) 


and then / = 1. There are intermediate minima of 
intensity, of amount : 


J = 



(89) 


The differences between the maxima and minima 
become the greater the more /x approximates to 1 — that is, 
the greater the angle of incidence 9. 

But what characterizes the phenomenon here considered 
most sharply in its bearing on the measurements is the 
steepness of these maxima. This is due to the fact that 
the parameter £ in (8G), and consequently also p in (88), 
is in general a very great number, on account of the 
difference in the orders of magnitude of D and A x . Con- 
sequently the intensity of the light varies very markedly 
with the ray direction 9. 

Now we must bear in mind, as already emerges from the 
considerations of § 5, and as will be further explained in 
§ 36, that in optics we never observe a single definite ray 
direction or a single system of parallel wave-planes, but 
always only a cone of ray directions, which may be only 
very narrow. Thus if we ascribe to every ray direction 6 
the point of an infinitely distant screen at which all rays 
parallel to this direction aim, then there will always be a 
great number of such adjacent points 9 on which the light 
will impinge, and since J varies greatly with 9, a great 
number of maxima and minima of intensity will appear 
next to each other on the screen, the distances between 
them being regulated by the successive order numbers p 
in (88). Now if /x is nearly equal to 1, the intensity J, by 

(87) , is in general weak, and it is only when the condition 

(88) is exactly fulfilled that it becomes equal to 1 — that is, 
the maxima stand out in sharp contrast on the dark back- 
ground. 
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All these considerations refer to a perfectly definite 
frequency co of the wave — that is, to absolutely homo- 
geneous light. Since, as we have seen, real waves always 
contain many, although very adjacent, frequencies, the 
observed radiation always contains several systems of 
maxima and minima of the type considered, each system 
having the maxima and minima at different intervals ; 
and in every ray of the transmitted light or, respectively, 
at the point 6 on the infinitely distant screen several rays 
of different colours meet, and in general they are of 
different intensity. It may also happen that two different 
colours have their maximum intensity at a definite point 
d, so that by (88) : 

p x 1 = p'X\ 

or, referred to air : 

pX = p'X' ( 90 ) 

Then the maximum of order p for the wave-length A 
coincides with the maximum of order p' for the wave- 
length Ah If A and A' differ only slightly the nearest 
maxima of the two wave-lengths, whose order numbers 
p and p f differ by only a few units, will fall close together, 
and the two systems of maxima and minima will also 
nearly coincide in the immediate neighbourhood. Rut 
if we proceed to higher order numbers the little difference in 
the distances between the maximum and minimum will 
make itself felt, and the maxima of the two systems will 
move apart, so that the maxima of the greater wave- 
length, say A', will, on account of their greater separation, 
move ahead of those of the smaller wave-length A of the 
corresponding order number. Finally, when the order 
number is sufficiently increased, the maxima of A' with 
the order number p' -f r will, so to speak, catch up those 
of A with the next higher order number p + r + 1. Then 
the two systems will again coincide, but in such a way that 
the order number of the maximum of A will have increased 
by one more than the maximum of A', that is : 

{p + r + 1)A = (p f + r) A' 
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This, combined with (90) gives : 


A 7 - A 1 

A r 


(91) 


Thus by counting up the maxima between two suc- 
cessive coincidences of the two systems we have a very 
exact method for measuring small relative differences of 
wave-length. 

The intensity J' can be found for the reflected light just 
as has been done for the light transmitted from the plate. 
But it is unnecessary to carry out the calculation in detail. 
For the result comes out at once by applying the energy 
principle, which states that the sum of the intensity J' 
of the reflected light and of the intensity J of the trans- 
mitted light is equal to the intensity 1 of the incident 
light. Thus by (87) : 


J' = 1 — J = 


4/x 2 sin 2 1 


1 — 2/x 2 cos £ -f /x 4 


(92) 


The reflected light is, as we say, “ complementary 5 ’ 
to the transmitted light. The minima of intensity are 
zero and stand out sharply from the surroundings when 
/x is nearly equal to 1. In other respects, the same laws 
hold here as for transmitted light. 

Hitherto we have considered only the light which 
vibrates in the incident plane — namely, the /-wave. But 
it is evident that the same inferences may be drawn also 
for the light which vibrates perpendicularly to the incident 
plane, namely the £/-wave. The only difference in the 
results is that in the case of the (/-wave the coefficient 
a (§ 9) takes the place of the coefficient /x. The position 
of the maxima and minima remains precisely the same 
as in the case of the /-wave. 

§ 19. According to the equation (86) the angle £ which 
is characteristic of the intensity of the light depends not 
only on the ray-direction 6, but also on the thickness of 
plate D. This circumstance is made use of in interfero- 
meters of variable thickness (air plate) in order to measure 
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the dependence of the intensity of light J on the thickness 
of layer D for the case where the light is incident normally 
= 0). Then by (88) those thicknesses of layer for 
which the transmitted light has the maximum intensity 
1 are whole multiples of the half wave-length in the 
substance of the layer : 

D^p.r} ( 93 ) 

This is easy to understand if we reflect that when the 
thickness of layer is a multiple of the half wave-length 
the phases of the waves which are superposed on each 

other after several reflections 
within the layer before passing 
through it differ by 2tt , so 
that all the waves reinforce 
each other. 

A phenomenon which ex- 
hibits the effect of several 
thicknesses of layer simul- 
taneously is that of Newton's 
Rings. If the curved surface 
of a plano-convex glass lens 
is placed in close contact on 
the surface of a plane -parallel glass plate (Fig. 4) a thin 
air -film is formed, whose thickness increases from the 
value zero to the value D at the distance p from 0. If 
we now allow homogeneous light to fall in a parallel 
beam normally on to the plate, we see in the transmitted 
beam a system of light and dark rings with a bright 
spot in the centre, and in the reflected light we see the 
complementary pattern. 

An approximate theory of this phenomenon is obtained 
very simply from our preceding remarks. If we restrict 
ourselves to values of p which are small compared with 
the radius R of the sphere, we may write : 

D = AB = PO = it! — VR* - r 
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Consequently, by (93), we obtain for the maxima of 
intensity in the transmitted light or, respectively, for the 
minima in the reflected light, if X denotes the wave-length 

in air : 

P = VpXB (94) 


Thus the radii of the dark rings in the reflected light 
are in the ratio of the square roots of the series of natural 
numbers. 

For two colours which differ only slightly from each 
other, such as the two components of the yellow line 
of sodium (5890 A and 589G A) the ring systems differ 
slightly and move apart more and more as the distance 
from the central point 0 increases, so that they blur 
each other's effects until for, say, p = r the maxima 
and minima, respectively, coincide again, as can be 
seen from the renewed sharpness of the again common 
ring-system. This allows us to calculate, by (91), 
the relative difference of wave-length of the two sodium 
lines. 

§ 20. In the preceding section we have been led td 
consider the combined action of several trains of waves 
superposed on each other and moving in the same direc- 
tion ; we shall now treat the same question of the super- 
position of light waves systematically for the general 
case. 


Let f ± {t — and f 2 [t — be the wave -functions of 


two rays which vibrate in the same direction and which 
propagate themselves in the same direction. They com- 
bine to form a single ray whose wave-function is : 


/=/x+/ 2 (95) 

It follows from this that the intensity J of the resulting 
rays depends on the intensities J x and J 2 of the two com- 
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ponents in the following way, if we disregard a factor 
which is of no importance : 

J~J 1 -\-J 2 J r 2/ lt /2 (96) 

Thus the intensity of the composite ray deviates in 
general from the sum of the intensities of the individual 
rays. This phenomenon is called the “ interference 55 
of the two rays 1 and 2. Under certain circumstances 
no interference occurs — namely, when the wave-functions 
f 1 and f 2 are fully independent of one another — that is, 
when the rays 1 and 2 are “ non-coherent. 55 For then, 
as in (43), we have : 

/ 1/2 ~ J i -A ~ 9. 

This is always the case, for example, when the two rays 
come from different sources. Hence two such rays, even 
if they have exactly the same colour, never interfere 
with each other. But non-coherence is not a necessary 
.condition for the absence of interference. Rays can even 
be “ fully coherent 55 and yet not interfere : we take 
“fully coherent 55 to mean that one of the two wave- 
functions is completely determined by the other. This 
occurs, for example, in the case of homogeneous light, 
if the arguments of j\ and f 2 always differ by an odd 

multiple of Then, in s]>ite of complete coherence, 
£ 

/ 1/2 = G The general case is exemplified by cc partial 
coherence/ 5 where and / 2 are not independent of 
each other, and yet are not completely defined by each 
other. 

§ 21. A somewhat different way of regarding the same 
problem results if we write the two wave-functions f x 
and f 2 as Fourier series having the same fundamental 
period T. Then the resulting wave-function / — A -f/ 2 
may also be represented directly as a Fourier series by 
compounding each two partial waves of the two series, 
which have the same frequency, into a single partial 
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wave. We then have, for example, for the radian 
frequency co of the partial waves : 


a x cos | co (t — -) + 9j\ 
a 2 cos jco (t — + 9 2 j 


. (97) 


In summing these two expressions an important part 
is played by the phase-difference of the two waves : 


e 1 - e 2 = a ( 98 ) 

Instead of the phase-angle 9, sometimes the constant 
length d or the time constant 8 is introduced into the 
wave -function by writing : 


»cos(2^-i+i)} 

< *° 08 M t T 5 -x)} 


and we then speak of the Cf difference of path, 5 ’ d x — c? 2 , 
of the two weaves or of the “ retardation,” 8 X — S 2 , of the 
second wave with respect to the first. The phase differ- 
ence 2tt corresponds to a difference of path of one wave- 
length, A, or a retardation of one period of vibration, r. 
By adding the two expressions (97) we get for the resulting 
partial wave : 

a cos (t — + flj (99) 

where : 

a cos 6 = a x cos d x + a 2 cos # 2 
a sin, 8 = a x sin 8 X + 

Thus the amplitude and the phase-angle of the resultant 
wave are formed from those of its components according 
to the same rules as the absolute value and the direction 
of a vector which results from two other vectors. 

From (100) and (98) we obtain for the square of the 
amplitude a : 

a 2 = a x 2 + a 2 2 + 2 a x a 2 cos A . . . (101) 

E 


( 100 ) 
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This equation expresses that two absolutely homo- 
geneous waves of the same frequency moving in the 
same direction and with their vibrations in the same 
direction interfere in general. This does not contradict 
the law of § 20 that two rays coming from different 
sources of light and having the same colour never interfere. 
For, according to § 16, there is no absolutely homogeneous 
light in nature ; rather, even the most near homogeneous 
optical ray always contains very many partial waves of 
nearly the same vibration frequency ; each of these partial 
waves interferes with the partial waves of the same 
frequency in the other ray. Now if the two rays are 
non-coherent the phase-difference A changes in a quite 
irregular manner in passing from one pair of waves to 
another, so that no appreciable interference effect can 
occur. 

§ 22. We shall next consider the composition of two 
waves which propagate themselves in the same direction, 
but whose directions of vibration are at right angles to 
each other, as in § 1 1 except that now we start from the 
Fourier expansion. As in § 11 let z be the direction of 
propagation, x the direction of the electric field-strength, 
y that of the second wave. We can then represent the 
electric field-strengths of two partial waves of the same 
frequency co after the model of (97) by : 


a cos |o> [t — + 0 X j ~ x 
b cos | co (t — + 0 2 j = y 


. (102 


where we choose appropriate units for x and y. Without 
affecting the generality of our problem wo may take 
a^b^O. This method of representation enables us to 
form a very clear picture of the law, according to which 
the resultant electric field -strength of a definite wave- 
plane z varies in magnitude and direction with the time, 
since it is represented at every moment by the position 
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of the point x, y which moves about in the plane z = 
const. 

As for the orbit of this point, it is obtained by eliminating 
the time t from the equations (102), which is done most 


simply by calculating the values of cos 



sin a )(t — -) 9 then squaring and adding them and putting 


their sum equal to 1. The equation of the path then 
comes out as : 


^ , t _ 

a 2 ^ b 2 


2 xy 
ab 


cos A 


sin 2 A . 


(103) 


This is an ellipse which degenerates into a straight line 
when A = nrr ; whereas when A = (n + l)rr its axes coin- 
cide with the co-ordinate axes. In general the angle <j> , 
which one axis of the ellipse forms with the rr-axis, is 
determined by the relation : 

, . 2a&cosA / lru \ 

tan2 ^-^r w ' ' • (104) 

which can be obtained directly from (45) if we insert the 
following values, which correspond to our present case : 

„ 7 ~ ab . 

G ^jg = -w cos A. 

If, further, we inquire into the velocity with which the 
elliptic orbit is traversed, we get from (102) by differ- 
entiation : 

xdy — ydx = ojab sin, A . dt . . . (105) 

This signifies that the motion of the reference -point 
x, y occurs in accordance with the law of sectional areas 
(I, § 50), and moreover, since ab > 0, it occurs in the 
positive or negative sense according to the sign of sin A. 
Since the 2 -axis of our right-handed co-ordinate system 
(I, § 16) is the direction of propagation of the light, the 
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sense in which the rotation represented by the motion 
occurs has a definite physical meaning. 

In optics light is called dextrorotatory (or dextro- 
gyrous) if the motion of the reference -point x, y is clock’ 
wise to an observer looking towards the light. Por 
this reason we have to regard the sense of rotation of 
dextrorotatory light as negative, and that of lsevo- 
rotatory light as positive. 

We can obtain a clear conception of the way in which 
the motion of the reference-point depends on the con- 
stants a, b, A by keeping, say, a and b constant, but 
allowing the phase-difference A to increase from 0 to 2 tt. 
When A = 0 the elliptic motion degenerates into a recti- 
linear vibration, which by (104) or (102) is inclined to 
the x-axis and lies in the first and third quadrants, the 
angle of inclination being : 

6 = tan -1 - < ^ • • • (105a) 

and the amplitude Va 2 -j- b 2 . 

Now if the phase-difference A increases, the ellipse widens 
and is traversed in the positive sense. The light is then 
lsevorotatory. At the same time the angle <j> becomes 
smaller — that is, the major axis of the ellipse moves 

towards the £-axis. When A = ~ the ellipse attains 

its greatest curvature, the semi-axes are a and b , and 
their directions coincide with those of the co-ordinate 

axes. If A increases beyond the major axis passes i 

over into the second and fourth quadrants and the ellipse 
at the same time becomes flatter and flatter until, when 
A = 77, it again shrinks together to a straight line which 
is inclined at the angle given by (105a) to the x-axis. From 
then on, when A increases still further, the form of the 
ellipse is exactly repeated in the reverse direction, as 
we see from the fact that the orbital equation (103) 
remains unchanged if A is replaced by 2tt - A. The 
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ellipse again widens and its major axis again moves 

3 

towards the a;-axis, which it reaches when A = But 

Zi 

the orbit is now being continually traversed in the negative 
direction. Finally, when A = 2 t t, the original motion is 
again resumed. An equally instructive example is that 
obtained by considering the phase-difference A fixed and 
the amplitudes a and b of the components of vibration to 
change gradually. Among other tilings we see from (104) 
that then the direction <j> of the axes depends only on the 
ratio of the amplitudes a and 6, and that the same holds 
for the ratio of the lengths of the axes. 

If we enquire into the condition which makes the 
vibration circular, we obtain the answer most simply 
from the expression (104) for the direction of the axes, 
which must become indeterminate in the case of a circle 

— that is, must assume the form jj. This gives us the 

following two conditions for a circular vibration : 

a = b, A = (n + . . . (106) 

which must be fulfilled simultaneously. The circle is 
traversed in the positive or negative sense according as 
n is even or odd. 

§ 23. The elliptic vibration considered in the preceding 
section at the same time represents the most general case of 
the composition of any arbitrary number of plane waves 
moving in the ^-direction and having the same radian 
frequency to. For every linear vibration may be resolved 
into two linear components which vibrate in the x- and 
the ^-direction and have the phase-difference zero, and 
all the vibrations which occur in the x~ or the y-direction 
yield on being superposed, by § 21, a single vibration in 
the direction in question. 

We can also picture the general case of the elliptic 
vibration, in another manner besides that involving two 
linear vibrations — namely, by means of two circular 
vibrations with opposite senses of rotation. This is 
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obvious, if we reflect that every linear vibration may 
be regarded as the resultant of two circular vibrations 
having the same radius but being performed in opposite 
directions. For the two revolving points then always 
meet at the same points, and these points define the 
direction and half the amplitude of the resultant linear 
vibration. 

In this way we obtain for the general case your circular 
vibrations; of these there are two pairs of vibrations 
which occur in the same direction and so compound in 
each case into a single circular vibration in the same 
direction : the two single resultants vibrating in opposite 
directions can then be compounded. The sense of the 
final resultant circular vibration is, of course, that of 
the single resultant which has the greater radius. 

§ 24. Let us finally glance at the way in which the 
electric field-strengths x and y in (102) depend on the 
wave-normal z. If t is constant, these equations repre- 
sent the so-called “ wave-line, 55 in the present case an 
elliptical helix : and by (105), if we replace the dt in it 
dz 

by — — > we get a right-handed screw (II, § 32) if 

sin A < 0 — that is, if the sense of rotation is negative. 
The pitch of the screw is then equal to the wave-length A. 

If we take into consideration what was said in § 23 about 
the sense of traverse of the dextrorotatory light, it follows 
that dextrorotatory light is represented by a right- 
handed screw, and laevorotatory light by a left-handed 
screw. The motion of the reference-point x, y in the 
plane z is obtained if we displace the helix in the direction 2 
with the velocity q and fix our attention on its instantaneous 
point of intersection with the plane z = const. If we 
reverse the direction of displacement of the helix, a wave 
results which propagates itself in the negative direction 
of the 2 -axis and at the same time turns in the reverse 
sense, so that a right-handed screw again represents 
dextrorotatory light and a left-handed screw lsevo- 
rotatory light. 
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§ 25. We shall next enquire wliat definite phase- 
differences A can be realized between the two vibration 
.components x and y. If wo start with a linearly polarized 
wave having its direction of vibration in the first and 
third quadrants, then in (102) the ratio of the amplitudes 
is constant : 

~ = tan<£<l (107) 


just as in (105a) ; and the phase-difference is 9 X — d 2 = 
A = 0. We shall cause this wave to be reflected at 
some arbitrary angle of incidence, so that the electric 
field-strength x vibrates in the plane of incidence. For 
ordinary reflection (sin 6 < n) no abrupt change of phase 
occurs, so the phase-difference A remains equal to zero, 
and the two reflected waves again combine to form a 
linearly polarized wave, except that, by (105a), the 
direction of vibration becomes changed owing to the 
difference in the reflection coefficients /x and cr. This has 
already been remarked in § 11. 

But it is quite different in the case where the wave is 
totally reflected (sin#>n). For here, conversely, the 
amplitudes a and b of the two wave-components remain 
unchanged, whereas a phase-difference A = S — r now 
introduces itself, the value of which is obtained from 
(57) and (GO) : 




9 ' . cos 9 
sin 9 . sin 9-, 


(108) 


and S' is determined by (56). 

This phase-difference A is, as we have already found at 
the end of § 13, always positive and vanishes at the two 
extreme cases of the limiting angle — namely, when 0=0 

and when 9 = - (grazing incidence). Thus between these 

Ji " 

values it attains a maximum. This occurs, by (108), 
at the angle of incidence : 




( 109 ) 
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and the maximum value is : 

A = tt — 4 tan -1 M .... (no) 
0 < tan _1 w < 

The smaller the index of refraction n, the greater the 
phase-difference A that can be obtained. After the 
reflection the two vibration components combine into an 
elliptic vibration, the direction of whose axes, <f>', is by 
(104) and (107), given by : 

tan 2</>' = tan, 2<f> cos A . . . (m) 

To obtain circularly vibrating light we must, by (106), 
produce a phase-difference of If we enquire what is 
the least value that n can have in order that the phase 
difference ^ may be obtained, we find that it follow's from 
(110) that : 

n = tan| = V2-l. . . . (112) 

The reciprocal of this value is \/2 + 1 = 2-414. A 
substance must therefore have an index of refraction of 
at least this value if it is to be possible for the phase- 

difference | to be produced in it between the f- and the 

(/-waves by total reflection. Of the known substances 
only diamond fulfils this condition. 

Now if a phase-difference of | cannot be practically 

obtained by means of a single reflection, wo can achieve 
this fairly conveniently by means of two successive 
reflections, as by a Fresnel rhomb. This is an oblique- 
angled parallelepiped of transparent glass (Fig. 5) whose 
angle is chosen so that a ray of light which passes through 
the glass and impinges on the air at the angle of incidence 
d = a. undergoes total reflection with a phase-difference 

A = | between the two vibration components. By (108) 
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t his occurs for any definite sort of glass at two angles, 
either of which may he chosen. 

If we make the ray of light which was assumed at the 
beginning of this section to vibrate in the direction cf> 
with respect to the ce-axis, fall normally from below on the 
base of the surface at A (Fig. 5), the equation (107) will 
hold for the ratio of the amplitudes before and after the 
entrance of the ray into the parallelepiped. The first 
reflection at 0 on the oblique bounding surface then occurs 

at the angle of incidence a, with a phase-difference of - 

and unaltered amplitudes; the second reflection occurs 
at P on the opposite bounding surface 

with a further phase-difference of ~ 

4: 

and again unaltered amplitudes ; finally, 
at B the ray emerges normally into 
the air with uniformly diminished 
amplitudes and unchanged phase- 

difference ~ + j — * 5 . Hence the ray 

which emerges at the top at B per- 
forms ellipotic vibrations, whose princi- 
pal axes coincide with the co-ordinate 
axes x (towards the right) and y 
(towards the back of the figure). 

Since the light which is incident at A is assumed to 

be vibrating in the first and third quadrant [ ^ 

the light which emerges at B is lmvorotatory, by (105). 
But if the light vibrates in the second and fourth quadrants 

— ^ <f> < tt^ the emergent light is dextrorotatory. When 

77 377 

(j) = - or — we have tan cf> — ± l and a = b, and then the 

light executes left-handed or right-handed circular 
vibrations. 

§ 26. Our discussion so far has been entirely restricted 
to a definite radian frequency co — that is, to a single 


z 
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partial vibration in the Fourier series of the optical 
vibration. Now since even the most nearly homogeneous 
light wave in nature includes a great number of partial 
vibrations, we have in all cases to sum up over all these 
partial vibrations and so arrive at the relevant laws for 
the resultant vibration. 

Instead of this, however, we may also, as in (73), 
regard every homogeneous optical vibration as a single 
elliptic vibration of the corresponding frequency with 
the lengths and directions of the axes changing slowly 
and irregularly, or as consisting of two circular vibrations 
in opposite directions whose radii change slowly and 
irregularly. 

But since the laws which have been deduced in the last 
sections, from § 22 onwards, are quite independent of the 
particular value of the frequency co, they also apply 
unaltered when many partial vibrations of almost the 
same frequency are superposed. This holds, in particular, 
for the ratio of the amplitudes (107) of the two vibration 
components in a partial wave of linearly polarized light, 
as well as for the phase -difference A between the two 
vibration components, which is given by (108), and hence 
also for the direction (111) and the ratio of the lengths of 
the two principal axes of the resultant elliptic vibration 
(cf . end of § 22). All these quantities are the same for the 
different partial waves — that is, they are subject to no 
fluctuations. Optical vibrations whose direction of axes 
and ratio of axes are constant are said to be “ elliptically 
polarized.” This does not only express, then, that the 
optical vibrations are elliptic — that is obvious — but also 
that the vibration ellipse, even if it changes slowly and 
irregularly, nevertheless retains the directions of its axes, 
the ratio of its axes, and its sense of traverse unchanged 
throughout. In the same way “ circularly polarized 55 
light is represented by circular vibrations whose radius 
alters slowly and irregularly. Elliptically or circularly 
polarized light may be produced from linearly polarized 
light by means of Fresnels rhomb, as above described. 
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The direction of the axes is determined, as we saw, by the 
plane of incidence at total reflection, and the ratio of the 
axes by the azimuth of polarization. Prom this we 
easily see that the two vibration components of light of 
this kind are completely coherent (§ 20). 

§ 27. Taking these results as a basis we are now in a 
position to analyse up to a certain point any plane wave 
that we may encounter. Let us take any light-wave 
present in nature, whose origin is unknown to us, and 
undertake the problem of determining its form by making 
certain measurements, that is, to specify as far as possible 
the properties of the two wave-functions f and g of which 
it is composed. 

First we resolve the light spectrally by one of the 
methods of interferometry above described. We are 
then able to deal only with, homogeneous light. 

We then determine, by § 11, the two principal directions 
of vibration and the two principal intensities for a homo- 
geneous light-wave of this kind. Three cases can occur. 
The first is that one of the principal intensities vanishes 
entirely. This means that the light is linearly polarized. 
The second is the opposite extreme where the two principal 
intensities are equal to each other. The light is then either 
circularly polarized or is natural light or a mixture of both 
(partially circularly polarized light). To distinguish 
between these three possibilities we allow the light to 
fall normally on a Fresnel rhomb (Fig. 5). If the light is 
circularly polarized, the amplitudes of the two vibration 
components / and g are equal to each other and their 

Stt 

phase-difference then amounts to or according as 

the light is kevorotatory or dextrorotatory. After 
emergence from the glass the amplitudes are still equal 
to each other but the phase-difference now amounts to tt 
or 27 r, respectively — that is, the emergent light is linearly 
polarized and vibrates in the second and fourth or first 

and third quadrants, respectively, at the angle <f> = or 
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to the £-axis. Thus here the double reflection supple- 
4 

ments the original phase-difference of the two components 
to form a whole multiple of n — that is, the Fresnel rhomb 
acts as a “ compensator.” 

On the other hand, if we are dealing with natural light, 
no compensation occurs, and the emergent light, like the 
original light, shows no trace of differences in the different 
directions of vibration. 

If, lastly, the light is a mixture of circularly polarized 
and natural light, the natural part remains unchanged, 
whereas the circular part again yields linearly polarized light 

whose azimuth is ^ or respectively, and whose in- 
tensity is represented by the difference in the principal 
intensities of the emergent light. 

In the third and most general case in which the two 
principal intensities differ from zero and from each other, 
we allow the light to fall on the rhomb in such a way that 
the direction of the greater principal intensity is coincident 
with the r-axis (Fig. 5). If the light is elliptically polarized 
the two vibration components have the phase-difference 

~ or -jr, respectively, and the ratio of the amplitudes is 

A A 

equal to Owing to the action of the compensator, 

Cl 

linearly polarized light is produced, which vibrates in the 

1) TT 

direction <b = =F tan" 1 - and forms an angle < . with the 

a;-axis. But if the light is a mixture of linearly polarized 
and natural light, the greater principal axis also coincides 
with the £-axis in the case of the emergent light. If, 
finally, the light is partially elliptically polarized the greater 

principal axis in the emergent light forms an angle 

with the :r-axis, which indicates the ratio of the axes of 
the vibration ellipse, while the intensity of the polarized 
portion is obtained from the difference in the principal 
intensities. 
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In this way we obtain in every case that presents itself 
a certain insight into the conditions of polarization of a 
given ray. There is still, of course, a considerable re- 
mainder which cannot be defined and which is caused by 
the physical complexity of the elementary processes that 
occur in every light source ; it receives expression in the 
numerous terms of the Fourier expansion. 



CHAPTER III 
GEOMETRICAL OPTICS 

§ 28 . The laws of reflection and refraction which we 
have developed above owe their simple character essen- 
tially to our assumption that both the wave-points of the 
light (§ 5) and the boundary planes of the active body are 
infinite planes. In reality these surfaces are neither 
unlimited nor plane, and hence, strictly speaking, the 
above simple laws are not applicable to nature at all. 
Nevertheless, in practice they represent an extra- 
ordinarily close approximation to reality, when these 
assumptions are very nearly fulfilled, valid, that is, when 
both the cross-sections of the wave -fronts and the 
boundary planes, as well as the radii of curvature of these 
surfaces, are very great in comparison with the wave- 
length of the light under consideration. Now since in 
optics the wave-lengths involved are, in general, of a 
smaller order of magnitude than the dimensions of the 
bodies used, the optical laws of propagation, reflection and 
refraction assume a particularly simple form in contrast 
with acoustics where, it is true, the laws of wave-motion 
also hold, but where this assumption is not in general 
fulfilled. 

Geometrical optics comprises the relationships which 
are obtained when the wave -fronts of the light and also 
the surfaces of bodies are imagined to be divided up into 
many small parts, and the laws of propagation, reflection 
and refraction of plane waves of unlimited extent which 
fall on unlimited plane surfaces are applied to these 
individual small parts. Each of these small portions of 
wave-front, which are regarded as plane, is graphically 
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represented by the corresponding ray which is normal to it, 
andeveryraypropagatesitselfrectilinearly with the velocity 
characteristic for that body until it reaches another body 
where it is reflected and refracted as at a plane surface in 
accordance with the relative position of the normal. 
Hence geometric optics is also called “ray optics” in 
contrast with “wave optics,” which is more general but 
also more complicated. 

§ 29. We take as our first example the course of a ray 
of light through a prism — that is, through a transparent 
body which is bounded by two plane faces which are 
inclined to each other at an angle </>. After what has been 
said above, we must exclude processes which occur in the 
immediate neighbourhood of the edge A (which in Fig. 6 
we take perpendicular to the plane 
of the diagram), because the radius 
of the refracting surface becomes 
infinitely small there. 

A ray which is incident at B at 
an angle 9 in the plane of the 
diagram is refracted along BG and 
emerges into the air again at C at 
the angle 8 r to the normal at G. The following relations 
hold for these two refractions : 

sin 0 = n sin 8 X .... (113) 

sin O' = n sin, 8\ . . . . (114) 

Moreover, the sum of the three angles of the triangle 
ABC gives : 



that is : 

e 1 + 9\ =--</> (ii5) 

From these equations we can find how the angle of 
emergence 6' depends on the angle of incidence 0. 

The total deviation S of the ray caused by the prism — 
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that is, the angle between the emergent and the originally 
incident ray — is given by : 

s = (0 - e x ) + ( e ' - e\) 

or, in view of (115), by : 

§ = 0 + O' - cf> . . . . ( 116 ) 

If we investigate the manner in which the total deflec- 
tion 8 depends on the original angle of incidence 9 , we 
get, by using (113), (114) and (115) : 

clS _ _ cos 6 . cos 9\ ' ' ^ /-j^m 

d9 cos 9 X cos 9 ' 

This expression vanishes if 9\ == 9 ±i which simul- 
taneously makes 9' = 9 — that is, when the ray traverses 
the prism symmetrically — so that the lengths AB and AG 
are equal to one another. In this case the deflection S 0 is 
a maximum or a minimum, and we find out which it is 
by finding the value of the second differential coefficient. 
This is obtained by differentiating (117) with respect to 9 ; 
for the value 9 ' = 9 , which we are considering, it has the 
value : 

(d?d\ 2 sin (9 + fl x ) . sin (9 - 9 ± ) / llg x 

\d9 2 ) o ~ sin 9 cos 8 cos 2 9 L 

If the index of refraction of the substance of the prism 
n > 1, then in the symmetrical position 9' = 9 the second 
differential coefficient is positive and hence the deflection 
is a minimum. 

This limiting position can be determined fairly sharply 

experimentally by leaving the position of the incident ray 

unchanged, but rotating the prism to and fro on some axis 

parallel to its edge and observing the direction of the 

emergent ray. If the minimum deviation S 0 has been 

found by trial, the refractive index of the substance of the 

prism is calculated from (116), (113) and (115) and comes 

out as : f 

. + S 0 

sin, 0 

71 = .... (119) 
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§ 30. Before passing on to consider curved wave-fronts 
and curved surfaces of bodies we must make some 
remarks of a general character. The primary wave-front 
in an isotropic body is a spherical surface. Eor every 
point-like source of light sends out spherical waves of 
light which can be represented at a sufficient distance 
from the source by a beam of divergent rays. If these 
rays are refracted or reflected at the curved surface of 
another body, their directions will, in general, no longer 
intersect at a point — that is, they will no longer be homo- 
centric but astigmatic. If, in a special case, they are 
again homocentric, then the point at which they meet is 
called the “ optical image ” of the source of light, being 
real or virtual according as the rays converge towards or 



diverge from the point of intersection. Since light can 
always traverse its path in the reverse direction (end 
of § 10), the source and the image can exchange places, 
and are therefore also called “ conjugate ” points. The 
image of an infinitely distant point of light, or, what 
amounts to the same thing, the point of intersection of 
rays which result through refraction and reflection from 
an incident beam of parallel rays, is called a “ focus ” — 
that is, when it exists. If, conversely, we make the focus 
a source of light, we obtain a beam of parallel rays. It 
is on this circumstance that we base the experimental 
production of plane waves of light. 

§31. We shall now investigate the refraction of light 
at the spherical bounding surface between two bodies 
1 and 2 (Pig. 7). Let % and n 2 be the refractive indices 
of the two bodies, C the centre of the curved surface, and 
r = PC', the radius of the sphere. We first consider the 
image of a point of light A x situated on the straight line 
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PC , the axis of the system. Of the rays which emerge 
from A x only the axial ray A X P passes through the 
spherical surface without being deflected. Let another 
ray which emerges from A x at an angle 9 X to the axis 
meet the spherical surface at Q. It becomes refracted 
there according to the law : 

sin A X QC __ n 2 
sin A 2 QC n x 


and meets the axis at the point A 2 at an angle d 2 . If we 
calculate the two sines from the triangles A X QC and A 2 QC 
we get the relation : 


A 1 C 'A 2 C __ n 2 
A X Q A 2 Q u x 


( 120 ) 


In general the position of the point A 2 will change with 
the value of 6 X or with the point Q at which the ray 
impinges — that is, the beam which starts from A x becomes 
astigmatic after the refraction. But if 6 X is taken 
sufficiently small, we can replace the lengths A X Q and 
A 2 Q in (120), except for errors of the second order, by the 
distances A X P and A 2 P. We then obtain, if we denote 
the distances of the points A x and A 2 from P by e 1 and e 2 : 


or : 


e x 4- r . e 2 — r _ n 2 
h * e 2 ~ n x 


( 121 ) 


Va _j_ ^2 ^2 

% ^ H ~ r 


( 122 ) 


Since the point Q no longer occurs in the latter ex- 
pressions, the beam of rays which starts out from A x is 
homocentric even after the refraction, provided the 
conical angle of the beam is small ; and the point A x now 
has an optical image, the conjugate point A 2 . By (122) 
the law for the formation of the image is expressed by a 
linear equation between the reciprocals of the distances 
of the conjugate points from the refracting spherical 
surface. 
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If A x moves off to infinity on the left, A, moves to the 
focus F 2 , whose abscissa is : 


J? 

h ~n,-n, (123) 

If A x moves to the right, A s first moves off to infinity 
on the right— namely, until A x reaches the focus F v whose 
abscissa is : 

f 1 n» — w, r (124) 


If A x moves beyond F x , the refracted ray becomes 
divergent— that is, A., appears on the left as a virtual 
image, which follows the light A x as it moves further to 
the right, until it catches up A x at the point P. The 
point P is its own conjugate, just like the point 6, since 
the rays which start from C are not deflected at all. 

The abscissae f x and /, of the foci exhibit the simple 
relationships : 


/a “A = r 
fi _ % 

/a n 2 


(125) 

(126) 


If we substitute for r and n from them, 
image formation (122) becomes simply : 


the law of 


(127) 


All these lelations, as we shall see, are capable of being 
considerably generalized. We may add here that the 
formation of an image by reflection at the spherical 
surface may be included in the above by writing 
n x = — n 2 . 

If the rays which start out from A x intersect at the 
image point A 2 after refraction, the angular aperture or 
the “ divergence ” of a beam at A x is different from that 
of the conjugate beam at A 2 . Por from Pig. 7 we get 
for the angle of inclination of the two conjugate rays 
and A 2 Q to the central axis : 


sin : sin 0 2 = A 2 Q : A X Q 
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and to a corresponding degree of approximation : 

sin 9 1 : sin 6 2 = A 2 P : A X P = e 2 : e x . , (128) 

§ 32. We shall now investigate the image of a point of 
light which does not lie on the central axis. Then B x has 
an optical image, just as much as A x , namely the con- 
jugate-point B 2 which lies on the straight line B X C . The 
circle which results from rotating B x about the central 
axis A X A 2 forms an image which is the circle described 
by the rotation of B 2 , and the circular area of radius A X B X 
forms its image on a certain surface of rotation which is 
bounded by the circle B 2 and which has its centre at A 2 . 
If we restrict ourselves to distances which are not far 
from the central axis, we may regard this surface of 
rotation as a plane of circular contour (radius A 2 B 2 ) 
which is perpendicular to the central axis. The surface 
of light then forms a “ similar image on the image plane, 
and the following equation holds for the ratio of the 
length of the light-path A X B X = l x to the length of the 
image-line A 2 B 2 = l 2 : 

l x \l 2 = A X G :A 2 G 

= (e 1 + r) : (e 2 -r) 

and by (121) : 


Together with (128) this gives the following simple 
relationship between the ratio of the linear dimensions 
and the ratio of the ray directions in two conjugate planes : 

sin 6 X = n 2 l 2 sin d 2 (129) 

If the point A x coincides with the focus F lt l 2 = 00 and 
6 2 = 0. The two planes which pass through the two 
foci and are perpendicular to the axis are called the 
“ focal planes.” 

§ 33. All the above laws can easily be generalized for 
the case of any arbitrary number of refracting spherical 
surfaces. We only need to choose the measures for the 
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abscissae of two conjugate points A x and A 2 appropriately. 
Hitherto we have used for the abscissae their distances 
e x and e 2 from the refracting spherical surface, and we 
have reckoned the distances e x as positive towards the 
left and the distances e 2 as positive towards the right. 
The origin P is conjugate to itself. We shall now choose 
a special origin for each of the two abscissae, but in 
such a way that the two origins O x and 0 2 are again 
conjugate to each other, but otherwise arbitrary. Further, 
we shall assume the directions of the abscissae x x and x 2 
to be positive towards the right in both cases. If a x and 
a 2 then denote the new origins 0 1 and O z of the spherical 
surface, where by (127) : 

& +4 = i (130) 

a x a 2 v ' 

then the relations between the new abscissae x x and x 2 and 
the old abscissae e x and e 2 are (Fig. 7) : * 

O x A x = O x P — A X P = a x — e x = x x 
0 2*d 2 == ~P A- 2 — -P 0 2 == — Cl> 2 === X 2 

and we get for the law of image formation from (127), if 
we replace the e’s in it by the x’s, and take (130) into 
account : 

a . l . - /i + = i . . . (131) 

X ^ %2 

This is again a linear equation between the reciprocals of 
the abscissas of the two conjugate points. 

§ 34. We now pass on to investigate a system con- 
sisting of an arbitrary number of line-centred spherical 
surfaces , that is, spherical surfaces whose centres all lie on 
a straight line, the central axis of the system. If we 
again restrict our attention to such rays and such points 
of light as lie near the central axis, that is, if we consider 
only small values of 6 and l , the laws of image formation 
may be derived directly from those obtained above. Let 
us first take a point of light A situated on the axis and a 



70 


THEORY OF LIGHT 


CHAP. 


narrow beam of light which starts out from it in the first 
body. Corresponding to it there is a definite conjugate 
point on the axis and a definite conjugate beam of rays 
(real or virtual) in the second body ; and corresponding 
to these again there are definite conjugate quantities in 
the third body. Proceeding in this way we finally arrive 
at a definite conjugate point A' and a definite conjugate 
beam of rays in the last body. 

Now since the reciprocals of the abscissae of two succes- 
sive conjugate points, when referred to any two points as 
origins, depend linearly on each other, there must also 
be a linear relationship between the abscissae x and x' of 
the conjugate points A and A' in the first and last bodies, 
and the law of image-formation has a form similar to 
(131), namely : 

f f 

•h+S = 1 (132) 

For x = 0, x' becomes equal to 0. A and A' then 
coincide with the origins 0 and O'. The constants / and/' 
are clearly the abscissae of the foci in the first and the 
last bodies. If the foci lie on different sides of the origins, 
/and/' have opposite signs. 

By (129) the following relation holds for the lengths 
l and V of two conjugate lines that pass through A and A' 
and are normal to the central axis : 

nl sin 6 = n'V sin 8' (133) 

where n and n r denote the refractive indices, 8 and 6' the 
angles of inclination which any two conjugate rays that 
pass through A and A ' make with the central axis. 

In the sequel we shall reckon both the lengths l and V 
and also the angles 8 and 9' as positive in the same 
direction : this causes no change of sign in (133). 

§ 34$. Let us now see whether there are two conjugate 
points for which l = l\ that is, for which the representa- 
tion of the two surfaces that pass through them normally 
to the central axis is fully congruent. We shall call these 
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two points H and H' and shall try to find their abscissae 
x and x' . From (133) we then get : 

n sin 6 = n' sin 9' 


or, if we draw any two conjugate rays through H and H' 
which make the (small) angles 9 and S' with the central axis 
and meet the planes that pass through the origins 0 and O' 
at the distances l 0 and V 0 from the origins (Fig. 8), we have : 


n- 


y 

x 


From this and (132) we get : 


(134) 



O' * ^ 


Fig. 8. 


The points II and H' which are uniquely defined by 
these equations are called the “ principal points 5 5 and the 
planes which pass through them and are congruent 
representatives of each other are called the C£ principal 
planes 5? of the system. 

The results become particularly simple if we make the 
origins 0 and O' coincident with the principal points H 
and H'. For then x and x' become equal to zero in (134), 
and Z 0 = l' Q , consequently : 

/;/' = - TKtt' .... (135) 

that is, the foci lie on opposite sides within or without the 
principal points, and their distances from the principal 
points, the so-called ec focal distances/’ are in the ratio 
of the refractive indices. We then get for the law of 
images (132) : 

n n r __ n 
x x' ~~ / 


■ ( 136 ) 
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If the first body is of the same nature as the last, for 
example, air, then n r = n and the focal distances become 
equal, so that the law of images becomes simplified to : 

3. _ 1 1 

s x ' ~ f 

§ 35. Analogous to the two principal points H and H ' 
whose planes form congruent images there are two con- 
jugate points K and the nodal points, whose beams 
form congruent images, so that for each pair of conjugate 
rays 9 = 0'. Their position is obtained directly by 
observing that two such conjugate rays which pass 
through K and K' and are parallel to each other intersect 
the principal planes at the same distance l 0 = 1' 0 from the 
central axis; that is, the abscissae x and x' of the two 



Fig. 9. 



nodal points are also equal to each other. Hence it 
follows from (132) that : 

* = *'=/+/'=/( 1-| (13V) 

that is, the nodal points are at the same distance from the 
principal points, namely a distance equal to the difference 
of the focal lengths. 

The relative position of the principal points, foci and 
nodal points are depicted graphically in Fig. 9, in which 
the abscissae of the nodes are positive and the foci lie out- 
side the principal points, so that f +f> 0, /< 0, f > 0. 
Then : 

EK = H'K' = H'F' - FIF 

Consequently : 

H'F' = HK + HF = FK 
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that is, the point-pairs FF\ HE' and KIP , which, taken 
together, are also called the “ cardinal points/' have a 
common centre, the so-called optical centre of the system. 
If there is only a single refracting surface (Fig. 7), this 
.surface forms a congruent image on itself ; that is, the 
two principal planes coincide in it. Consequently the 
nodal points also coincide, namely with the other self- 
conjugate point, the centre G of the sphere, as is imme- 
diately clear, since every ray which passes through C is 
self-conjugate. This causes the equations (135) and (137) 
for the ratio and the difference of the focal distances to 
transform directly into the special equations (126) and 
(125) which were found earlier in § 31. 

If we have a system of glass lenses in air, the focal 
distances are equal to each other and the nodal points 
coincide with the principal points. If only a single thin 
lens is present, the principal points coincide at a point 
within the lens. 

In general, the position of the cardinal points leads 
directly to a simple geometrical method of determining 
the point P' in the last body which is conjugate to any 
point P in the first body (Fig. 9). We need only draw 
the conjugate rays to any two rays that pass through P. 
The simplest ray to use is PK which passes through the 
nodal point K . The parallel ray P'K' is conjugate to it. 
Besides this, we can also take the ray through P parallel 
to the axis. This ray has as its conjugate that ray in the 
last body which passes through the focus F' in it and 
intersects the principal plane in it at the same distance 
from the central axis as that at which the original ray 
intersected the principal plane in the first body. This 
serves to define Pb In addition, although this is un- 
necessary, we may also consider the ray PF. This ray 
has as its conjugate a ray, parallel to the axis, which 
intersects the last principal plane at the same distance 
from the axis as the first ray does the principal plane. 

§ 36. A complete description of the laws of optical 
image formation must deal not only with the course of 
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the rays but also with the intensity of the rays. We 
must above all bear in mind here that a finite quantity of 
radiant energy never starts out from a point but always 
from a surface, and never in a definite direction but 
always within a cone of directions. We can depict this 
circumstance graphically by representing the energy of a 
beam of rays by means of the number of rays contained 
in the beam. Every beam then contains a quadruply 
infinite number of rays, since the surface from which it 
starts out defines a doubly infinite number of points, and 
each point a cone comprising a doubly infinite number of 
directions. But there is a certain difference which must 
not be overlooked. The rays which start out in various 
directions from a luminous point are coherent with one 
another, corresponding to the circumstance that the point 
is the centre of a wave-surface, every point of which is in 
the same phase. But the rays which come from two 
different points of the initial surface are in general non- 
coherent, particularly when the initial surface acts as 
the source of light. Hence in general no interference 
occurs between such rays ; their intensities simply become 
added. But it may happen — and this is sometimes over- 
looked — that two different points of a surface emit 
coherent light; and again, conversely, non-coherent rays 
may start out from a single point of the surface. For 
example, when the surface (aperture or slit) is illuminated 
by light from another source ; in this case certain inter- 
ference phenomena can occur for which geometrical optics 
taken alone cannot account. We shall therefore always 
assume that the initial plane is self-luminous. 

If the source of light is placed in a focal plane, a beam 
of parallel rays results, but, in conformity with the above 
remarks, not in the sense that all rays of the beam are 
parallel, but rather that, corresponding to every indi- 
vidual point of the luminous surface there is a particular 
direction of the rays of the beam, whereas the size of the 
luminous surface determines the angle of aperture of the 
cone of rays which emerge in various directions. Con- 
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versely, rays which come from a practically infinitely dis- 
tant source, such as the sun, will not all become focused 
at one point ; rather, a small picture of the sun will be 
produced in the focal plane, since all the parallel rays 
which come from a definite point of the sun meet at a 
definite point in the focal plane. 

The intensity J of the radiation emitted normally by a 
small surface element / perpendicular to the central axis, 
the rays lying within a narrow cone of angular aperture 
Q is, after what has been said above : 

J = K .f.Q (138) 

where the finite quantity K is called the “ specific inten- 
sity of radiation 55 or the “ specific luminosity ” of the 
beam. 

In the same way the following relation holds for the 
conjugate beam which lies within the conjugate cone of 
rays Q ' and is incident on the con j ugate surface /' : 

J' = 1C ./' .Q' .... (139) 

By “the laws of image formation there is a general 
relationship between the quantities / and Q in the first 
and the last body, which is fixed by (133). Since, on 
account of the similarity of the image with the object, 
the form of the surface / is of no account, we may choose 
it to be a circle of radius l. We then have 

f = l'hr. 

If, further, we imagine Q to be a circular cone with a 
small angle of aperture 0, then. 

Q = 477 sin 2 ~ = 0 2 7t. 

In the same way : 

Q' = d'h r. 

These values, combined with (138) and (139), give the 
relation : 

n JQ = n j'Q' .... (140) 
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v?J »'*J' 

K ~ K' 


( 141 ) 


If the total intensity of radiation J is the same in the 
two conjugate beams — which will, however, never actually 
be the case, if only because the radiation suffers a loss at 
every refraction — we have : 

|-5° r Kf- K 'i' (142) 


That is, the intensities are in the inverse ratio of the 
squares of the velocities of transmission of the rays. 
For n = n r we get K = K'. Hence the intensity of a 
beam of rays can never be increased by optical image 
formation alone. It is true that the radiation is concen- 
trated on a smaller surface at the focus but then it 
diverges in various directions. 

§ 37. The simple laws above developed for optical image 
formation by means of a system of line-centred spherical 
surfaces mostly hold only as a first approximation to the 
real conditions. For in practical optics we deal both 
with surfaces of finite size as well as with finite cones of 
rays. In addition to the corrections necessitated by these 
circumstances there is another due to the fact that the 
refractive index depends on the colour of the light (§ 9). 
If all these influences are to be taken into account, 
problems of a difficult kind arise; their solution has 
played the predominant part in the highly developed 
technique of construction of optical instruments. But 
the laws of geometrical optics, such as we have described 
them in § 28, are not affected by these practical problems. 
They may all be condensed into a single simple graphical 
law, which is one of the first examples of the equivalence 
of an integral law with a differential law (cf. Ill, § 42), 
namely, Fermat's Principle of Least Time . 

This principle controls the course of a ray of light 
through any arbitrary number of different bodies with 
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any boundary surfaces by asserting that the time which 
light requires to pass from any definite point on a ray to 
another definite point on the same ray is less along the 
actual path of the ray than along any neighbouring 
path. 

A simple proof of this principle is obtained if we con- 
sider the optical wave-surface that starts out from any 
luminous point P as centre. According to § 5 there is a 
definite wave-front corresponding to every moment of 
time t. It signifies the boundary to which light which is 
emitted by P at the time t = 0 has advanced in all 
directions during the interval of time 0 to t, independently 
of the presence of any bodies whatsoever in the vicinity 
of P. Thus the wave-front comprises the end-points of 
all the rays which start out from P in all directions, each 
propagating itself along its own particular path. Now if 
Q is any point on the wave-surface, PQ denotes a definite 
ray, and the light takes the time t to propagate itself 
from P along the path of the ray to Q. Now if the light 
from P were to take a smaller time than t to pass from 
P along any neighbouring path to Q, then the propaga- 
tion of light from P in all directions would proceed 
beyond Q in the interval 0 to t, and this contradicts the 
definition of wave-surface. Consequently the time t is 
the minimum of all the times taken by the light to arrive 
from P to Q along any path. 

The whole of geometrical optics can be developed from 
Fermat’s principle, and even for non-homogeneous and 
anisotropic bodies. The result that a ray travels in a 
straight line in a homogeneous isotropic medium follows 
directly, according to Fermat’s jDrinciple, from the cir- 
cumstance that the straight line is the shortest distance 
between the two points P and Q. 

Let us consider, as a second example, the passage of 
light through the plane boundary surface of two homo- 
geneous isotropic bodies. If q and q x are the velocities 
of propagation, A and B the projections of P and Q on 
the boundary, 0 the point through which the ray PQ 
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passes through the boundary (Fig. 10), and 6 and 0 X the 
angles of incidence and refraction, then the time which 
light takes to pass from P to Q is : 


where : 


. PC CQ r r x 

f= T + ir“s + j; 

r 2 = AP 2 + AC 2 , and r x 2 = BQ 2 + BC 2 


(143) 

(144) 


If the point C through which the ray crosses the 
boundary is displaced, the distance AP, BQ and' AB 
remains constant, whereas AC = a and BC = b change, 
so that : 



So; + 8b = 0 


(145) 


Consequently, by (143) and 
(144) : 

f . a§a , b8b 

dt = t - 

rq r x q x 

If we set U = 0 and take (145) 
into account, it follows that : 

ah ■ a ■ 0 
q : q x = - : — = sm 0 : sin d x 
r r x 


which is identical with (20). 

Fermat’s principle may be used in a similar way to 
find the direction of propagation of light through any 
arbitrary bodies with any arbitrary surfaces. 

If the end-point Q of a ray which starts out from a 
luminous point P coincides with a point which is optically 
conjugate to P the ray-path PQ is not definite ; rather, 
there are an infinite number of rays from P to Q. Con- 
sequently the time which the light takes to pass from P 
to Q is the same along all these paths, and the optical 
wave-surface around P has a singular point in Q. 



CHAPTER IV 
DIFFRACTION 

§ 38. If the wave-fronts or the surfaces of the bodies 
may no longer be regarded as of infinite extent compared 
with the wave-lengths, the laws of geometrical optics 
fail, the conception of a ray and its propagation lose 
their meaning, and phenomena occur which are denoted 
by the term “ diffraction.” It is clear from the discus- 
sion in § 28 that the part played by diffraction becomes 
so much the greater the greater the. wave-length of the 
light used. An exact theory of diffraction is possible 
only on the basis of, firstly, the differential equations for 
electromagnetic waves in the medium under considera- 
tion, which we shall always take to be a vacuum space 
in which light has the velocity q = c, and, secondly, the 
boundary conditions which hold at the surfaces of the 
diffracting bodies. 

The differential equations for the interior can always 
be reduced to a single differential equation (III, § 87), 
which must be satisfied by every component <f> of every 
vector. It is the wave-equation : 

= c 2 . A</> (146) 

The boundary conditions differ according to the nature 
of the diffracting body. They come out most simply for 
so-called absolute conductors (III, § 92), in which the 
electric intensity of field vanishes and which act as 
perfect mirrors. The opposite extreme is given by the 
so-called black bodies which absorb all the incident light. 
Diffraction problems are among the most difficult encoun- 
tered in optics. The exact solution of a problem of this 
kind was first given by Arnold Sommerfeld in 1895. 

79 
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Since the wave-function <j> can always be represented 
as a Fourier series, the general solution of (146) can be 
reduced to a sum of particular solutions, each of which 
corresponds to a definite frequency co. We may write as 
a particular solution of this kind : 

4> = ^ (146a) 

where ifs denotes a certain complex function of the space- 
co-ordinates, whose absolute value (modulus) character- 
izes the amplitude and whose argument is the phase- 
constant of the periodic motion (cf. § 18 above). 

The wave-equation (146) then leads to the following 
differential equation of the function iff in space-co-ordi- 
nates : 

A^ + ^V = 0 

or : 

Ai A + ^ = 0 • • • • (1466) 

which is to be integrated with due regard to the boundary 
conditions. 

Now as the rigorous solutions of the diffraction problem 
formulated in this way- are of a very intricate character 
and demand a relatively heavy mathematical equipment, 
it is the more important that approximate solutions, 
which are fully sufficient for the short optical waves, 
should be obtained by introducing Huygen’s Principle. 
This principle is founded on the idea that every point of 
space on which a light-wave impinges itself becomes the 
centre of a new light-wave which spreads out from it in 
concentric spherical surfaces. G. Kirchhoff was the 
first to succeed in formulating Huygen’s principle exactly ; 
it forms an extension of Green’s theorem in the theory of 
potential, which we proceed to deduce here. 

§ 39. Let U and V be two uniform functions of the 
space-co-ordinates x, y, z. These functions, as well as 
their differential coefficients, are to be continuous. Then 
for a space whose volume-element is dr, surface- 
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element da and inwardly directed normal v we have, 
by II, § 18, equation (80) : 


dUdV , dUdV ,dUdV\ 


Consequently : 

■ l«- 


V 


dr — — Jv-da-J VAUdr 
- j U da — J UAVdr. 

dU ^da = j (VAU - UAV)dr . ( 147 ) 


dv 


If we set TJ = where r denotes the distance (taken 

as positive) of the point x, y , 2 : from any point A which 
lies outside the region of integration, then V fulfils the 
conditions of uniformity and continuity; further, we 
have Laplace’s equation I (129) A TJ = 0. Consequently 
(147) becomes : 


If the point A lies inside the assumed space, this equa- 
tion may also be applied, provided that we exclude a 
very small sphere with its centre at A from the integra- 
tion. This does not appreciably alter the right-hand side 
(cf. I, § 33), but the left-hand side acquires an additional 
term in the form of an integral over the surface-elements 
da 1 of the sphere wdiose normal, directed towards the 
interior of the integration space, coincides with r, thus : 


IdV 
r dr 


+ F- 


X 

Zj> 


If the radius of the sphere is taken sufficiently small, 

dV 

this integral, since — and V are finite, reduces to : 




do 1 
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where F 0 denotes the value of F for r = 0, that is, at the 
point A, and when this is added to the left-hand side of 
the equation (148) we arrive at Green’s theorem : 



which enables us to obtain the expression of a function 
F at every point A of the assumed space, if the value of 
AF is known at every point inside and the values of F 


3F 

and are known at every point on its surface. 

If the function F satisfies Laplace’s equation AF = 0 
the space integral drops out, and we require only the 
data for the surface. It is particularly in this form that 
Green’s theorem is often useful, but its importance is 


essentially restricted in that the values of F and may 

not be chosen independently of each other. For we know 
that the expression for F is everywhere completely deter- 
mined by the values of F at the surface (III, § 19) and 

r) V 

simultaneously by the values of — at the surface (II, 


§ 71). Hence if either only V or only ~ is given at the 

surface, Green’s law can be successfully applied to deter- 
mining V in the interior only when we have in some 
special way acquired a knowledge of the missing surface 
values. 

§ 40. Huygen’s principle follows directly from Green’s 
law if we substitute for F a certain wave-function which 
is chosen in such a way that the space-integral in (149) 
can be reduced to a surface-integral. This is not, of 
course, made possible by simply substituting for V a 
solution <f> of the wave-equation (146). We shall see, on 
the other hand, that it is always possible to reduce the 
space -integral to a surface -integral if, in any function (j> 
of %, y, z, t which satisfies the wave-equation (146), we 
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T 

substitute t — - for t — where r, as before, denotes the 
o 

distance of the point x, y, z from the point A, which is 
situated in the interior of the integration space — and 
then identify the expression so obtained with V in (149). 

We must therefore he careful to note that V now 
depends in two ways on the co-ordinates x, y, z: in the 
first place, explicitly, since cj> is a space-function, and, 
secondly, implicitly, since the argument : 


t - ‘- = t . . . . (150) 


contains the distance r. To call attention to this we set : 

V = cj> r {x, y, z) (151) 

in which we write the space co-ordinates in brackets on 
the line but attach the time-co-ordinate as a suffix. There 
is no reason to fear confusion between the time-argument 
r and the space-element dr, particularly as the latter will 
soon vanish from the formulae. Hence we get from (149) : 

4^(0) = - / ' v 1 1 ch + j a7 ~ 7^7 da • ( 152) 

This is the value of the wave -function </> at the time t 
at any point A of the integration space enclosed by the 
surface a which has the inwardly directed normal v. To 
reduce the space-integral to a surface -integral we first 
calculate the value of A </> T , which may not, of course, be 
set equal to (A <j>) t . For this we have : 


dfc _ (d£\ (d£\ m dr (ty\ _ im\ dr 
dx \ dx/ T \di/ r dx \dx) r c\dt) x dx 


and, further : 

dx 2 \dxV r c \dxdt) x dx c 2 \ dt 2 / r \dxs c\dt/ z dx 2 


Corresponding expressions hold for 


d^r 

dif 


d*c/> 
dz 2 



84 


THEORY OF LIGHT 


CHAP. 


By addition we get : 


= (A (f)) t 


2 / d 2 cj> 3 r d 2 cf) dr 

c \dxdt dx + dydt dy 




2 (Q) 

cr \ dt / . 


If we now substitute : 



from the wave-equation (146), the expression for A<£ 7 
may be written as a differential coefficient with respect 
to the time : 


2 9 fdcf) dr d(f> dr dcf) dr ^ 1 d<j> (f> 

~ c dt\dx dx dydy dzdz c dt r 


or : 


A^ = 


2 3 / 3 ^ 
cdt\dr 



Using this expression we get for the space-integral in 
(152) : 

[Ml dT= _2d fdr^ M 

Jr cdtJr\dr r / 

and if we express the space-element dr in terms of polar 
co-ordinates with the origin A and the angle of aperture 
da) of an elementary cone whose vertex is at A : 

dr = r 2 drdco 



- Ut! drdw -( r w + * 


(153) 


where we have now to insert in the definite integral con- 
taining the square bracket the distance r of the surface- 
element da from A as the upper limit and r = 0 as the 
lower limit. Since the latter value vanishes, we are left 
with only the term referring to the surface, and this is 
related to the size of the surface-element da, cut out 
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from the integration space by the elementary cone da, as 
follows : 

r 2 da > — — da . cos (r, v) = — da ■ • (154) 

ov 

since the directions of increasing r and increasing v form 
an obtuse angle with each other. Hence it follows that : 





<fh dr 
r dv 


<-> 


If the elementary cone daj in question cuts the surface 
of the integration space more than once — whichcan happen 
only an odd number of times, the integral (-153) in r 
resolves into several separate parts, and we can then 
easily convince ourselves that the relation (155) remains 
valid so long as we perform the integration with respect 
to dcr over the whole of the surface which encloses the 
integration space. 

If, finally, we remember that in (152) : 


dcj > T __ __ 1 /dj>\ dr 

dv \dv / T c \ dt ) x dv 


(156) 


Green’s theorem (152) passes over into Huygen’s principle : 

4r7r<f>t(0) = j Q v dcr .... (157) 

where, in accordance with (155) and (156), we have used 
the abbreviation : 

*1 

Q v ( 15S ) 

\ dv r dv crdt dvj x 
Instead of this we sometimes find written : 




!(£•)- m 

dv \ r / r \ov 


(159) 


which is, however, wrong in this form. For this does 
not express the fact that in the first term the differentia- 
tion of <f> r (x, y, z) with respect to v must be performed 
only with respect to the functional dependence on the 
argument r and not on that of the co-ordinates x, y, z. 
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The importance of Huygen’s principle is due to the 
fact that, according to it, the wave-function <f> at any 
point A of space at any time t is completely determined 

by the values of <f>, and ~ at the surface of any space 

which contains A and in which the wave-equation (146) 
holds everywhere, these values being for certain moments 
of time r which are different for every surface-element da 
and which can easily be pictured owing to the circum- 
stance that they precede the moment t by the time which 
the light takes to travel from da to A. 

Accordingly, every surface-element da continually sends 
out, in consonance with the phenomena which occur in 
it, certain spherical waves with the velocity of light in 

Q all directions, and the value of the 
wave-function at A is made up of the 
sum of the contributions which arrive 
there at the time t from all directions. 
2 We must not forget, however, that 
this theorem is valid only if the 
fig ii surface a is completely closed. There 

is no definite physical meaning in 
speaking of the effect of individual surface-elements da. 
For the integral (157) can be resolved into its differentials 
in an infinite number of ways. 

The significance of the restriction just stated becomes 
particularly clear in the case where, among the straight 
lines that start from A, there are some which intersect 
the surface a in more than one point, say in the three 
points P, Q , R (Fig. 11). The length PQ does not then 
belong to the integration space and so the wave-equation 
(146) need not hold between P and Q ; indeed, any bodies 
opaque to light may be situated there. Nevertheless 
Huygen’s principle holds in the form (157) ; and we have 
to perform the surface integration over the whole of the 
surface a and may not omit the surface -elements at Q 
and R on the plea that the light from them cannot by 
travelling in a straight line reach A . 
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§ 41. If we divide the whole integration space into two 
parts 1 and 2, so that the part 1 contains the point A, 
we obtain two spaces with the surfaces do v da 2 and dcr 0 ; 
we shall assume the latter two to denote the surface- 
elements of the surface of separation. We can then apply 
Huygen’s principle in two different forms ; firstly to the 
original total surface, with the surface-elements dcr 1 and 
d<j%, so that : 

4tt</>/(0) = J Q Vl d<j x + J Q Vi <d<y 2 

and secondly to the surface of the space 1, with the 
surface -elements da x and d< t 0 , so that : 

47t^(0) = J Q Vl dcr t + J Q, ol dcr 0 

where v 01 denotes the normal of dcr Q , directed inwardly 
towards 1. 

By subtracting the last equation from the preceding 
one we get : 

0 = j £2 v „dcr 2 J Q,, 0l dcr q 

or, since v Ql and v Q2 are opposite : 

0 = J Q Vl dcr 2 + J Q v da 0 . . . (160) 

Now, the elements do 2 and dcr 0 form the closed surface 
of space 2. Hence we have the general law that Huygen’s 
integral, when applied to a point A outside the integra- 
tion space, has the value zero — which corresponds exactly 
with Green’s theorem in the form (148). 

§ 42. We shall now apply Huygen’s principle to a 
space which is bounded on the one hand by a surface 
which comprises all the bodies that are involved in the 
process in question, namely, sources of light, screens and 
so forth, on the other hand by an immensely great 
spherical surface of radius R and with its centre in finite 
space. The equation (157) is then also applicable in this 
case provided only that the equation (146) holds in the 
whole integration space. But Huygen’s integral then 
resolves into two separate partial integrals : one which 
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stretches over the surface in finite space, having as its 
element da and as its normal v, directed towards the 
integration space, that is, towards the point A ; the 
other stretching over the spherical surface at a distance R 
and having as its element dS. We cannot make this 
second integral vanish by simply choosing R very great. 
For even if the integrand Q decreases to an indefinite 
extent as R increases, the magnitude of the surface S 
increases beyond all limits. Nevertheless we can always 
arrange that this integral can be neglected by making 
an appropriate assumption about the form of the wave- 
function <f > — an assumption which must be always admis- 
sible and in no way restrict the actual course of the 
process. Namely, we need only assume that for all times 
t < T, where T denotes an immensely long period of time, 
the wave-function (j> with all its derivatives everywhere 
vanishes; in other words, that all the sources of light 
started into activity at a far distant moment of time. 
For then the quantity and its derivatives vanish for 
all surface-elements dS, so long as : 

r =*- -< - T 
c 

that is, so long as we choose : 

R > (T + t) . c. 

It is true that the variable time t is included in this 
condition. But that introduces no difficulty. For we 
can choose R from the outset so that the inequality is 
also satisfied for even the greatest possible time that can 
occur in the measurements. In physical language, R 
must be chosen so great that an effect which starts out 
from the distant spherical surface cannot make itself felt 
at A at the time t. 

After these remarks we may also regard the equation 
(157) as valid for the case where the point A lies outside 
the surface o and where the bodies producing the effects 
all lie inside the surface. 
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Summarising, we may say that Huygen’s principle holds 
in the form (157) for every case where the surface a 
completely shuts off the point A from all the active 
bodies, so that there is no path from any of the bodies to 
A which does not somewhere pierce the surface cr. 

§ 43. We shall formulate the general problem of 
diffraction which concerns us here as follows. A given 
source of light is placed opposite a screen which is opaque 
to light and which has one or more apertures of definite 
shapes at definite points. We enquire what is the 
intensity of light at any point A behind the screen. 

We can simplify the problem in the first place by 
restricting our attention to a single point C of the light- 
source. For the intensities of light which are due to the 
remaining points have merely to be added to the intensity 
which is due to G , in accordance with the results of § 36. 
We can imagine the wave-function which starts out from 
the point of light C as a Fourier series and restrict the 
whole calculation to a single term of the series. For in 
the case of natural light of constant intensity the indi- 
vidual terms of the series are independent of one another 
and lead to no appreciable interference effects. 

To be able to apply Huygen’s principle to the presentcase 
we place the Huygen plane in such a way that it completely 
shuts off both the source of light C and also the screen 
from the “point of diffraction ” A, namely closely at the 
surface of the screen and in front of the opening p, on 
the side facing the point A . According to (157) and 
(158) we can then calculate the value of <f> at the point 
A by integrating over the Huygen surface. 

But, in order to be able to proceed, we now have to 
make a certain sacrifice in accuracy. The conditions 
are exactly the same as those encountered in applying 
Green’s theorem (end of § 39). For, to introduce the 
quantity Q in (158) into the calculation, we must know 
the values of </> and all its derivatives with respect to v 
and t at all points of the Huygen surface; and this 
assumes that the problem of calculating <f> is already 
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partially solved. This difficulty is accentuated by the 

fact that the values of <f>> ^ and at a definite point 

of the surface are not independent of one another, and that 
in general we cannot avoid arriving at an inherent con- 


tradiction if we introduce any approximate values into 


the expression for Q. Of course, it is always possible to 
establish afterwards whether such a contradiction presents 
itself and what its nature is, by making the diffraction 


point A move into the Huygen plane and comparing the 
values calculated for cf> and its derivatives from (157) with 
the values which were originally inserted. 


In spite of this we shall here be able to proceed along 
the lines proposed. For innumerable applications have 


shown that Huygen’s principle, not only when formulated 


rigorously but also when certain only approximately 
correct values of Q are used, leads to results which are 


fully sufficient for practical optics. These approximate 
values of Q are simply obtained by assuming the value 
of the expression for <f> and its derivatives at all apertures 
of the screen to be the same as it would be if the screen 


were not present at all, and by assuming the value zero 
at all other points of the surface a. The integral (157) 
is then only to be taken over the apertures of the screen, 
and in (158) we can ignore the screen altogether. 

Of course, (157) resolves into just as many component 
integrals as there are apertures, and each component 
integral refers to a particular aperture. Moreover, it is 
easy to see that the value of the integral depends only 
on the contour of the aperture. For there is nothing to 
prevent the Huygen surface o* from assuming different 
positions which have the same contour. Hence besides 


speaking of a “ diffracting aperture 55 we may also 
speak of a “ diffracting contour or edge. 55 

§ 44. In view of the above we are now concerned only 
with calculating the wave -function (j> and its derivatives 
for all points of an aperture, without paying attention 
to the screen ; a simple integration of the wave -equation 
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(146) serves for this purpose. We use as the simplest 
form for a particular solution of this equation, which 
corresponds to an absolutely homogeneous wave which 
starts out from the point-source of light G, by II (230), 
the expression : 


4 > = 


■(*-?) 


(161) 


where r ± denotes the distance of the reference point from 
the source 8, to distinguish it from r, the distance of the 
diffraction point from the surface-element da. A is any 
complex constant, the other terms have the same meaning 
as before. Cf., for example, (83). 

If we substitute the expression for cf> given by (161) in 
(158) and perform the various differentiations, the quantity 
Q v resolves into a number of terms which are of quite 
different orders of magnitude. For, since in the case of 

2ttc 

optical processes the wave-length A = — is small com- 
pared with the distances r and r 1 of the diffraction point 
A and the point-source G from the diffracting aperture, 
all those terms in Q v which contain co as a factor are 
large compared with the rest, so that the latter can be 
neglected (cf. Ill, § 88). 

Finally, if we write r = t in place of t 3 we get : 

c 



(162) 


We shall simplify the further treatment of this problem 
by introducing a specialization which does not affect the 
characteristic features of diffraction phenomena. We 
shall first assume that the diffraction screen is plane, 
and secondly, that the light is incident normally, in 
parallel rays and wave-planes; that is, that r ± has a 
constant value at all points of the diffracting apertures, 
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this value being great compared with the other dimensions 


involved. Then 


dr x 


1, whereas, since the direction of 


increasing v forms an obtuse angle with that of increasing 

r, we may set equal to — cos 8, where 6 denotes the 

“ diffraction angle.” 

Hence (162) becomes : 


*(0) = — e 


iA i'<u (< - O [d(T 


cos 2 H-®0*e A 
A r 2 


£(0 )AA e ia ^ ) (G-iS) 
where we have used the abbreviations : 


C = -A 


1 [da 


— COS 2 - COS — r- 
r 2 A 


.[ da 0 9 , 2ttT 
- / — cos 2 k sm — r- 
i Jr 2 A 


The real part of (163) represents a periodic vibration 
in t im e, a measure of whose intensity may be obtained, 
according to § 18, by multiplying the complex quantity 
(163) with its conjugate imaginary. Since, however, we 
are not concerned with the absolute value of the intensity 
of the light, but only with the ratio of the diffracted 
light to the incident light, we divide the quantity so 
found by the corresponding quantity for the incident 
light, namely by the product of the complex expression 

\A \ 2 

(161) and its conjugate imaginary : i 

r i 

In this way we arrive at the following expression for 
the diffracted light at the diffraction point A : 

J = C 2 + £ 2 (165) 

This reduces the diffraction problem to calculating the 
two integrals C and S } in which da denotes any surface- 
element of the diffracting aperture, r its distance from 
the diffraction point A , and 9 the angle of diffraction, 
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§ 45. The most characteristic peculiarity in the expres- 
sions for G and. 8 is the different manner in which the 
quantities contained in them depend on the integration 
variables. For whereas 6 and r change comparatively 
slowly with da, the trigonometrical functions, on account 

T 

of the great value of exhibit a very rapid fluctuation 

to and fro between a maximum and a minimum, with 
the sole exception of those surface-elements da for which 
the distance r from the reference point is a minimum, 
because then r changes only slightly with da. This is 
the case for all those surface-elements da which lie in 
the direct neighbourhood of the straight line which 
connects the point of light C with the diffraction point 
A and which is perpendicular to the plane of the screen. 

Hence this straight line AC, which we may call the 
“ central line,” and in particular its point of intersection 
P with the Huy gen plane, plays a characteristic part 
in determining the value of the intensity of the light 
at A. 

But there is a special case for which the central line 
does not exist at all, namely, when the point A is situated 
at infinity, as well as C. Then AO denotes no particular 
straight line, but only a definite direction; no point 
of the Huygen plane plays a particular part and the 
integrals G and 8 change their character entirely. This 
is the case to which the diff raction phenomena discovered 
by and named after Fraunhofer belong. 

§ 46. We shall first deal with the general case, the 
so-called Fresnel diffraction phenomena , starting with the 
question of the diffraction by a plane screen which is 
bounded by a single straight edge. The answer to this 
question is also of interest for the general case of any 
arbitrary edge, that is, for the general problem of shadow 
boundaries, because an edge of any shape whatsoever, 
provided that its curvature does not compare with the 
order of magnitude of the light-waves, may be imagined 
to consist solely of parts of straight lines. 
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We take as our plane of representation in Fig. 12 the 
plane through the diffraction point A and perpendicular 
to the edge of the screen ; further, the line of intersection 
of the figure with the screen is the :r-axis, the origin 
0 in the edge of the screen, a plane through the edge 
of the screen and perpendicular to the plane of the 
figure is the y-axis, and the 2-axis lies in the plane of 
the figure as shown. We designate the co-ordinates of 
the diffraction point A by x, y — 0, 2 ; that of a point in 
the plane of the screen by £, 77, £ = 0. If the aperture 
of the screen lies on the positive side of the #-axis, the 
, integration is to be per- 

formed over all the surface- 
elements dcr = d£, dr) of 
the xy - plane with positive 
values of £, that is, in the 
„ >Jt case of | from 0 to 00 and 
in the case of 77 from — co 
to + co . 

To calculate the integrals ( 164 ) we refer back to the 
above remarks on their nature and imagine the central 
line AG to be drawn, which cuts the screen-plane in 
the point P. We must now resolve the surface-integrals 
( 164 ), which are to be performed over £ and 77, into two 
parts, namely firstly an integral over those points g, rj 
which lie in the neighbourhood of the point P, so that 
their distance from P is small compared with AP — 2, 
and secondly, an integral over all the remaining points 
of the region of integration. We shall reserve a more 
definite delimitation of the first component integral till 
later. After what has been said above we can replace 
the slowly variable factors r and 6 in this integral by 
2 and 0 respectively, and so obtain : 


A , 4 

! A A, 

1 

j 

i 

i ..... 


0 P 
Fig. 12. 


P, 


0-i//<!£J,co S ?Hr + . , . | 

s -kll ded ’>* in? T + • ■ ■ j 


(166) 
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where the second component integral is denoted by the 
dots To perform the integration first for y we write : 


or, since : 


where : 


(I - x) 2 + {7f- y) 2 + a - *) 2 

y — 0 and £ = 0 : 

r 2 = p 2 + rf . 

p 2 = {$ — x) 2 + z 2 . 


(167) 

(168) 


remains constant during the integration. We integrate 
over rj from - if to + if, where : 

if < < z (169) 


Since according to our assumptions y < < p, we may 
expand (167) in a power series : 

T — p + (170) 

P 

and omit the remaining terms. 

If we substitute this expression in (166) and write : 

2-nr 2rrp rry 2 . 2 up . irrf 

cos -y = cos -m- cos - sm sin (17 1) 


the problem reduces to calculating the following two 
integrals : 

j + ^ dy cos and j dy sin y . . (172) 

or, if we use the abbreviation : 

(173 > 


we may write the integrals as : 






cos uHu and 




sin u 2 du 


The quantity : 


u ’ = V a/’ 


(174) 

(175) 
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may be regarded as very great compared with unity, in 
spite of the condition (169), so that : 

rf > > -\/Az (176) 

and, in writing this, we are supplementing our above 
convention concerning the domain of integration of C 
and S, respectively. 

The two integrals in (174) then run : 



[ + w 

cos u 2 du and / sin u 2 du 

— CO 


(177) 


These two “ Fresnel Integrals ” have, as we shall see 
in the next paragraph, a finite value— namely in both 

cases. If we substitute this value in (174), both integrals 
(172) assume the value for which, on account of 

(169), we may also write J-, and so we obtain from 
(166), in view of (171), the expressions : 

C = -^(c-s)+ . ,S = ~^(c + s) + . . . 

where we have used the abbreviations : 

l . . {m) 

’~Vzl ain? T de 

and, finally, from (165) : 

J = c 2 + s 2 -f (179) 

This reduces the problem to calculating the line- 
integrals c and s in (178), whose similarity to the integrals 
0 and 8 in (166) is immediately obvious. Of course 
c and s, 1 ike G and S, are pure numbers. 

A particularly striking feature of Fresnel’s integrals 
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( 177 ) is the fact that these integrals have a finite value 
although they stretch to infinity in both directions and 
although the functions that are to be integrated do not 
vanish at infinity, as is usually the case with such 
integrals, but retain finite values by fluctuating to and fro. 
The finite value of the integral results in spite of this 
because the maxima and minima + 1 and — 1 approach 
continually nearer each other as u increases, so that the 
fluctuations in the value of the integral become smaller 
and smaller. This peculiar behaviour is the mathematical 
counterpart of the physical circumstance that the con- 
tribution to the wave-function at the diffraction point A 
made by a strip clrj of the Huygen surface does indeed 
retain its order of magnitude, but changes its sign the 
more often the further the strip of surface is from the 
central line. Now what holds for the ineffective part 
played by the surface-elements da whose co-ordinate 77 is 
of the order of magnitude 77' in ( 176 ), must apply in a still 
greater measure to those surface-elements which are still 
further removed from the central line — that is, for all 
those surface-elements whose co-ordinates 77 are still 
greater than 77'. Hence it follows that the second com- 
ponent integral ( 164 ) of C, whose limits lie between 77 = 77' 
and 77 = 00 , and also that of S, can make no appreciable 
contribution to these quantities. The fact that the factors 
r and 8 in front of the cosine and sine can change appreciably 
does not come into question here. For their range of 
variation with 77 is so small comparatively that they can 
be neglected entirely for any interval of 77, which comprises 
many maxima and minima of the cosines and sines. 

Hence we may now omit the dots that were added in 
our formula and may also, in dealing with the Fresnel 
diffraction phenomena in the sequel, restrict the integra- 
tion over the Huygen plane to such elements da as lie 
near the central line, as determined by the conditions 
( 169 ) and ( 176 ). From now onwards the following 
relation holds : 

J = c 2 + s 2 . . . . (179a) 



98 


THEORY OF LIGHT 


chap. 


§ 47. We shall make a digression here by discussing 
Fresnels integrals. We have first to fill in a gap by 
calculating the definite integrals (177). For this purpose 
we start from Laplace’s integral : 



e ~ ax2 dx 3 where oc > 0 


( 180 ) 


We may write it in the form : 

/ + » 

e - “s' 3 dy. 

- oo 

Multiplying the two integrals together, we get : 

r + 00 r + 00 


y+oo y + oo 

L 2 ~ / dxdy . 

J — oo j — 00 


g ~ a(* 2 + 1 /) 


which, if we introduce polar co-ordinates : 

dxdy = p . dp . dcj) 
becomes transformed into : 







which works out to : 


Hence : 



( 181 ) 


Fresnel’s integrals are related to the Laplace integral 
in that i replaces oc in the latter. For the sake of brevity 
we shall forgo the proof that the equation oc =i also 
remains valid. Assuming this, we obtain from (180) : 

f + » 

L = / (cos x 2 — i sin, x 2 ) dx 

— CO 

and from (181) : 
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If we equate the last two expressions we arrive at the 
values already used in the preceding section : 


H- °° 


+ 1:0 


cos u 2 du = 
sin uhlu = 



(182) 


Since we also require the values of Fresnel’s integrals 
for variable limits, we shall carry the investigation a 
little further. Since ( + u ) 2 = ( — u ) 2 we easily obtain 
the integral with the limit 0 : 

/ oosu 2 du — / cos uhlu = / sin^ 2 ^ = (183) 

A oo A) 'o 2 \ ^ 

On the other hand, it is not possible to reduce the 
integral with one variable limit u to terms of elementary 
functions. We are compelled to resort to expansions in 
series. 

An expansion which is valid for the integral taken 
from 0 to u and which is always convergent is derived by 
applying the formulae of integration by parts : 


rU I U 

/ cos uhlu = u cos u 2 + 2 / sin u 2 du 
Jo Jo 


fsi 

Jo 


sin w sin u 2 


fU 

2 / u 9 
Jo 


COS llhlu 


If we repeat the operation of integration by parts on 
the integrals on the right, we ultimately get : 


/; 


cos u 2 du = cos u 2 . U + sin u 2 . T 
sin uhlu = sin u 2 . U — cos . T 7 


where we have used the abbreviations : 

9 , 9 , 2 . 2 . 2 . 2 


u 


rs u ° 


3. 5. 7. 9 


r =r 3 


2 . 2.2 . , 2 . 2 . 2 . 2. 2 u 
r « 7 + ,-,- r — 


3.5.7 


3.5.7.9.11 


(184) 

(185) 

(186) 
(187) 


These series converge for all positive and n 
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of u, because the quotient of two successive terms 
decreases to an unlimited extent as the order of the term 
increases. But on account of the rapid increase of the 
power indices they are quite unsuitable for calculation 
for u’ s of any considerable magnitude. 

In the case of greater values of u we find it more 
advantageous to use expansions in decreasing powers of 
u expansions which are often divergent but, as we shall 
see, are sometimes of great service in the calculations. 
Since a series with negative powers of u approximates to 
the value zero as u increases, it is expedient to start out 
from the integral which has u — co as its limit. We shall 
therefore now apply integration by parts in the following 
form : 

« 1 sin^ 2 /"sin uHu 

f cos v?du = « Q 0 su*--du = - + L ~2vF~ 

Ju> 0 Ju 

/ a . o 1 7 cos u 2 f 00 cos uHu 

I sin u 2 du =1 usinu 2 -- da = - J, 2 « 2 

Ju Ju 

and we obtain by applying the same process to the integrals 
on the right-hand side : 

[ cos u 2 du = cos u 2 . y,i — sin u 2 . a „ + R c ( 188 ) 
Ju 

r sin u 2 du = sin it 2 . y,i + cos u 2 . a* + R« ( 189 ) 
J U 


The series : 

„ 1 1.31 

2<7n ~ u~ 272 ' up 


. . . (- l )”" 1 


1 1 


1.3 



2y, ‘~2\t s 2.2.2 m. 7 

1.3 . 


. . (- 1 ) 


n-1 . 


(4:71 — 5 ) 1 


,,/dn-S 


( 4:71 — 3 ) 1 


2.2 .. . 


( 190 ) 


( 191 ) 


are divergent and hence have been broken off at a definite 
order number n. 
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The remainder terms are : 

& = (- l) ? ‘" Hr--- 


1.3. 

. . (4 n — 

J) 

r cos u 2 du 

2.2 . 

2 


/ i / 4 ti 

J u 

1.3. 

. . ( 4:71 — 

1) 

f a sin u 2 du 

2.2. 

2 


J„ 


They both increase beyond all limits as n increases. 
Nevertheless these series may often be used successfully 
in calculations. For closer inspection shows that for 
certain values of the order number n the remainder terms 
become very small and hence can be entirely neglected 
without introducing an appreciable error. 

For in the first place it is easy to see that both R c and 
R s are less than the expression R n which is obtained if we 
replace the cosine and sine in the integrals (192) and (193) 
respectively, by 1, that is less than : 


1.3 . . 

. (4 n - 1) 

/■“ <lu 

2.2.. 

2 

L 

1.3 . . 

CO 

1 

1 

2.2". . . 

2 

2 u l 


This is precisely the last term in the series for y n , a 
number which, in certain circumstances, is extremely 
small, particularly when u is great. The essential step 
is, of course, the choice of the order number n. If we 
allow n to increase from 1 onwards, R n first decreases, on 
account of the power (4?& — 1) in the denominator, but 
afterwards increases to an unlimited extent owing to the 
factorial term in the numerator. 

The minimum value of R n is that at the order number n, 
for which the ratio R n +i : R approaches most nearly to 1. 
This is the case when : 

Rn+ 1 __ (4 n — 1 ) . (4 n +1) 1 i __ •, 

‘ ~x~ “ — 2 — : — 2 ' ^ nearly ™ L 

That is, n nearly = . . . (195) 

The greater the value of u , the more the number of 
terms that must be given to the series y n and a n in order 
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that the remainder terms shall become as small as possible 
and the smaller these remainder terms become. If we 
neglect them we obtain approximate values for Fresnels 
integrals from (188) and (189). Such series, which, 
although divergent, approach a definite limit where a 
definite number of terms are taken but afterwards recede 
from this limit again, are said to be “ semi-convergent. 55 
They can often be used to better purpose than convergent 
series, but have the fundamental disadvantage that the 
approximation cannot be carried as far as we please, as in 
the case of convergent series, but comes to an end at a 
definite point. 

Nevertheless, even for moderate values of u the semi- 
convergent series y n and or* are much more convenient 
than the convergent series F and Z for determining the 
approximate value of Fresners integrals. For u = 2, for 
example, there are already very many terms to be taken 
into account in the convergent series (186) and (187) if the 
fluctuations in the sum are to be made even moderately 
small, whereas the remainder term R of the divergent 
expansion, in which, by (195), n has been put equal to 2, 
has, by (194), the value : 




which, compared with 1, is already a fairly small number 
and can in most circumstances be neglected without 
danger. 

§ 48. Reverting to our problem of determining the 
intensity of light J at the diffraction point A , we proceed 
to calculate the decisive integrals c and s in (178), by first 
expressing p in terms of the integration variable £ in 
accordance with (168). Since, as we showed in § 46, we 
need take into consideration only the points near the 
central line, we may assume £ — x to be small compared 
with z and write the following approximate value for p : 


( 196 ) 
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Thus : 


2 itp 2tj -z 

cosY = cos -)f cos 


- sin sin 
Az A Az 


If we substitute this in the expression (178) for c, we 
obtain the two integrals in it as : 


j cos ~ and j sin ^ ■ dg 


which are both to be taken between the limits | = 0 and 
| = oo. By introducing the integration variables : 

«-(*-£)• VS • ' • • (197) 


we may replace them by the integrals : 



cos u 2 du and 



sin u 2 du 


where u , the upper limit, denotes the value of (197) for 
£ = 0, that is : 


u = 



(198) 


This transforms the expression (178) for c into : 

1 ( 2ttZ n . 2772 a \ 

c -?;l 0 o ’-r'° _sm x- s ) 

where we have now used the abbreviations : 

rU rU 

C = l gos u 2 dic and$ = / sin^dw . (199) 

J— CO J —CO 


And in the same way we get from (178) : 

1 ( . 27 TZ „ 27TZ _A 

5 = V^v sm “- C7 + cos ^-' s i 

from which, by (179a), we get for the desired value of the 
intensity of light : 


J = ~ (C 2 + £ 2 ) 

77 


( 200 ) 



104 


THEORY OE LIGHT 


CHAP, 


This reduces the whole problem to the calculation of the 
Fresnel integrals G and S, in which u has the value given 
in (198). 

Since the intensity of light J depends only on u, it is 
quite independent of y, as is obvious, and we can restrict 
our attention to the plane of Fig. 12, y = 0. Moreover, 
on every curve : 

oc 

— const. 

yz 

we have that the intensity of the light is constant. This 
is a branch of a parabola which touches the a;-axis at its 
end and vertex 0 and is more or less strongly curved at 0. 
If we make the constant assume all values from — oo to 
4- oo and bear in mind that Vz is always positive, but x, 
the abscissa of the diffraction point A, is either positive or 
negative, we first get 2 = 0, x < 0 — that is, the points which 
lie immediately behind the screen, and then the parabolic 
branch which lies on the side of the screen ; further, for 
const. ~ 0 we have the positive 2 -axis, and then the 
branches which lie on the side of the aperture, as far as 
the positive £-axis : 2 — 0, x > 0 . 

To obtain a survey of the relation between the intensity 
of the light J on the position of the point A, it is therefore 
sufficient to displace the diffraction point A along any line 
2 = const, parallel to the #-axis. By (198) u is then a 
direct measure of the abscissa x of the diffraction point 
or of the distance OP of the central line AP from the edge 
of the screen. 

We shall now allow the diffraction point A to move along 
such a line (Fig. 12) from the right to the left across the 
whole plane of the figure, and shall inquire as to the 
intensity at such point as given by (200) and (199). 

1. For very great positive values of u (for which A is 
very far to the right, say at A^ we can set the upper 
limits of the integrals (199) equal to co , and we obtain 
from (182) and (200) that J = 1 — that is, at a great 
distance from the screen the intensity of the light is 
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just as great as if the screen were not present, as is to be 

anticipated. 

2. If u has moderate positive values, then, as we have 
seen, the expansions in semi-convergent series give us a 
useful approximation for estimating the values of Fresnel’s 
integrals. In this case we can even go down as far as about 
u = 1-5. For then, when n = 1, the remainder terms 
Ec and B„ which are to be neglected are, by (194) always 
smaller than : 

^ = i • i ‘ ra® = °' 074 - 


To be able to apply the formulae (188) and (189) we 
write instead of (199) : 


C 


/ + CO 


cos u~du 


,-foo 

L c 


cos u 2 du 


r -\~ 00 

S = / sin uklu 
7 — 00 


.-1-00 

/ sin u 2 clu. 
Ju 


On account of (182) we obtain the approximate values : 
C = — y n cos u 2 + a,i sin U‘ 


( 201 ) 


S = aJ ~ — y n sin u 2 — (Ju cos u 2 


Consequently, by (200), we get : 


/ = - 
77 


t r 


— 2\/ iryn COS [u 2 — 2J + 


2 V 77 a n sin ( u 2 — + yn 2 + cm 2 (202) 


The question arises whether the intensity of light J 
has maxima and minima in the region of u under con- 
sideration here. By (200) such points must result from 
the equation : 

C^ + S^ = 0 .... (203) 

du du 

or, since by (199) : 


dO 

du 


„ . dS . „ 

cos m 2 and = sm u 2 
du 


(204) 
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we have : 

«(*-*)-$ 

1 

2 V 77 

The expression on the right-hand side is very small 
for all the values of u that come into question ; conse- 
quently the equation will be satisfied by such values of 

the argument u 2 — ~ as approximate very closely to 

one of the values . . . and this gives the 

following approximate values of u : 

% = Vt* u * = Vt 1 u ° = “ 4 = V 1 ?’ • • • (2oe ) 

The first and smallest of these values % is equal to 
1*535 and lies just within the range of u which we are 
considering. 

In this range we thus actually have an infinite number 
of maxima and minima of intensity. We determine 
whether a value is a maximum or a minimum by substi- 
tuting in (202) . Since y n and a n are small positive numbers, 
the value of J is essentially influenced by the term with 

ovzsin (^ 2 — and the sign of the sine decides whether 

<7 is a maximum or a minimum, greater or less than 1. 
So we find that J is a maximum for the points u v u 3 , 
% . . . and a minimum for the points u 2 , u 4 , u 6) 

. . . From (202) and (190) we get as approximate values 
for the maxima and minima : 




13 5 
2 * 2 * 2 ' 


77 + 


■ ■■) 


(205) 


J = 1 ± 

the maxima being, by (206) : 

2 , _ 


wV TT 
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and the minima : 

T 1 J _ I _ ^ 

J '~ tVi ! 4 WT5’ • r 

As « increases both these points and the values of 
the maxima and minima continually approach each other 
until they finally blur each other. 

3. If the diffraction point A approaches so near to 
the 2 -axis or the shadow boundary of the screen that the 
abscissa u has rather small positive or negative values, 
say between — 1*5 and + 1-5, we find it advantageous, in 
calculating the Fresnel integrals, to use the expansions 
in convergent series. We therefore write (199) in the 
form : 

/0 /« 

C = / cos u 2 dw + I cos u 2 du 
J — « Jo 


/0 . /« 

8 ~ / sin u 2 du + / sin u 2 du 

J~ co J 0 

and by (183), (184), (185) and (200) : 


+ V^r . 2 . cos ( u 2 — 


+ Vw.r.sin(«*- j) +A 2 + F 2 . 


(207) 


If, to find whether there are maxima and minima of 
J in this region we again proceed as above, we arrive at 
the condition : 

cos (it 2 - ~) + E = 0 

or, by (186) : 

co s(« 9 - -§•!««+ = (208) 

In this equation the first summand is always positive 
within the region under consideration (— 1*5 <u< + 1-5); 
the second changes its sign with u, but its absolute 
value never reaches that of the cosines. To see this 
it is sufficient to calculate the series 2 for the extreme 
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case u = 1-5. We can then find the degree of approxi- 
mation which corresponds to the number of terms used 
by narrowing the neglected remainder term to the extent 
corresponding to (194). 

The question of the roots of the equation ^ = 0 may 


be illustrated still more strikingly if we represent geo- 
metrically the dependence of the function J on u ; this 
dependence is fully fixed by the simple equations (200) 
and (204) taken in conjunction with the boundary con- 
ditions for u = 0. This is done by taking G and S as 
rectangular co-ordinates of a point in a plane and inter- 
preting J as the square of the distance of this point from 
the origin. 

In accordance with what has been stated above, the 
expression (208) has the sign of the cosine in the whole 
region of u under consideration — that is, it is positive, 

or > o, and the intensity of the light increases steadily 

du 

as the diffraction point A passes from the negative to 
the positive abscissae, until it reaches the first maximum 
at u v In the immediate vicinity of the point A 0 of 
the boundary of the geometric shadow we have by (207) 


and (186) : 


J = 4 + 


u 


Vs 


(209) 


In A 0 itself we have J = -g . 

4. If, lastly, the diffraction point A moves a greater 
distance from the shadow -limit (u < — 1-5) towards the 
left, say to A 2 , we can again fall back on our semi-con- 
vergent series by writing : 

u — — u'(u' > 1-5) 
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and by using the formulae (188) and (189) for these expres- 
sions, omitting the remainder terms, this gives : 

J = i (y /t 2 + an 2 ) .... (210) 


whereas the condition for a maximum or a minimum 
reduces, by (203) and (204), to y a = 0, which is fulfilled 
nowhere. Hence the intensity of the light decreases 
steadily as u' increases. For very great values of u f 


we finally get : 



( 211 ) 


and for u = -oowe get J = 0. 

The above results also fix the intensity J for all other 
points on the xz-p\ane. For a definite parabolic branch 


of the family 


; const, passes through every point A 


of the straight line A 2 A 0 A x in question, and the intensity 
of light at A is the same at all points of any particular 
parabolic arc as far as the singular point O , where all 
the curves touch. Thus when the straight line A 2 A 1 
approaches the ar-axis, always remaining parallel to it, 
the maxima and the minima of J draw closer and closer 
together, without altering their values, and when the 
straight line coincides with the a*-axis, they all con- 
centrate in O, while the negative axis of x represents the 
parabolic branch of intensity zero and the positive axis 
of x represents that of intensity unity. This accords fully 
with the assumptions with which we started in § 43 in 
applying Huygen’s principle. 

§ 49. Having treated the general case of Fresnel’s 
diffraction phenomena, we shall now link up with the 
remarks made at the conclusion of § 45 and discuss the 
special case of Fraunhofer's diffraction phenomena, which 
have incomparably greater significance for practical optics 
and which are fortunately much more amenable to mathe- 
matical .treatment. We shall again start from the equa- 


tions (164) and (165), which hold for a plane diffraction 
screen on which light falls normally. A further special 
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assumption is that not only the point -source C but also 
the diffraction point A lie at infinity. We shall use the 
same notation as earlier — that is,, we shall choose the 
plane of the screen as the #?/-plane and shall take x, y } z 
as the co-ordinates of the diffraction point, and rj, £ == o 
as the co-ordinates of a point of the aperture. We shall 
take as our origin any point within the aperture. 

Since the distance of the diffraction point A from all 
points of the diffracting aperture is to be infinitely great 
compared with the linear dimensions of the aperture 
we shall assume the latter to be entirely in finite regions. 
Then the rays which travel from points in the aperture 
to A are all parallel and their direction is characteristic 
of the position of A. For the distance of A plays no 
part physically. Hence in this case it is sufficient to 
speak of the. direction of diffraction instead of the diffrac- 
tion point. This direction is defined by the two polar 
angles 9 and f in the equations : 


™ = sin 8 cos ” = sin 9 sin </>, ™ = cos 9 . (212) 

where x, y, z and r 0 , the distance OA, are infinitely great. 

When substituting these values in (104) we must note 
that in the denominator r is to be regarded as constant, 

as well as 8 , in the integration, while in the ratio ^ the 

variability of the numerator becomes of decisive import- 
ance. Thus : 

= (* - if + (y- Vf + (z~ if 

= r 0 2 - 2xg - 2 tjrj 

fr + yy 


If we substitute this value in (164), the expressions for 
0 and S may be written in the following form : 


G 


cos 


2^0 


c -f sin 


2rrr n 


S = sin, 


277T n 


cos 


27rr„ 
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where : 


1 9 f 

= X^ cos2 2' J da cos ^ + N 
p cos 2 -■ J da sin (oc£ + (3rj) 


An, 


(213) 


(214) 


2^ 2t t sin 0 cos cf> _ 2rry _ 27r sin 9 sin ^ n 
__ x -_ a>Ar; - A - P . 

This, combined with (165), gives ns for the intensity/ 
in the direction of diffraction (0, <j>) the value : 

J = c 2 + 5 2 (215) 

Since oc and (J remain constant during the integration, 
there is no difficulty in performing it. 

For the diffraction angle 9=0 (diffraction point A 
on the z-axis) oc = (J = 0, thus 5=0 also, and we obtain 
the “ axial ” intensity of light from (215) and (213) as : 


J a — 


AV 


(216) 


where a denotes the area of the aperture. What strikes 
us at once as unusual about this formula is that it is 
proportional to a 2 and not to ex; for the energy which 
passes from the source of light through the aperture is 
certainly proportional to cr. This apparent paradox is 
explained by the fact that J 0 does not directly denote the 
energy of radiation ; rather, a finite amount of energy is 
radiated only within a finite cone from the directions of 
diffraction (6, <f>), as we established in § 18. 

Since we are not concerned with absolute values in 
finding the intensity of the diffracted light, it is more 
convenient to reduce it to terms of the axial intensity 
and to write instead of (215) : 


J = J 0 .(O 2 + £ 2 ) . . . . (217) 

where, by (216) and (213), we have the newer meaning : 

C = cos 2 1 • - [dor cos (ag + $rj) 

2 a J 

S = cos 2 ~ 1 f da sin (ocg + ffy) . 

2 cr J 


( 218 ) 
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§ 50. As an illustration let us calculate the diffraction 
by a rectangular aperture of length l and width b, so 
that a = bl. If we take the centre of the rectangle as 
the origin 0 , we have in the integration . 


0 > 
“2 <f < 

and we get from (218) : 


<i 7 < 


0 sm ~2 Sm 2 


and from (217) : 


J = J 0 ■ cos 4 ^ ■ 


This, in conjunction with (214), is an explicit relation 
between the intensity J of the light and the direction 
of diffraction ( 9 , <j>). 

for 0 = 0 we get J = J 0 , as should be. This axial 
intensity is the greatest that occurs, which agrees with the 

circumstance that the function has its greatest value 

when x = 0. As a; increases this function fluctuates to 
and fro between a constantly decreasing maximum and 
the minimum zero at regular intervals, until it ultimately 
vanishes. Moreover, J is doubly variable in that two 
systems of equi-distant zero-points become superposed, 
The one is denoted by the values : 

2 ttX _ 2n 4y G 77 

a = X “TP TP 1P 

the other, by : 

. 2ttW 277 477 077 

P-Vp-PT- T' 
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Between these two zero points there is a maximum 
which decreases as the order number increases. To this 
there is added the gradual decrease of J as the angle of 
diffraction 9 increases. 

These variations of intensity may be observed by allow- 
ing the diffracted light to fall on a condensing lens. A 
definite point in the focal plane of the lens then corresponds 
to all the rays diffracted in a definite direction ( 9 , <j>) and 
our above discussion shows that the diffraction pattern 
that appears will be a bright cross with two systems of 
dark bands on both sides, situated symmetrically with 
respect to the edges of the rectangle. The narrower the 
breadth of the rectangular aperture and the longer the 
wave-length used, the more the bands are separated. 
If we have a long and narrow rectangle (slit), the bands 
of the one system merge together and only the other 
system, parallel to the length of the slit, remains. 

This diffraction picture does not, of course, represent 
an image of the diffracting aperture in the sense of 
geometric optics ; for we are not dealing with optically 
conjugate points here. Hence we must not expect a 
geometrical similarity between the diffraction image and 
the diffracting aperture. 

§ 50a. The practical importance of Fraunhofer’s dif- 
fraction phenomena becomes apparent when we use a 
great number of diffracting apertures congruent to one 
another. We assume a number N of such congruent 
apertures, of arbitrary form and arbitrary distribution 
but similarly placed on the screen. The intensity J of 
the diffracting light can then be reduced simply to the 
intensity J ± of the diffracting light of a single aperture. 

For we refer the co-ordinates £, rj of the points of the 
nth aperture to any arbitrarily chosen point On of the 
aperture as origin, by writing : 

i = in + rj = rjn -f r\ . . . (221) 

where and 7] n are the co-ordinates of On with respect 
to the general origin of co-ordinates 0, which may be 
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chosen anywhere, and we take the points On in all the 
other apertures at the corresponding point. Then each 
of the two integrals (218) resolves into N separate integrals 
each of which is to be taken over one aperture, with the 
given limits f and V, which are the same for each integral. 
We obtain : 

(7 = cos 2|.I^ j J dg' dr]' COS {d(gn + g') +P(p« + V)} • ( 222 ) 

The summation refers to the index n, the integration 
to the co-ordinates g’, V- These tw0 operations are 
quite independent of one another, and so can have their 
order reversed at will. 

If, analogously to (218), we use the abbreviations : 


cos 2 1 -^11 « cos («f + P 1 ?') = °i ^ 

c °s 2 f V // d ^' dr1 ' sin ^ + ^ = Sl - 


where, from now onwards, we use the index 1 to denote 
the quantities that refer to a single aperture, and likewise : 

cr = N . cq 

and : 

n=N 

Z COS (oc in + firjn) = On | 

n=l 
n-iV # 

2 sin (oc£n + ffyn) = Sn \ 


n = 1 

we get, by (222) : 


Similarly : 


G = {C N C 1 - SxSj 


S = ~ (S N C 1 + CxSj) 


and, by (217) : 


J n 


J = p ( G i 2 + s i 2 )( c » 2 + &»*)■ 


(224) 

(225) 
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Lastly, by (216) and (224) : 

J = J 0 i(Ci 2 + *5i 2 )((7^ + &v 2 ) . . (226) 

for which we may now write, by (217) : 

J ^JAOifi + Sit 2 ) .... (227) 

This simple relationship reduces the intensity of the 
diffracted light of the system of apertures to terms of 
that of a single aperture. It shows that the use of 
several apertures does not cause an extension and broaden- 
ing of the diffraction pattern in space, but only an increase 
in the intensity of light of the diffraction pattern due to 
a single aperture. But we shall see that the intensifica- 
tion factor assumes totally different values according as 
the apertures are distributed regularly or irregularly over 
the screen. For by (225) we have : 

N 

On 2 + Sn 2 = 2 cos 2 (ocf« + fan) 

i 

JV-l N 

+ 2 2 2 cos (oc in + (Jij«) cos (a£ n - + (Lj 

7i=ln'=n-\-l 

N 

+ 2 sin 2 (oc t;n + P^n) 

ft— i 

N~1 JV 

+ 2 2 2 sin (a£« + 1 3ij«) sin (a + firjn’) 

n—ln'=n-\-l 

= N + 22 2 cos{a(^-^)+P(ij»— ij*)} (228) 

7i=l ft' = ?t-{-l 

If the apertures are distributed quite irregularly— as, 
for example, in the case of a card which has been pierced 
by a pin indiscriminately, then the terms after the double 
summation sign have partly positive and partly negative 
values; this causes the sum to fluctuate violently for 
different values of a and (3, but it acquires no finite mean 
value. The result is that, by (227), the intensity of the 
diffracted light amounts, except for local fluctuations, to 
N - times the intensity of the diffracted light for a single 
aperture. 

The result that emerges is quite different if we use 
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regularly distributed apertures. If, for example, the 
points On lie so that for certain values of a and p all 
angles that occur in (228) are exact multiples of 2 t r, the 
double sum becomes equal to N (N 1)? and so, for 
great values of N, it becomes of the order N 2 . 

This peculiar behaviour is, of course, due to the fact 
that in the case of irregular distribution the rays that are 
diffracted in a certain direction by the JV-apertures do 
not interfere appreciably with one another, so that their 
intensities simply become added together, whereas in 
the case of regular distribution these N-rays can all he 
in the same phase, so that the absolute values of the 
amplitudes of the field-strengths add themselves together, 
and not the intensities. Now since the resultant intensity 
is represented by the square of the resultant amplitude, 
it amounts to N 2 times the intensity due to a single 
aperture. It is on this enormous increase in the intensity 
of the light that the great practical importance of Fraun- 
hofer’s diffraction phenomena depends. 

§ 51. We shall perform the calculation for the case of 
a grating. 

Let b be the width of aperture of a single slit, c the 
grating-constant — that is, the distance between two corre- 
sponding points of two neighbouring slits or the sum of 
the width of a slit and of the space between two neigh- 
bouring slits. Then we may set all the tj’s in (225) equal 
to zero, whereas : 

$1 = o, £ 2 = c, £3 = 2 c, ... = (N - l)c. 

We may sim plify the summation of the two trigono- 
metric series (225) considerably by recollecting that the 
expression. Cjfi + with which we are alone concerned 
here, is no other than the square of the absolute value of 
the complex quantity : 

QiaN f ]_ 

Cs + iSs = Z = Z e**"- 1 * = ; — y 

n= 1 n=l e — 1 

By multiplying this quantity with its conjugate imagin- 
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ary we obtain the desired quantity and from this, by 
(227), the intensity of the diffracted light : 


J = 


sin 2 JV 


a c 


Jv 


snr 


, ac 


(229) 


The maxima of intensity lie at : 


2ttx a 

oc = "> — = 0, 

Ar 0 


2tT 4w 
c ’ c ’ 


2mr 


(230) 


that is, at equal intervals which are directly proportional 
to the wave-length and inversely proportional to the 
grating constant. For axial light (0 = 0) the maximum 
is common to all wave-lengths ; a scattering of colours 
first shows itself at the next maximum, for n = 1, which 
gives the spectrum of the first order. 

The value of the maximum is the same for all orders 
and is equal to N 2 J V We must take care to observe 
that the factor J x in (229), the diffracted intensity of a 
single aperture, is not constant but depends in its turn 
on the angle of diffraction and the width of slit b in the 
manner given by (219). Hence gratings with the same 
grating constants often differ considerably from one 
another in their effects. If, for example, the interval 
between two neighbouring slits is exactly as great as the 
slit width (c = 26), the zero-points of J x in (220) coincide 
with the maxima of even order number n in (230), and 
hence all spectra of even order numbers are missing in 
the diffraction pattern of such a grating. 
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§ 52. The optics of crystals is founded on the same 
electromagnetic field equations (1), (2) and (3) as the 
optics of isotropic bodies. The only difference is that 
the relationship between the electric induction D and the 
electric intensity of field E, which is expressed in the case 
of isotropic substances by the simple proportionality (4), 
has a more general form in the case of crystals, in that 
the components of electric induction are certain linear 
homogeneous functions of the components of the electric 
intensity of field, that is : 


Ex — £]j Ex A £i 2 E,J -f- Ez 

Ey = £21 Ex + £ 2 2 Ey + £-23 Ez - 

Ez = £ 31 E x + e 32 Ey + e 33 Ez 


(231) 


where the constants £ depend on the nature of the crystal. 
The higher the degrees of symmetry the crystal has (II, 
§ 26), the more relationships there are between these 
constants. For the limiting case of an isotropic sub- 
stance £ u = £ 22 = £ 33 and all the other e’s are equal to 
zero, because (231) then becomes transformed into (4). 

But even when there is no symmetry at all only six of 
the nine quantities € are independent of one another. 
This follows if we apply the energy principle to any 
electromagnetic process that occurs within the crystal. 
For if we fix our attention on any portion of space in the 
interior of the crystal and apply the energy principle to 
it, exactly as was done in III, § 9, we get for the amount 
of energy that flows through all the elements da of the 
surface into space in the time dt, as in III (25) : 

dt. [ da S v — - dt [ dr . div 8 
J 121 J 
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and, by applying (3) and (1), this equation is converted 
into : 

^ ■ I dr (Exih + Eyby + EJD Z + H X H X + Hyily + UM Z ) (232) 

Hence by the energy principle this expression simul- 
taneously also represents the change in the time dt of 
the electromagnetic energy contained in the part of the 
crystal under consideration. It consists, as we see, of 
two parts : electric and magnetic energy, as in the case 
of isotropic bodies. The time differential of the electric 
energy-density is : 

~ {EdlDx + E,jdD y + EjD z ). . . (233) 

that of the magnetic energy-density is : 

- 1 - {HxdHx + HydHy + HzdHz) . . (234) 

4*77 

Now whereas the magnetic energy-density, as a com- 
parison with III, 3 shows, is exactly equal to that of ail 
isotropic body of permeability [jl = 1, the electric energy- 
density has a considerably more involved form. For if 
we substitute (231) in (233) and then integrate we get a 
quadratic form in the components of the electric intensity 
of field, whose coefficients are composed of the constants 
€. But if the integration is to be possible at all the 
coefficient of E x dE y must be equal to that of E y dE x — 
that is, we must have : 

€ i 2 ~ e 2i> likewise e 2 3 = €32, = e 13 . 

Then the electric energy-density is : 

(^n^ 2 d* e 22-®V 2 ~b CggE-s 2 + Qe^oExEy + 

2<t 23 E y Ez + 2<z 31 E z Ez) . (235) 

This expression may be considerably simplified by 
choosing as co-ordinate axes the three mutually perpen- 
dicular directions which form the axes of the ellipsoid 
represented by the tensor (235). These are the three 
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“principal axes” of the crystal. The electric energy 
then assumes the form : 

8 *. {ejEf + eoE,f + e 3 £ 2 2) (236) 

and the relation between the electric induction and 
electric intensity becomes more simple : 

Dx = Ex, Dy — D z = eJEz . . (237) 

The three constants e are called the “ principal dielec- 
tric constants 55 ; they are always positive. 

In the sequel we shall use the principal axes as the 
co-ordinate axes throughout. 

§ 53 . We next propose the general question : are plane 
linearly polarized waves possible in a crystal? And if 
so, what laws do they obey? To answer this question 
we form the expressions for the field-components of such 
a wave and substitute them in the field-equations ( 1 ) and 
( 2 ). 

Let n be the normal to the plane wave, and a, (3, y its 
direction cosines, thus : 

n = olx + py + yz . . . . (238) 

Further, let D be the value of the electric induction, and 
97 , £ its direction cosines, so that : 

B x = D . f , By = D . 77 , Dz = D . £ . . (239) 

and E the value of the electric intensity of field, and £', 
7 ] ', £' its direction cosines, so that : 

E x = E . f , E y =E .7]', E Z =E . £' . . (240) 

Lastly, let i? be the value of the magnetic field-strength 
and A, / a, v its direction cosines, so that : 

Hx = H .A, Bv~H.fl, rn^E.v . . (241) 

The condition for a plane linearly polarized wave is 
then that the quantities D : E and H, besides depending 
on the time t depend only on the length n, whereas the 
direction cosines that have been introduced are all 
constant. 



124 


THEORY OE LIGHT 


CHAP. 


With this assumption the field-equations (1) become 
transformed into the conditions : 

: dH ' N 

0 = c 0 ((3v - W -) ' -fa 

tS / * ^ dH 

t]D = c 0 ( r A - av) • -fa ^ 

0 = c 0 ( a/ x ~ (SA) • || 

AH = c o 0 J 'y-rp)-§f' 

aH = p,(r« - f y) • If » 

^ = Co (rp-^a).|| 

which contain the complete answer to the questions aboye 
proposed. In them we have denoted the velocity of 
light in vacuo by c 0 as the symbol c will be used a little 
later (§ 55) with another significance. 

§ 54. Concerning the directions of the various vectors 
which characterize a plane linearly polarized wave it 
follows from (242) that the electric induction is perpen- 
dicular both to the wave-normals and to the magnetic 
intensity of field : 

Din and D 1 H . . . (244) 

Further, it follows from (243) that the magnetic inten- 
sity of field is perpendicular both to the wave-normals 
and the electric intensity of field : 

H 1 n and H 1 E . . . (245) 

These theorems also contain the results which are arrived 
at by substituting the assumed values for the field-com- 
ponents in the field-equations (2). 

If we add to these the fact that, by (3), the flux- vector 
S is perpendicular to both the electric intensity of field 
and the magnetic intensity of field : 

S ± E and £ 1 JET 


(242) 


(243) 


. ( 246 ) 
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then all these relationships can be expressed in the single 
theorem that the magnetic intensity of field H is per- 
pendicular to all the other vectors that have been men- 
tioned ; that is, the latter all lie in a plane, the electric 
induction D being perpendicular to the wave-normal n 
and the electric intensity of field E being normal to the 
ray S. These relationships are represented graphically 
in Fig. 13, the plane of which has the magnetic intensity 
of field for its normal, its direction being towards the 
observer. This plane is called the “ vibration plane, 55 as 
in the case of isotropic bodies (§ 9). Thus the plane of 
vibration contains the normals to the wave, the ray, and 
directions of electric induction and electric intensity of 
field, which are perpendicular to the ray. The distinctive 


part played by the magnetic 
intensity of field in this relation- 
ship is, of course, due to the fact 
that it coincides with the magnetic 
induction, whereas this is not, 
in general, the case with the 
electric intensity and the electric 
induction. 


E 9 



Fig 13. 


If we denote the angle between 


the electric intensity of field and the electric induction 


8, then 8 is also the angle between the ray and the 


normal. Since the principal dielectric constants in (237) 


are positive, 8 < -. 

The physical meaning of the difference between the 
directions of the ray and the normal may be made clear 
by the following considerations. If we imagine the 
assumed plane-wave to be laterally limited — which is 
effected by selecting any finite surface area OA (Fig. 13) 
in a wave-plane, this area represents the cross-section of 
a cylinder of light which propagates itself onwards from it. 
It is formed by rays which proceed from all points of the 
cross-section in an oblique direction, since they deviate 
from the direction of the normal by the angle S. In 
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Fig. 13 the cross-section of the wave at different points 
of the cylinder, from OA to PB , is indicated by broken 
parallel lines. 

The value of § comes out from its definition as : 

cos 8 = ££' + yy + ££' > 0 . . . (247) 

or, if we also denote the direction cosines of the ray S by 
oc'py : 

cos S = oca' + pp' + yy . . . (248) 

Both the unaccented and the accented direction cosines 
form an orthogonal triplet with those of the magnetic 
intensity of field. 

§ 55. We shall now pass from the directions to the 
actual magnitudes of the vectors in question. For this 
purpose we eliminate the direction cosines by squaring 
and adding the .equations (242) and (243). The sum of 
the squares of the bracketed quantities in (242) is 1, since 
H±n. The same result is not, however, obtained for the 
corresponding sum in (243), since E is not ±n. Rather : 

(yy - m 2 + (£'« - f r ) 2 + (i'P - V *) 2 

= (f 2 + y 2 + n . (a 2 + p 2 + y 2 ) - (fa + f p + C'y) 2 

= 1 — sin 2 8 = cos 2 8. 

Consequently we obtain the relationships : 

1)2 = 0 o 2 (H ) 2 and ^ 2 = Co 2 (II) 2 - 0^-8. (249) 

In taking the square root, attention must be paid to 
the signs. They may be obtained from any particular 
case, for example, for that where the normal n lies in the 
#-axis, the electric induction D in the y-axis, and the 
magnetic intensity of field H in the 2 -axis. Then a = 1, 
P = 0, y = 0, £ = 0, y = 1, £ — 0, A = 0, y, = 0, v = 1, 
y = cos 8. Hence the equations (242) and (243) become : 

JD== - Co ^ i?= - C o• e0s8 ^^• • (250) 
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and, by (249), these relationships also hold in the general 


case. 

If -we eliminate H by differentiating the first equation 
with respect to t and the second with respect to n, we 


get: 


0 2 L 

dt 2 


c n 2 cos S 


02 E 
dn 2 


(251) 


On the other hand, we get from (237) : 

D ■ £ = ei ■ D . t) = e 2 . Erj', D . £ = e 3 . E£'. 


If we take the value of D obtained from one of these 
three equations and substitute it in (251) we obtain the 
well-known wave-equation (6) in the form : 


0 ^_ z d^E 

dt 2 q dn 2 ' 


(252) 


where q 2 , the square of the velocity of propagation, has 
the following value : 

q 2 = — cos 8 = ~ — cos S = y, — cos 8 


or, if we introduce the “ principal velocities of propaga- 
tion 55 : 

r 2 r 2 r 2 

-°- - a 2 , = Z> 2 , = c 2 . . . (253) 

€i 63 

gq 2 — £a 2 cos 8 , rj'q 2 = rjb 2 cos 8 , %'q 2 = £c 2 cos 8 . (254) 


The differential equation (252) is satisfied, as we know 
from II, § 35, by the expression : 


( 255 ) 

where /, the form of the wave, represents any arbitrary 

7b 

function of the single argument t . From this equa- 

Q 

tion we further get : 

3 E 10 E /occ\ 
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and, by substituting in (250) : 

E = ~ ■ cos S . E . . . . (257) 

D = ~ ■ cos 8 . E . . . . (258) 

Hence, as soon, as we have found a value for q which 
satisfies all three equations (254), we know that it is 
possible to have a wave of any form in the crystal with 
thecorresponding direction of propagation and polarization. 

§ 56. Let us now see whether it is possible to find 
expressions for the different constants which fulfil all the 
necessary conditions, above all, the equations (254). 

In the first place, it is easy to see that if the direction 
of the electric induction D is assumed in any arbitrary 
way, all the other quantities result uniquely. 

For if £, 77 £ are known then by multiplying the indi- 
vidual equations (254) by £, 77 , £ and adding, we get, in 
view of (247), that : 

q2 _ a 2g2 _|_ £ 2^2 ^ 2£2 . . . (259) 

and, further, from (254), the direction of the intensity of 
field E : 

£:rj':V =a^:b* v :c% . . . (260) 

But the directions of D and E fix the vibration-plane 
and hence also the normal direction H, as well as that of 
n and S. So the wave is completely determined except 
for its form, which remains arbitrary. 

In addition, we shall here derive only the expressions 
for the direction cosines oc, (3, y of the wave-normals. 

Since the wave-normal n lies in the same plane as the 
intensity of field E and the induction D, the ratio of 
their direction cosines is : 

a : P : y = (£' + x£) : W + xrj) : (£' + x£) 

where x denotes a certain number. But, on the other 
hand, the wave-normal is perpendicular to the induction. 
Hence : 


of + (fy + y£ = 0 . . . . (261) 
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whence we get, in view of (247) : 

x — — cos S 

Thus : 

a : [3 : y — (£' — £ eos §) : (V — V cos 8) :(£' — £ cos 8) 
from which we get, by (254) : 

a:(3:y = (a 2 - g 2 ) . £ : (Z> 2 - g 2 ) . t, : (e 2 - f) . £ . (262) 

§ 57. Let us next inquire as to how far a plane wave is 
determined if the direction of the wave-normal n is given 
in any way. Our problem is then to express all quan- 
tities in terms of oo, (3, y, and for this it is sufficient, as we 
have seen, to determine the direction cosines £, rj. £. 
These come out by (262) as : 




_ ft ■ V 

q 2 ' b 2 — q 2 ' c 2 — q 2 


(263) 


The velocity of propagation q, which still remains 
undetermined in this equation, is calculated according to 
(261) from the equation : 


— q 2 ^b 2 q 2 


+ 


- q 2 


= 0 


(264) 


As an exactly similar investigation in III, § 25, shows, 
this equation has two real positive roots in q 2 } and hence 
also in q. They lie in the two intervals that are bounded 
by the values of the constants a, 6, c. If we call them q 1 
and q 2 , we can, without restricting the generality of our 
case, agree that : 

a^q x ^b^q 2 ^c . . . . (265) 


With each of the two values q x and q 2 there is associated, 
by (263), a definite direction of the electric induction D 1 
andD 2 , and hence, by § 56, all the other directions are fixed. 

The two directions of induction D x and D 2 that belong 
to a definite wave-normal a, (3, y, v are very simply related 
to each other. For if in the equation (261) we substitute 
for rj, £ the values r) v £ x , but for oc, (3, y the values 
(262) with q 2 and f 2 , rj 2 , £ 2 , we get : 

(* ~ qi)i if* + (& 2 - 2 2 2 bi*?2 + (c 2 - = 0. 
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Similarly : 

(a 2 - gy 2 )^ + (& 2 - f li~)'<h r lz + (° 2 - S , i 2 )£i £2 = 0. 
From these we get by subtraction : 

+ ViVz + £ 1^2 ~ 0 ... (266) 

and : 

a2 ii^2 + WyiVz + c 2 £ i £ 2 = 0 * * • ( 2 67) 

The first equation states that B 1 _L D 2 , the second that 
D x J_ E 2 and that D 2 ± E v From this we also obtain 
further that H x || D x and H 2 || D v 

§ 58. The relationship between all these directions and 
the quantities of the corresponding velocities of propaga- 
tion q 1 and q 2 are realized graphically by considering the 
so-called Cauchy “ ellipsoid of polarization. 55 This is the 
ellipsoid : 

a 2 x 2 + b 2 y 2 + c 2 z 2 — 1 . . . (268) 

which has the semi-axes 

a o' c 

If we wish to find the two waves that correspond to 
an arbitrary given normal direction n , we draw the 
diametral plane of the ellipsoid perpendicular to n. The 
directions of the two axes of the elliptic section then 
represent the directions of the electric inductions D x and 
D 2 , and the reciprocal lengths of the corresponding semi- 
axes represent the corresponding velocities of propagation 
q 1 and q 2 . 

For if instead of x, y, z we introduce a new rectangular 
co-ordinate system x r y r z r having the same origin and with 
its s'-axis coinciding with the given normal direction n 
and its x'~ and y '~ axes coinciding with the directions of 
the electric inductions D x and J D 2 , the equations of trans- 
formation run : 

^ = £ x x' + + (xz r 

V = 7]!%' 4- rj 2 y' + §z r 

z = t x x f + l 2 y’ + yz f 

If we substitute these expressions in (268) and set 
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/ - 0 we get the equation of the elliptic section in view 
of (267) in the form : 

(a%2 + &V + c 2 £i 2 )*' 2 + + b\ z 2 + C 2 W 2 = j 

from which by (259) the values gq and g 2 result for the 
reciprocals of the semi-axes of the ellipse. 

Since the magnetic intensity of field II is perpendicular 
to both n and D each of the two axes of the elliptic 
section simultaneously represents the direction of the 
magnetic intensity of field, which belongs to the direc- 
tion of electric induction represented by the other axis, 
whereas the direction of the electric field-strength E is, 
by (260) given by the normal of the ellipsoid at the end- 
point of the axis which represents the electric induction. 
Finally the ray-direction S is represented by that tangent 
of the ellipsoid which lies in the vibration-plane defined 
by n 3 D and E. 

If a principal axis of the crystal, say the £-axis, is 
a normal to the wave, then the other two axes y and z are 
the corresponding directions of the induction with the 
velocities of propagation 6 and c. 

All these theorems, of course, remain valid for the 
case where a = b = c. The body is then isotropic 
optically and the polarization ellipsoid is a sphere, the 
elliptic section is a circle, every diameter is an axis 
and the velocity of propagation q = a. 

§ 59. Whereas hitherto we have assumed oc, (3, y as given 
and q x and q % as determined by them, the case may also 
arise where q x and g 2 are prescribed initially, of course 
within the limits (265), and the corresponding values 
oc, (3, y of the wave-normals are to be found. 

To solve this problem it is useful to write the equation 
(264) in the form : 

a 2 (6 2 - q*) (c 2 - q 2 ) + (3 2 (c 2 - g 2 )(a 2 - g 2 ) + y 2 (^ 2 - q 2 ) (& 2 ~ q 2 ) 
^(q 2 -qi)-(q 2 -q2 2 ) • • . • M 

This equation is an identity in so far as it holds, for 
any arbitrary value of q, provided only that q t and q % 
are the roots of (264). 
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a* 4 = 


According as we now set q 
(269) : 

(a 2 — (?i 2 ) • (o 3 — q^) 'i 
(6 2 - a~) • (c 3 - a") 
(& 2 - I 

"(c^Y j*) . (a* - W) 1 

(c 2 - q-i) • (c 2 - ? 2 2 ) 
(a 2 - c 2 ) . (d 2 - c 2 ) , 


a, & or c, we get from 


P 2 = 


(270) 


These relationships enable a 2 , (3 2 , y 2 to be determined 
as soon as q ± and q, are given. It is easy to convince 
oneself that by (265) the direction cosines a, (3, y are 
always real. They denote, according to their signs, eight 
different directions, each of which lies in a particular 
octant of the co-ordinate system of the principal axes, 
or, if we take any two in opposite directions as standing 
for one, we get four different directions. 

Then, by § 58, the values of a, (3, y also denote the 
direction of the electric induction and so forth that 
belong to q 1 and to q 2 . 

The special case where q 1 = q 2 is of particular interest. 
For by (265) we then have : 

?! = g 2 = & m) 


and, by (270), the direction cosines of the corresponding 
normals come out as : 


a n 


2 


a 2 - & 2 
r/,2 - c 2 ’ 


yz __ c 2 

Po = 0’ 7o 2 = a 2 _ c 2 


(272) 


Of these four directions that lie in the ccz-plane those 
two which form an acute angle with the z-axis — -namely, 
the axis of the smallest principal velocity of propagation, 
for which the following values hold : 


a„ = ± 


' a 2 & 2 

a 2 — c 2 ' 


Po = 0» Yo 


lb 2 - c 2 
\ a 2 — c 2 


(273) 


are called the “ optic axes ” of the crystal. The crystal 
is called “ optically positive ” if the two optic axes form 
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an acute angle with each other — that is, when y 0 > 


6 2 > 


a 2 + c 2 
2 ' 



(274) 


In the contrary case, the crystal is called “ optically 
negative.” 

In Cauchy’s ellipsoid of polarization the optic axes 
denote the normals of the two diametral planes which 

cut the ellipsoid in a circle of radius i Since the circle 

has an infinite number of axial directions, a wave 
whose normals coincide with that of the optic axis has an 
infinite number of directions for the electric induction, 
just as in the case of an isotropic body. Consequently 
we may say that with respect to an optic axis as 
wave-normal a crystal behaves just like an isotropic body. 
But an essential difference is that in an isotropic body the 
electric intensity of field coincides with the electric induc- 
tion and consequently also the ray with the wave-normal 
in this direction, whereas in the case of a crystal the 
electric intensity of field — that is, the normal of the 
ellipsoid of polarization at a point of a circular inter- 
section — in general by no means coincides with the electric 
induction — that is, with the radius of the circle at this 
point. Rather, each of the infinite number of directions 
of induction has its own particular field-strength and its 
own particular ray direction perpendicular to its field- 
strength, so that for each definite optic axis regarded as 
a wave-normal there are an infinite number of different 
rays, constituting a whole cone of rays. Among the 
generators of this cone there is also to be found the 
optic axis itself. For at the point where the circle of 
intersection meets the y-axis the normal of the ellipsoid 
coincides with the radius of the circle, and hence the 
electric intensity of field with the electric induction, and 
consequently also the ray with the wave-normal. 

If two of the three principal velocities of propagation 
a, 6, c are equal the crystal is called uniaxial, being 
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“positive, 55 by (274), if a = b>c, and “negative 55 if 
a>b = c. 

In numerical calculations it is customary to use, not 
the values of the principal velocities of propagation, but 
those of the “ principal refractive indices. 55 Bef erred to 
a vacuum these are : 


~ — % 


. (275) 


or, by (253) : _ 

n x = V*i w* w. = Ve 

where, by (205) : 


72 , ^ 72 « ^ 72 * 


Thus for 72 x = 22 2 the crystal is positively uniaxial, and 
for 72 2 = 22 s it is negatively uniaxial. 

The optical properties of a crystal are closely related 
to its elastic structure; for every crystallographic sym- 
metry also entails an optical symmetry. Hence the 
crystals of the non-symmetrical, monosymmetrical and 
rhombic system (II, § 26) are optically biaxial (topaz 
positive; aragonite and mica negative); those of the 
hexagonal and tetragonal systems are optically uniaxial 
(quartz, ice, zirconium positive, calcspar, turmaline 
negative) ; those of the regular system are optically 
isotropic (rock-salt, sylvine). 

§ 60. The physical significance of the velocity of 
propagation q is best seen graphically in Fig. 13 (§ 54), 
in which the wave-plane OA is depicted as having moved 
to PB in unit time. The magnitude of q is then repre- 
sented by the corresponding increase in the normal n; 
that is, by the distance OQ. If instead of measuring the 
velocity of propagation along the normal direction n we 
measure it along the ray-direction S, we get the distance 
OP, which we appropriately designate by q'. *• 

The general relation : 


g QQ 
q' ~ OP 


cos 8 


(277) 


then holds. 
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The appropriateness of our nomenclature is shown in 
the circumstance that each of the different relationships 
that have been set up between the direction cosines 
introduced and the velocities of propagation retain their 
validity if we exchange the accented and the unaccented 
direction cosines with each other and in addition replace 

q by q' by and a, &, c by ~, The direction of 

the magnetic intensity of field A, /x, v and the angle § then 
remain unchanged in the transformation. 

We have immediate proof of the correctness of this 
theorem if we consider that by § 54 the same laws hold 
for the accented direction cosines as for the unaccented, 
and that the equations (254) to (277), which serve to 
determine the velocity of propagation q , remain correct 
if we execute the specified exchanges in them. 

In this behaviour we see a dualism, a general law of 
reciprocity, which is valid throughout and which may 
often be profitably used to derive new relationships 
between the properties of a wave. With its help we can 
immediately enunciate a number of theorems which are 
no less important than those which have preceded. We 
shall here mention a few of the most important of them. 

Corresponding to every ray-direction oc', (3', y there are 
two velocities of propagation q\ and g' 2 , where : 


a^q\ ^ b ^ q\^c . 

They are the two roots of the equation : 

aV 2 6 2 (y 2 c 2 y' 2 

a 2 - q'* + 6 2 - g' 2 + c 2 - q'* U ’ 


(278) 

(279) 


These lengths represent the semi- axes of the ellipse 
which is the intersection of the diametral plane perpen- 
dicular to oc', (V, y with the ellipsoid : 


& , t + ?_ 2 - i 
a 2 ^ 6 2 ^ c 2 


. . (280) 


The directions of the two axes give the directions of the 
corresponding electric intensities of field, whereas the 
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direction of the electric induction is indicated by the 
normal of the ellipsoid at the end of an axis, and the 
direction of the -wave-normal by the tangent at the point 
of the ellipsoid in the vibration-plane defined by the ray 
and the electric intensity of field. The magnetic intensity 
of field forms the normal of the vibration plane, and so 
coincides with the other axis of the ellipse. 

Conversely we have for any two arbitrarily given values 
q\ and g' 2 , which fulfil the conditions (278), eight ray 
directions which are oppositely directed in pairs. They 
are represented by the values : 


,/2 — 

6 2 c 2 

(a* - 

. (a 2 

- q' 2 2 ) ' 



CC = 

' n' % n > 2 

~L 1 Z 2 

(a 2 - 6 2 ) 

.(a 2 

- c 2 ) 



p'2 = 

c 2 a 2 

2 ri' 2 
( 1 1 C 1 2 

(V - q’ i 2 ) 

(& 2 - c 2 ) 

. (5 2 
. (6 2 

<M 

l 1 


(281) 

'2 _ 

aW 

(c 2 - ?V) 

. (c 2 ■ 

-S' **) 



y 

' ?'iYa 2 

(c 2 - a 2 ) . 

. (e 2 ■ 

-& 2 ) J 



If q\ and 

q' 2 can 

be set equal to 

each other, that 

is, if 

by (278) : 









q'i = q'n 

= b 

. 

• 

(282) 


we get from (281) the two singular ray-directions which 
are called the “ secondary optic axes 57 of the crystal : 


, , c la 2 - 

a °~ ± b‘ 

[3'o = 0 = (3 0 
, _ a jb 2 - c 2 

7 0 b " V a * - c 2 



a 

= lYo 


(283) 


These are the normals of the two diametral planes 
that pass through the y-axis and which cut the ellipsoid 
(280) in circles of radius b. Comparison with the direc- 
tion cosines a 0 , (3 0 , y 0 of the primary optic axes shows 
that their angle with the 2 -axis is smaller — that is, it is 
enclosed by primary optic axes. In uniaxial crystals 
they coincide with the primary axes, for positive crystals 
(a = b) with the 2 -axis, for negative crystals {b = c) with 



PLANE WAVES 


137 


the positive and negative ;r-axis. To a secondary optic 
axis as a ray there belong a whole cone of wave-normals, 
the generators of which include the ray itself. Corre- 
sponding to each of these normals there is a definite 
vibration-plane which passes through it and the ray and 
which contains the electric intensity of field and the 
electric induction. 

All these theorems can, of course, also be deduced 
directly without the use of the law of reciprocity if we 
repeat the line of reasoning used from § 56 onwards by 
considering the unaccented quantities to be everywhere 
replaced by the accented quantities. 



CHAPTER II 
WAVE SURFACE 

§ 61 . The realisation of a plane wave is impossible in 
nature if only on account of its unlimited cross-section. 
But by applying the same considerations that were used 
in dealing with isotropic bodies (§ 5) we are able to realize 
plane waves to a degree of approximation that can be 
carried as far as we please. We take a point-source of 
light somewhere in the interior of a crystal which we 
suppose of any sufficiently great extent and assume that 
it begins to emit light at the moment t = 0 and we call 
the surface to which the light has advanced after the 
time t (which entirely surrounds 0) the wave-surface 
belonging to the point 0 as centre. If we choose t 
sufficiently large we can mark off on the wave-surface a 
portion of any magnitude whose dimensions will yet be 
small compared with its distance from 0 and which may 
therefore be regarded as nearly plane. The whole wave- 
surface is then composed of facets, as it were, of these 
plane portions, for each of which the laws of propagation 
of plane-waves hold. 

If there is difficulty in imagining the point-source of 
light in the interior of the crystal we may suppose the 
source to be situated at a point in the surface of the 
crystal, from which it propagates light into the interior. 
But in this case we do not obtain the completely closed 
wave-surfaces, but only a part cut out by the surface of 
the crystal; this does not, however, affect our subsequent 
arguments, as they refer essentially only to the considera- 
tion of a single plane portion of the wave. 

Such a plane portion of a wave propagates itself in the 
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costal, in accordance with the law which we deduced in 
the preceding chapter, with a velocity and a polarization 
which correspond with its orientation with respect to 
the principal axes of the crystal, exactly like the limited 
transverse section of an absolutely plane wave such as 
we indicated by PB in Fig, 13 (§ 54). In this case OP or 
OB, since PB < <0P , denotes the direction of the ray 
and OQ the direction of the wave normal ; and we have 
OP = q' . t, OQ = q . t. The fact that in the neighbour- 
hood of the light source 0 the wave surface can no longer 
be supposed to consist of sufficiently great plane portions 
causes no difficulty since we can always choose the time t 
so great that the distances OP and OQ can be set pro- 
portional to the time t. 

The above reflections furnish two different methods of 
determining completely the wave-surface corresponding 
to a sufficiently great time, say t = 1. 

The first method uses the ray velocity. Let us draw 
the distance q' in any arbitrary direction a'p'y'. The 
end-points P of all these distances form the required wave- 
surface. The second method uses the normal velocity. 
We draw the distance q in any arbitrary direction a(3y 
and construct the plane perpendicular to this distance 
through its end-point Q. Then all the planes that are 
constructed in this way form the tangential planes of the 
required wave-surface. The first method gives us the 
equation of the wave-surface in point-co-ordinates, the 
second method gives us this equation in plane-co-ordinates. 

It is possible to establish some properties of the wave- 
surface on the basis of these constructions, without having 
to make particular calculations. Since for every ray there 
are two values of the velocity of propagation, namely q\ 
and q' 2 , every straight line drawn from 0 as starting-point 
cuts the wave-surface in two points P 1 and P 2 ; that is, 
the wave-surface consists of two shells, the outer of which 
is composed of the points P x and the inner of the points 
P 2 . The tangential planes of the wave-surface at P x and 
P 2 are the two wave-planes that belong to the ray OP^P^, 
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their normals being the wave-normals. But there are 
four singular directions, namely the secondary optical 
axes (283) together with their opposite directions, which 
give only a single point of intersection with the wave- 
surface, at the distance b from 0 . Thus the wave-surface 
has four singular points situated in the ^-plane, at which 
the two shells meet. Since corresponding to each of 
the singular rays there are an infinite number of normal 
directions, the wave surface has an infinite number of 
tangential planes at each singular point n, that is, a 
whole cone of tangents, which give the surface in the 
immediate neighbourhood of the point a funnel-like 
appearance. This causes the outer shell to appear dented 
inwards at the point in question, and the inner shell to 
appear to have an outward bulge. The secondary optic 
axis itself is also a generator of the tangential cone. 

On the other hand, corresponding to any arbitrary 
direction of the wave-normals there are two tangential 
planes of the wave-surface at the distances q x and g 2 
from 0 ; these planes are perpendicular to the given 
direction. The lines connecting the point of contact with 
0 represent the directions and the velocities of propaga- 
tion of the two associated rays. But there are four 
singular tangential planes, perpendicular to the optic axes 
(273) and at the distance b from 0, for which the distances 
q x and g 2 both coincide with b. Corresponding to each 
of such singular tangential planes there are an infinite 
number of rays, and hence also an infinite number of 
points of contact with the wave-surface ; that is, the 
plane touches the wave-surface in a curve. Hence the 
above-mentioned funnel has a plane edge; it can be 
closed completely by a plane sheet. The optic axis itself 
is also a member of the generators of the cone of rays 
marked out by the curve of contact. 

All the above theorems, which follow of necessity from 
the developments of the last section, will be confirmed 
and illustrated graphically in the following special applica- 
tions. 
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§ 62. We shall now build up the equation of the optical 
wave-surface in the point-co-ordinates x, y, z for t = 1 ; 
for this purpose we shall use the first of the two above- 
mentioned methods. According to this a point of the 
wave-surface is represented by the co-ordinates : 

x = ra' , y = r(3' , z = ry 

where : 

r = va^ + y 2 + z2 = • 


If we substitute these values of a', (3', y', q' in equation 
(279) we obtain as the equation of the optical wave- 


surface : 


0, 2 X 2 lf-yi C 2 z 2 

a- — r 2 ' b‘ 2 — r 2 ft2 — r ' 2 


(284) 


The surface has 0, of course, as centre and divides into 
eight symmetrical octants. To obtain its order we mul- 
tiply (284) by the common denominator. We then get 
an expression of the sixth order, but easily see that it 
contains the factor r 2 . If we omit this factor, we get 
from (284) : 

« 2 & 2 C 2 _ { a 2(&2 + c 2 )x 2 + b 2 {c 2 + a 2 )y 2 + c 2 (a 2 + b 2 )z 2 } 

+ (x 2 + y 2 + z 2 )(a 2 x 2 + b 2 y 2 + c 2 z 2 ) = 0 . . (285) 

The surface is therefore of the fourth order and this 
agrees with the fact that a ray which starts from 0 cuts 
the surface in four points. 

The form of the optical wave-surface is best shown 
graphically by considering its sections with the three 
principal planes. Let us first take the xy -plane. For 
z = 0 (285) becomes : 

(a 2 b 2 - a 2 x 2 - b 2 y 2 ){c 2 - x 2 - y 2 ) = 0 . . (286) 

The curve of intersection of the wave-surface with the 
Ky-plane thus resolves into the ellipse : 

^ 4-^=1 

b 2 + a 2 

whose semi-axes are b and a, and into the inscribed circle . 


x 2 + y 2 ■■ c 2 
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of radius c. It is represented in Eig. 14. The points 
(P ) of the ellipse form the trace of the outer shell, the 
points (P a ) of the circle form the trace of the inner shell. 
The wave-surfaces corresponding to a ray S which has 
the points of intersection P 1 and P a are the tangential 
planes of the surface at P x and P 2 , which are, of course, 
perpendicular to the xy- plane. This simultaneously deter- 
mines the two mutually perpendicular planes of vibration 
which belong to the ray S and which pass through it and 
the wave-normal n. The vibration-plane at P x is the 
ry-plane, that at P 2 , where the ray and the normal 



coincide, is the perpendicular plane. The vibration-plane 
contains the electric intensity of field E and the electric 
induction D, the first being perpendicular to S and the 
second to «. At P x they both lie in the plane of the 
diagram (see figure). At P 2 they are both perpendicular 
to this plane and coincide. Conversely, the magnetic 
intensity of field H at P x is perpendicular to the diagram, 
whereas at P 2 it lies in the plane of the diagram. 

Exactly corresponding remarks may be made about 
the curves of intersection of the wave-surface with the 
other two principal planes. Let us first consider the 
zy-plane. Eor x = 0 (285) becomes : 

(b 2 c 2 - b 2 y 2 - c 2 z 2 ) (a 2 - y 2 - z 2 ) = 0 . . (287) 
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We therefore obtain as the trace of the external shell 
the circle : 

y 2 + z 2 = a 2 

of radius a , and as the trace of the inner shell the ellipse : 


whose semi-axes are c and b. The corresponding diagram 
is shown in Fig. 15. It differs from that shown in Fig. 14 
only in this respect, that for the points P 2 of the inner 
shell the vibration-plane falls in the plane of the diagram, 
whereas for the points P x of the 
outer shell it is perpendicular to 
that plane. 

§ 63. The most interesting cross- 
section of the wave-surface is that 
made by the #z-plane, that with 
the greatest and least principal 
velocity of propagation. For the 
two spherical sections which the 
plane y = 0 cuts out of the wave- 
surface (285) : 

(< c 2 a 2 — c 2 z 2 — a 2 x 2 ) 

( b 2 -z 2 -x 2 ) = 0 (288) 

that is, the ellipse with the semi-axes c and a, and the 
circle of radius 5, intersects in four real points, as is shown 
in Fig. 16. Since for every ray drawn from 0 to any 
such point the two velocities of propagation q\ and q' 2 
corresponding to it coincide with each other and with 6, 
those two of these four directions which make an acute 
angle with the z-axis form the two secondary optic axes 
of the crystal, and the trace of the outer shell of the 
wave-surface is partly represented by the two elliptic 
arcs which project beyond the circle and partly by the 
two circular arcs that project beyond the ellipse. The 
converse holds for the inner shell. 

Like the secondary optic axes so the primary optic 



Fig. 16 . 
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axes- also lie in the xz- plane, and, in fact, they form the 
normals of those tangential planes which touch the wave- 
surface along a whole curve. Their distance from 0 
amounts to q 1 = g 2 = b, the angle between the primary 
optic axis is, in accordance with (283), greater than that 
between the secondary optic axes. The curve of contact 
represents the plane edge of the funnel, described at the 
end of § 61, which intersects the plane of the diagram at 
the points A and B and for the rest runs perpendicularly 
to this plane; for this reason it is depicted only by a 
dotted line in Fig. 16. The straight lines drawn from 0 
to the points of the curve of contact constitute the cone 
of rays that belong to the optical axis OA ; the generators 
of this cone also include the optic axis. 

To find the form of the curve of contact we imagine 
the point G on the curve in the figure to be movable and 
starting from A to pass along the front arc of the curve 
to B and then to travel along the back arc of the curve 
to A again. In any position of G the corresponding ray 
is denoted by OG and hence the associated vibration- 
plane is denoted by the plane OA C and the corresponding 
electric induction by the direction AG which is perpen- 
dicular to OA and lies in the plane OAG. 

The length of the distance AG is obtained from the 
right-angled triangle OAG and comes out as : 


AG = OA. tan. A00 


= b . tan § 


sin § 
cos 8 


C °S % - s) 
cos 8 


or, by (247) and since the wave-normal a 0 , (3 0 , y 0 is 
perpendicular to the electric induction g, rj, £ : 


AC = 


h . 1 gof + ftp V + r<>£' 1 

££' + m + tt 


where : 


a o£ + fio 7 ) + 7oi = 0 


(289) 
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From this it follows by (260) that : 

ac _ h l a2a o £ + b2 %i + c2 ro£l 

a^TWTW~ 

and, by (259) and (289), we get, taking into account that 
here q = b and p 0 = 0 : 


AC = 


. . (290) 


The required curve of contact is completely determined 
by this relationship. For from it we get for every value 
of the direction cosine £ and hence also for every direction 
of the straight line that rotates about A the corresponding 
value for the distance AC. If C coincides with A then 
AG = 0 and £ = 0 ; that is, the tangent of the curve at 
A is perpendicular to the 2 -axis, and since it is also 
perpendicular to the optic axis OA it lies in the direction 
of the y-axis, perpendicular to the plane of the diagram. 
If C then moves towards B the angle between AC and 
the 2 -axis decreases and so £ increases up to a maximum 
which is attained when C coincides with B. For then, 
since AB is perpendicular to OA, we have £ = a 0 , and 


by (290) : 


a* — C“ 

~—j , — a o7o 


. (291) 


When the point C returns along the other arc of the 
curve the same values are traversed in the reverse direc- 
tion until C again coincides with A. The relation (290) 
becomes still simpler if we introduce the angle CAB in 
place of £. This is the angle between the direction AC , 
which has the direction cosines £, rj £ and the direction 
AB, which has the direction cosines — y Q , 0, oc x . So : 

cos CAB = — y 0 £ + a 0 £ 

and by (289) : 

cos CAB = (te 2 + Oo) { = I 
\« 0 V oc 0 

Combined with (290) and (291) this gives : 

AC = AB . cos CAB 
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Hence it follows that the triangle ABC is right angled 
at C and that the curve of contact AGB is a circle. But 
this does not imply that the cone of rays OC is a circular 
cone For the circular cross-section AGB is perpendicular 
to a generator OA of the cone. We get as the simplest 
expression for the size of the aperture of the cone of 

rayK AB a 2 - c 2 

tan a) = -Jr] ~ ' — p a °^° 


where u denotes the angle between the rays OA and OB. 
And by (273) we get : 


tan co 



c 2 ' 


(292) 


Exactly corresponding theorems may he derived from 
the law of reciprocity (§ 60) for the secondary optic axes 
and for the singular points of the optic wave-surface. 
At every singular point there are, corresponding to the 
singular ray connecting 0 to it, an infinite number of 
wave-planes and wave-normals, which form a cone of 
aperture co where, analogously to (292) . 



or : 

, V(a 2 - b 2 )(b 2 - c 2 ) 

tan co = " ~ 


ac 


(293) 



CHAPTER III 
NORMAL INCIDENCE 

§ 64. Since the reflection and the refraction of light 
at the surface of a crystal obey considerably more com- 
plicated laws than in the case of isotropic bodies it seems 
advisable to begin here by considering some rather simple 
cases. We shall first take the simplest case, where a mono- 
chromatic plane light-wave in air falls normally on the 
surface of a crystal which we shall suppose to have any 
arbitrary direction with respect to its principal axes. 

Since according to our assumption all points of the 
surface of the crystal are excited simultaneously in the 
same phase by the incident wave, the wave-plane and 
the wave-normal remain preserved also in the interior 
of the crystal. A peculiar feature, however, presents 
itself in that within the crystal there are two waves that 
propagate themselves with quite different velocities of 
propagation q 1 and q 2 , and they also have quite different 
directions for their field-strengths and inductions as 
obtained from the laws derived above for a definite 
direction of the wave-normals. On emerging from the 
crystal each of the two waves obeys its own laws of 
refraction and reflection. Hence we may use the differ- 
ence in the two velocities of propagation in the crystal 
to separate the waves entirely from one another, for 
example, by allowing one of them to be totally reflected 
when it impinges on the opposite surface of the crystal, 
while the other emerges refracted in the usual way. This 
is effected in a Nicol prism, which furnishes the simplest 
method of obtaining linearly polarized light from natural 
light. 
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Eor our discussion we shall now take a polished plane 
parallel crystal plate, so that even when the wave-planes 
emerge from the crystal they remain parallel to the 
surface. Since the electric inductions of the two waves 
are perpendicular to each other and to the wave-normals, 
they lie in the surface of the crystal and are also called 
the principal sections (optical) of the plate. They repre- 
sent the wave-planes of the two waves and are denoted 
by I and II in the plane diagram of Fig. 17, which is to 
be imagined as parallel to the face of the plate. We shall 
suppose the light-waves to enter the plate from the rear 
and to pass out in the direction pointing towards the 
observer. 

If the incident light is natural light the transmitted 
^ light will also be natural light if, as we 

A shall do here, we disregard losses due to 
/ reflection. But the result is different if 

/ we allow polarized light to fall on the 

/ plate, say by allowing the light first to 

f pass through a Nicol prism, which is in 
° Fig 17 this ease called a polarizer. The vibration- 
plane of the light which has been polarized 
in this way and is incident on the crystal is denoted in 
Fig. 17 by OP; let it form an angle cf> with the principal 
section I. On entering the crystal the light-wave resolves 
into two waves whose vibration-planes are I and II , 
whose velocities of propagation are q ± and q%, whose 
phase-difference is zero, and whose amplitudes have the 
ratio cos </> : sin <j> (§ 22). When these two waves have 
passed through the crystalline plate, of thickness D , they 
have the following difference of phase on emerging from 
the plate : 

■ • < 294 > 


where a> denotes the frequency, A 0 the wave-length in 
air, n 1 and n % the indices of refraction. Hence in the air 
they in general combine to form an elliptically polarized 
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wave (§ 26), the directions of whose axes depend essentially 
on the quantity A. It is only when A is a whole multiple 
of 7 t that linearly polarized light again results. If A is 

an odd multiple of the axes of the ellipse fall in the 

directions of the principal sections I and II ; if, in addition, 

i _ !I or — the light is circularly polarized. Hence 
r 4 4 

a crystalline plate of this kind and of suitable thickness 
(for example, a thin sheet of mica) provides us with a 
simple means of transforming linearly polarized light into 
circularly polarized light, or vice versa by being used as 
a compensator, and of analysing elliptically polarized 
light into its rectilinear components. 

If we make the light which emerges from the crystal 
plate in the forward direction pass through a second Nicol 
prism acting as an analyser , its direction of vibration OA 
(Fig. 17) making an angle iff with the principal section I, 
then of the wave I with amplitude cos <j> only the com- 
ponent cos <f> . cos iff and of the wave II with amplitude 
sin cf> only the component sin <j > . sin iff pass through the 
analyser; and these two waves, both of which vibrate 
in the direction OA, give the following intensity, according 
to (101), the light delivered by the polarizer being taken 
to be of unit intensity : 

J = cos 2 <f> cos 2 i/j + sin 2 <f> sin 2 iff + 2 cos <f> cos iff sin <f> sin ip cos A 
or : 

J = cos 2 (<f> — iff) — sin 2cf> sin 2 i/j sin 2 ^ • ( 295 ) 

A particularly interesting case is that where the direc- 
tions of vibration OA and OP are mutually perpendicular, 
that is, where the polarizer and the analyser are in the 

<c crossed 55 position. For then i/j = <f> -f and from 

(295) we get : 

J = sin 2 2 cf > . sin 2 — • 


(296) 
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In this case J vanishes simultaneously with A. Hence 
this proves that whenever light can pass through a plate 
placed between two crossed Nicols the phase-difference 
A and hence also - q 2 differ from zero. This gives 
us one of the simplest means of testing whether a sub- 
stance is anisotropic or not. The directions of the prin- 
cipal sections I and II can also be directly determined 
then, since by (296) the transmitted light has its maximum 

3 ^ 

intensity for <f> — -* and <f) = -r~, whereas for </> = 0 and <f> = - 

it vanishes entirely, as is quite natural. 

§ 65. So far in dealing with light-waves in the crystal 
we have considered only the wave-normals and not the 
rays. This is allowable so long as the 
cross-section of the wave may be assumed 
to be unlimited. But so soon as the 
cross-section of the wave is limited the 
circumstance asserts itself that the energy 
of the wave does not propagate itself in 
the direction of the wave-normal but in 
that of the ray as has been explained 
earlier, in § 54, by means of Fig. 13. 

Now since in general there are two different rays for a 
definite wave-normal the two waves into which the plane 
beam of rays that falls normally into the plate resolves 
itself wall sooner or later, depending on the size of its cross- 
section, separate and emerge from the opposite side of 
the crystal as two distinct beams of rays, as is indicated 
in Fig. 18. In this sense the crystal plate is doubly refract- 
ing even when the light is normally incident, although 
the wave-planes in all the rays, inside and outside the 
crystal, always retain the same direction. Here then the 
double-refraction is not, as in the case of isotropic bodies, 
determined by the difference between the substances in 
contact but by the nature of the crystal alone. 

§ 66. A particularly interesting case occurs when the 
crystal plate, say of arragonite, is cut perpendicularly 
to the optic axis. For then an infinite number of rays 



Fig. 13. 
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belong to the one wave-plane parallel to the surface 
of the plate, and the light, on entering the crystal, 
propagates itself not in two, but, if it was originally 
unpolarized, in all the directions of the generators of the 
ray-cone which corresponds to the optic axis as normal. 
The optic axis is also included among them. The plane 
of the diagram of Eig. 19, like that of Fig. 16 in § 36, contains 
besides the optic axis also the axis of the greatest and the 
smallest principal velocity of propagation (az-plane). The 
directions OA and OB denote, as in Fig. 16, the optic -axis 
and the ray of the cone that diverges most from it. Since 
the section of the surface of the plate, that is, the singular 
wave-plane with the ray-cone, is a circle, the rays that 
pass out into the air form the mantle 
of a circular cylinder, whose diameter 
AB is obtained, according to Fig. 19, 

from the angle of aperture cv of the g 

ray-cone and the thickness OA ± D ; /f \ / \ 

its value is D tan oj, where the value \ L Li 

(292) must be inserted for tan co. 

The phenomenon which has just been Fig. is. 

described is one of the most brilliant 

confirmations of theoretical crystalline optics ; it is called 
“ conical refraction/’ and in the particular case considered 
“ internal conical refraction,” because the conical resolution 
of the light occurs when the light enters into the crystal. 

§ 67. As in the case of every law that refers to the 
ray-direction in the crystal there is, corresponding to the 
phenomenon of conical refraction just considered, a 
reciprocal law which refers to the normal direction and 
which is founded on the circumstance that certain singular 
rays have not two but an infinite number of different 
wave-normals corresponding to them. We have then 
only to endeavour to produce a ray whose direction 
coincides with a secondary optic axis. For this purpose 
we may very well use the same arragonite plate as in 
the experiment described in the preceding section, but 
in the following way. 
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We place a point-source of light 0 at the surface of 
the crystal, say by covering the surface by an opaque 
screen having a fine hole 0 and then illuminating the 
screen from all directions. Rays from 0 will then enter 
the crystal in all directions. If we now also cover the 
opposite surface of the plate likewise by an opaque 
screen having a hole at P, then of all the entrant rays 
only a single one, namely the one OP, can emerge from 
the crystal. As we know, there are in general two 
normals, that is, two wave-planes corresponding to this 
ray, which make certain angles with them and pass 
out into the air according to a definite law of refraction 
which we shall learn more about in the next chapter. 

Thus at P we obtain two rays of light 



Fig. 20. 


emerging in different directions. 

But in the special case where the direction 
OP coincides with a secondary optic axis 
of the crystal, there are an infinite number 
of differently directed normals and wave- 
planes belonging to the ray OP in the 
crystal, and since each of them is refracted 
in a particular way on emerging from 


the crystal, the light from P propagates itself conically 
in an infinite number of different directions into the air 


(Eig. 20), a phenomenon which is known as “ external 
conical refraction ” because here the resolution of the light 


occurs on leaving the crystal. 

To find the appropriate direction P we must, of course, 
arrange to be able to move one of the screens, and must 
move it until the straight line OP connecting the two 


holes coincides with a secondary optical axis. 


Of course in the case of conical refraction, both internal 


and external, each of the rays of the cylinder or cone 
that emerges into the air has a perfectly definite plane 
of vibration and polarization, which may be obtained 
from the well-known laws of propagation and refraction. 



CHAPTER IV 
OBLIQUE INCIDENCE 

§ 68. The laws governing the passage of a plane wave 
of light from an isotropic body into a crystal, or the 
converse, may be derived by the same method as was 
used in the first part of this volume : we first form the 
expressions for plane waves in the interior of the two 
bodies in contact and then set up the boundary conditions 
for the surface of separation which express the con- 
tinuity of the tangential components of the electric and 
magnetic intensity of field. We get from this, if the 
incident wave is arbitrarily given, definite expressions 
for the refracted and the reflected wave. 

Since, in order to be able to retain the relationships 
that have hitherto been used, we have again taken the 
principal axes of the crystal as co-ordinate axes, we cannot 
here, as we could in the case of isotropic bodies, take the 
normal of the incident wave as the a;-axis; rather, we 
shall call the direction cosines of this normal k 0 , (3 0 , y 0 
and the velocity of propagation in the first body, which 
is isotropic, q 0 . Then, as in (8), an incident wave of 
definite polarization will be characterized by its wave- 
function (magnetic intensity of field) : 

f(t— . , . ( 297 ) 

where / represents an arbitrarily given function of a 
single argument. Let the normal of the plane of separa- 
tion, the incident normal, have the direction cosines u, v. w. 
We wish to find the wave-functions of the refracted wave : 

f t (t ca + < ?>y + y* . . ( 298 ) 
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and the reflected wave : 

j, ^ + (3' 0 y + y' 0 g 


(299) 


As in the case of isotropic bodies so here the problem 
falls into two parts. Firstly we have to determine the 
directions of the refracted and the reflected wave ; 
secondly we have to determine the wave-functions / x and 
f. Both these objects are accomplished by setting up 
the boundary conditions which express the equality of 
the electric and magnetic intensity of field at both sides 
of the plane of separation : 

ux + vy + wz = 0 (300) 


As in (15) and (16), this again gives us two equations 
which are linear and homogeneous in the quantities /, 
f v f, from which we find that f 1 and /' are proportional 
to/. 

§ 69. Let us first consider the directions of the refracted 
and the reflected wave. Since the proportionality of j\ 
and /' with / is to hold for all times t and for all points 
of the boundary surface (300), therefore, just as in § 7, 
the arguments of these functions must be equal to each 
other ; that is, at the boundary surface we must have : 

oiffc + (3 0 y + y 0 z _ mo + $y + yz _ + ft'o y + y' 0 z , am , 

To 2 2o 

The first equation leads to : 



If this were to hold for any arbitrary values of x, y, z 
the bracketed quantities would all have to vanish. But 

x, y, 2 are connected together by the condition (300). 
So we may also proceed by expressing z in terms of x and 

y, substituting in (302) and then setting the coefficients 
of x and y both equal to zero. It is more expedient to 
use Lagrange’s method of undetermined multipliers as in 
I, § 97. This consists in multiplying the expression (300) 
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by a certain constant A, adding the result to (302), and 
then setting the coefficients of x y, z individually equal 
to zero. This gives us the three equations : 


- - % + Xu = 0 I 

2 2o 

£_Po + Au=0 . . . . (303) 

2 % 

¥ - ^ + Xw = 0 

<1 % J 

These equations allow of a simple geometrical inter- 
pretation. In the first place the three directions repre- 
sented by the direction cosines all lie in a plane, that is, 
the normal of the refracted wave lies in the plane of 
incidence formed by the normal of the incident wave and 
the normal to the surface ; and secondly, if we multiply 
the individual equations with the corresponding direction 
cosines of the line of intersection of the incident plane 
and the boundary plane and then add, we get : 


sin 8 ___ sin 8 0 _ 


(304) 


where 9 and 8 0 denote the angles between the two wave- 
normals and the normal to the plane. Hence Snell’s 
law of refraction also holds for refraction at a crystal 
surface so far as it refers to the normal of the refracted 
wave. But here it has not the same significance as in the 
case of isotropic bodies, because the velocity of propaga- 
tion q is not known and itself depends on the desired 
angle of refraction 8 . Hence in order to calculate the 
direction of the refracted wave we must use the relation- 
ship between q and a, (3, y expressed by (264). This is 
done most simply and intelligibly by introducing the 
abbreviations : 




x 2 + y r + z 2 = t 2 



( 305 ) 
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where, of course, sc, y, z now have a different meaning from 
before. For then it follows from (303) that : 



and from (264) : 

= 0 (307) 


The first equation states that the point P represented by 
the co-ordinates x, y, z lies on the straight line which is 
parallel to the normal u, v, w to the plane of the crystal 
and which passes through the point A whose co-ordinates 

I are -°, Here OA 

[ — ^ y o #o ?o 






/ 

0 


Fig. 21. 


and OP (Fig. 21) denote 
the directions of the 
normals of the incident 
and the refracted wave and 
the reciprocal values of 
the corresponding veloci- 
ties of propagation q 0 and 
q. The second equation 


states that the point P lies on the surface (307). Com- 
parison with (284) shows that this surface may be regarded 
as the optical wave-surface of a hypothetical crystal whose 

111 , 

principal velocities of propagation are — , p — . This 


surface is also called the “ index surface ” of the real 
crystal whose principal velocities of propagation are 
a, b, c, because the principal indices of refraction n v 
n n 3 are, by (255), inversely proportional to the quantities 


Ur j l/j V. 

The index surface bears the same relationship to the 
wave-surface as the ellipsoid whose semi-axes are a, b, c 
(§ 60) bears to the Cauchy ellipsoid of polarization whose 

semi-axes are \ J, - (§ 58). Hence all the consequences 

Or 0 C 
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that follow from the general law of reciprocity also apply 
to the properties of the index surface. For example, 
the primary optic axes of the index surface are at the 
same time the secondary optic axes of the wave-surface 
and vice versa. 

§ 70. The preceding theorems give a simple geometrical 
method for finding the direction and the velocity of 
propagation of the refracted wave. Through the central 
point 0 of the index surface of the crystal draw a plane 
parallel to the plane of incidence and take it as the plane 
of the diagram (Fig. 21). This plane contains the normal 
ocoPo/o incident wave and the normal of incidence 

uvw. The intersection with the two shells of the index 
plane is denoted in the figure by two curves, one of which 

encloses the other. Now draw the distance OA = — in 

Qo 

the direction a 0 (3 0 y 0 and through A draw the parallel to 
the normal of incidence uvw. Every point of intersection 
P x or P 2 index surface then gives, when joined 

to 0 (that is, the line 0P l5 or 0P 2 ), the direction a|3y and 

the reciprocal velocity of propagation * of the normal 

of the refracted wave. So in general when the incident 
wave enters the crystal it resolves into two different 
refracted waves. Corresponding to each of these two 
wave-normals there is (1) a definite ray whose direction 
is represented by the normal to the index plane at the 
point P, and so in general leaves the plane of incidence 
and (2) a definite plane of vibration which passes through 
the normal and the ray, so that the direction of the 
electric and magnetic intensity of field and of the induc- 
tion is determined. 

If q Q is not too great compared with the principal 
velocities of propagation a, b, c it can happen that the 
point A lies outside the inner or even the outer shell of 
the index surface. In this case the point P 2 or both 
points P x and P 2 possibly become imaginary, and the 
conditions for total reflection obtain. 
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The method of construction just described is graphically 
verified if we apply it to the refraction of light by an 
isotropic body, which may always be regarded as a special 
case of a crystal. In the case of an isotropic body the 
two curves of intersection of the index-surface coincide 


with the circle of radius — , and from the relationships 

between the angles 8 0 and 9, which the directions OA and 
OP form with the normal of incidence uvw, and the 


lengths — and - of these distances we obtain Snell’s law 
% <1 

of refraction sin 9 0 : sin 9 = q 0 : q. 

§ 70 a. For the reflected wave the second equation (301) 
gives : 

(oc'o - oc 0 )tf + (P' 0 - Po )V + (/ o - 7o) z = 0 


and it follows from this, by reasoning fully analogous to 
that given above, that the normal of the reflected wave 
lies in the plane of incidence, and that its angle 8' 0 with 
the incident normal is determined by : 


Hence : 


sin 9' 0 = sin 9 0 


= 7 T 


do- 


That is, the angle of reflection is equal to the angle of 
incidence, just as in the case of isotropic bodies. 

§ 71. In the same way as we have here treated the 
refraction and reflection of a wave that emerges from 
an isotropic body and falls on a crystal we may also treat 
the converse case where a wave which advances within a 
crystal impinges on the surface of an isotropic body. 
We then obtain by the same method the corresponding 
laws of refraction and reflection, which distinguish them- 
selves in a characteristic way from those obtained earlier. 
In general, for example, there are two different waves 
which are reflected back into the crystal from the surface, 
whereas only a single wave which can be determined 
directly from Snell’s law escapes into the isotropic medium. 
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§ 72. We shall select for treatment only a few of the 
many applications which can be made of the laws which 
we have derived for the reflection and refraction of 
crystals. To simplify matters as much as possible we 
shall restrict ourselves to uniaxial ‘positive crystals, so 
that by § 59 : 

a = b > c (308) 

For this the equation (264) becomes : 

{(a 2 + p 2 )(c 2 - <? 2 ) + y 2 (a 2 - g 2 )} . (a 2 - g 2 ) = 0 
whose roots are : 

= a 1 . . (309) 

g 2 2 = (a 2 + (3 2 )c 2 + y 2 a 2 j 


The wave whose constant velocity of propagation is 
a = q 1 is called the “ ordinary ” wave and that with the 
variable velocity of propagation g 2 is called the “ extra- 
ordinary ” wave. 

If we call the angle which the wave-normal makes with 
the optic axis, the z-axis, 9, that is, y = cos 9 , and if, for 
simplicity, we omit the suffix 2 in q %i then we get for the 
velocity of propagation of the extraordinary wave, from 
(309) : 

q 2 — a 2 cos 2 9 + c 2 sin 2 9 (3 10) 


The equation (285) of the optical wave-surface becomes, 
for a = b : 

[a 2 x 2 + cfiy 2 + c 2 z 2 — a 2 c 2 ) . (x 2 + y 2 + z 2 — a 2 ) = 0 (311) 

So the surface resolves into the elongated ellipsoid of 
rotation whose semi-axes are c, c 9 a : 


x 


>2 n / 2 «2 

+ L 4 - — = 1 
c* ^ c 2 ^ a 2 


(312) 


and into the sphere of radius a, where : 

x 2 + y 2 + z 2 = a 2 

which touches the ellipsoid at the ends of its axis. 

The ellipsoid forms the inner shell, the sphere the 
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outer shell. The primary and the secondary optic axes 
coincide with the 2 - axis. 

In Fig. 22 the two wave-planes are drawn which corre- 
spond to a definite wave-normal OQ . They touch the 
wave-surface at the point P of the internal and Q of 
the external shell. The corresponding rays are OP and 
OQ . The vibration-plane of the extraordinary ray OP 
coincides with the plane of the figure, namely the prin- 
cipal section of the crystal that contains the optic axis; 
the vibration-plane of the ordinary ray OQ is perpen- 
dicular to the latter. Accordingly the electric induction 
at P also lies in the plane of the diagram, whereas that 
at Q is at right angles to it (cf. Figs. 14 and 15). ' 


z 



Fig. 22. 


The angle which the ordinary ray 
OQ makes with the 2 -axis is the same 
as that between the wave-normals, 
namely 9 ; that which the extra- 
ordinary ray OP makes with the 2 -axis 
is 6'. The value of O' can easily be 
obtained by reflecting that in Fig. 22 
the angle 9 denotes the direction of 
the normal to the ellipse at the ex- 
tremity of the diameter OP drawn 


in the direction 6'. Thus : 


sm V : cos U 


sin 6' cos 9' 


From this we obtain : 

a, a 2 cos 9 /q ia\ 

cos 9 = . A - = • • (o!4) 

V a 4 cos 2 9 + c 4 sin, 2 9 

For the angle 8 between the ray and the normal of the 
extraordinary wave we get : 

* w an & 2 cos2 9 + c 2 sin 2 9 / o I 

cos 8 = cos (9 - O') = _ _ - 7 — — W 15 ) 

V (P cos 2 9 + c 4 sm 2 9 

and for the velocity of propagation measured in the 
direction of the ray : 


(P cos 2 9 + c 4 sin 2 9 
a 2 cos 2 6 + c 2 sin 2 9 


( 316 ) 
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In calculating these quantities we may of course use 
the relationships which arise from the general law of 
reciprocity (§ 60). For example, from equation (310) we 
may immediately derive the following relationship : 

1 — " (317) 


l'2 


cos 2 O' sin 2 O' 
~r 


a 2 


which may easily be verified subsequently. 

§ 73. We shall now consider the passage of a plane 
monochromatic wave through a plane-parallel plate of 
a uniaxial positive crystal which is cut perpendicularly 
to the optic axis. Let the angle of incidence be <9 0 , and 
let the plane of incidence, which at the same time repre- 
sents a principal section of the crystal, be taken as the 
plane of the diagram in Fig. 23, 7 

in which the z - axis represents 
the optic axis and the incident 
ray is imagined as passing 
from below on the left through 
the plate up towards the right. 

Since we shall assume the 
wave, like the plate, to be 
unlimited laterally we need not take the ray-direction 
into special consideration and may restrict ourselves to 
considering the wave-normals. 

A ray of natural light which falls from below on the 
point A of the plate there resolves according to Snell’s 
law into an ordinary wave : 

sin On 


A 



r 


A 


Pig. 23. 


sin ~~ 

and an extraordinary wave : 

sin 6 0 
sin 0 o 


<h 


9o 

a 


.So 

% 


(318) 


(319) 


where q 2 and 0 2 are related to each other as in (310). 
Each of these two waves traverses the crystal and in 
passing out through the opposite face, which contains 
the point A , each yields a ray which is directed upwards 

M 
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towards the right which is parallel to the ray which 
extends from below at A , since the second refraction occurs 
in exactly the opposite direction sense to the first. Hence 
the two waves subsequently reunite again to form a 
single wave which is parallel to the incident wave. But 
since the path measured along the wave-normal is different 
for the two waves the phase -difference between them 
when they leave the plate is different from that on enter- 
ing. We proceed to calculate this phase-displacement 
A which arises from transmission through the plate. 

The normal of the ordinary wave in the crystal is 
represented in Fig. 23 by the distance AO, and the angle 
GAA' = 9 V At the point 0 the second refraction occurs 
which restores the original direction which made an angle 
6 0 with the 2 -axis. 

The normal of the extraordinary wave in the crystal 
is represented by the distance AB, the angle BAA' being 
equal to 0 2 . At the point B the ray again escapes into 
the air in the direction 0 Q . The wave-plane of the trans- 
mitted light is then GO' which is perpendicular to BC', 
and the wave-normal is BC'. The desired phase-dis- 
placement A is obtained if we fix our attention on two 
definite wave-planes, the one before its entrance into the 
plate, say that which passes through A, the other after 
its emergence from the plate, say that through G , and 
compare the changes which have occurred in the wave- 

function co(t ~ in the case of each of the two waves in 

the intervening space between the two wave-planes. In 
this intervening space the normal of the ordinary wave 
is the straight line AG, that of the extraordinary wave 
is the refracted straight line ABC', of which the part AB 
is traversed with the velocity q 2> and the part BC' with 
the velocity q 0 , no sudden change of phase occurring in 
the process of refraction. From this we get as the required 
phase-displacement : 


A = oj 


AB BC' AG\ 
^2 + ?0 * ' 


(320) 
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The distances AB, BO', AO as well as the velocity of 
propagation q 2 are obtained from the angle of incidence 
9 0 . For simplicity we assume 9 0 to be small and restrict 
our attention to terms involving d 0 K It then follows 
from (318) that : 


% 

and from (319) that : 


o, cos e 1 = 1 _ 


d 2 = ^. 
2 % 


and by combining these results with (310) we get 
q 2 2 = a? - (a* - c 2 )0 2 2 


2 2o 


^ = ~e 0 = e v 


Further, if D is the thickness of the plate : 

2 9o ! 


AB = -^3- = d(\ + ^« 2 

COS 0 9 V 1 ^ 9.n . 2 


AO = = AB 

cos 6 1 


whereas the distance : 


BC' = BC.6 0 


is of a smaller order than D . 0 O 2 and may therefore be 
neglected. 

The substitution of all these values in (320) leads to : 


A = 



h = — g2 ~ c2 o 2 
aJ 2 a ‘ g 0 2 ' 0 


or, if we introduce the principal coefficients of refraction 
% and w 2 and the wave-length A 0 in air we get : 


A - 1 1\ 

A 0 V% 2 w 2 2 / 


V • . . (321) 
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The extraordinary wave has its vibration plane and 
hence also its electric induction in the principal section, 
that is, in the plane of incidence, the ordinary wave 
being perpendicular to it. Hence the extraordinary wave 
is given by that component of the incident ray which 
vibrates in the plane of incidence -and the ordinary wave 
is given by the component which is perpendicular to it. 

If the incident light is linearly polarized, say by a 
Nicol prism used as a polarizer, whose plane of vibration 
makes the angle <f> with the plane of incidence, the phase- 
difference between the two waves on entering the plate 
is zero and their amplitudes are in the ratio cos <f> : sin <f> 9 
if we discard losses due to reflection. On leaving the 
plate the phase-difference is, as we calculated, A, so that 
the emergent light is in general elliptically polarized. If 
we subsequently allow the light to pass through a second 
Nieol prism, acting as an analyser, whose plane of vibra- 
tion makes the angle ijs with the plane of incidence, we 
obtain just as in (295) the intensity of the transmitted 
light referred to the light delivered by the polarizer as 
unity : 

J = cos 2 (<f> —xfs) — sin 2 <f > . sin 2ifj . sin 2 ^ • (322) 

where A is now given by (321). 

If the two Nicols are crossed, then we have as in (296) : 

J = sin 2 2 cf > . sin 2 ~ (323) 

If, finally, the light passes through a focusing lens 
which is parallel to the crystal plate the light rays meet 
at a single point of the focal plane of the lens, whose 
position is denoted by the angle 9 0 , which is proportional 
to the distance of the point from the axial point, that is, 
the focus of the rays that pass perpendicularly through 
the plate without being refracted. 

From this we can derive the phenomena that occur 
when a whole beam of nearly normally incident rays pass 
through a plane-parallel plate of a uniaxial crystal cut 
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perpendicularly to the optic axis and situated between 
two crossed Nicols. Here we have an infinite number of 
planes of incidence adjacent to one another, each corre- 
sponding to a different value of cf>. Hence the intensity 
of light J , besides depending on 0 O , also depends on the 
azimuthal angle <f> of the plane of incidence, and we 
obtain in the focal plane of the lens a black cross (c/> = 0, 

<£ = with concentric black rings superposed on it, whose 

radii are given by those values of A which are integral 
multiples of 2 tt. Eor other positions of the Nicol prism 
we obtain the value of J directly from (322). 

§ 74. So far we have dealt only with the direction and 
the velocity of propagation of the waves refracted and 
reflected by a crystal. 

If we now revert to the remarks made in §68 and 
inquire into the values of the wave-functions we must 
take into account, just as in the case of isotropic bodies, 
the special form of the boundary conditions. But the 
problem will be considerably more complicated than for 
isotropic bodies because here we are no longer able to 
treat the refraction and the reflection of the component 
vibrating in the plane of incidence completely separately 
from the component that vibrates perpendicularly to the 
plane of incidence. In general, rather, each of these two 
components will undergo double refraction when it en- 
counters the crystal. 

Nevertheless a moment’s consideration will show that 
there is a way of separating the problem into two inde- 
pendent parts and so to simplify the calculation con- 
siderably. Only, in doing so, we must take care not to 
start with a definite component of the wave incident 
from the isotropic body, but must fix our attention 
primarily on one of the two waves advancing in the 
crystal. Eor if, with the help of the theorems of the last 
section, the directions of propagation and the polarization, 
as w r ell as the velocities of propagation of the two refracted 
waves, are completely determined we may propose the 
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question : of what nature must the incident wave be in 
order that only one of the two refracted waves is formed ? 

So we now consider as known the direction of the 
normal n , of the electric and magnetic intensity of field 
and of the induction, as well as of the velocity of propaga- 
tion q of one of the two refracted waves, and we deter- 
mine the corresponding incident and reflected wave. This 
is accomplished by Setting up the boundary conditions, 
which express the continuity of the tangential components 
of the field-strengths. 

If we denote the magnetic intensity of field of the 
refracted wave by : 

H-fft-f ( 324 ) 


then the electric intensity of field of the wave is, by 
( 257 ): 


E = 


<1 


c 0 cos 8 


•A(< 


( 325 ) 


The direction of the electric intensity of field, which 
we also take as known, may be conveniently expressed 
by means of its direction cosines, £ rj, £, referred to the 
three following mutually-perpendicular directions : the 
normal of incidence (|), the line of intersection of the 
incident plane with the boundary face ( rj ) and the normal 
of the plane of incidence (£) (cf. Fig. 1). 

Then, of the two components of the electric intensity 
of field which are tangential to the boundary surface, 
that in the plane of incidence is 


c 0 cos 8 


‘A* 7 ? 


( 326 ) 


and that perpendicular to the plane of incidence is : 


c 0 cos § 


( 327 ) 


The direction of the magnetic intensity of field is per- 
pendicular firstly to the electric intensity of field which 
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has the direction cosines £, rj, £, and secondly to the 
wave-normal which lies in the plane of incidence which has 
the direction cosines cos 6, sin 6,0. Hence the direction 
cosines of the magnetic intensity are : 

£ sin 6 — £ cos 9 rj cos 9 — f sin 6 (328) 

cos 8 ’ cos 8 ’ cos 8 


The correct signs are obtained by considering a special 
case, such as that, for example, for which 8 =0, 6=0, 7j=l. 
Hence, by (324) and (328), of the two components of the 
magnetic field-strength which are tangential to the bound- 
ary surface that which lies in the plane of incidence is : 


-A' 


£ cos 9 
cos 8 


(329) 


and that perpendicular to the plane of incidence is : 

f Veosd-Jsind _ _ (330) 

Jl cos 8 

We next form the corresponding quantities for the 
incident wave, which we, of course, also assume to be 
linearly polarized. If we denote its magnetic field- 
strength by/, then by (8) and (7) its electric field-strength 

is / . — . If the vibration-plane makes the angle j> with 
c o 

the plane of incidence, then of the two components of 
the electric field-strength which are tangential to the 
boundary surface that in the incident plane is : 

f — • cos <j> . cos 6 0 . (331) 

c o 

and that perpendicular to the incident plane is : 

/f° sin^ (332) 

c o 

while the corresponding quantities for the magnetic field- 
strength are : 

— / sin <f> cos 6 0 .... (3 33) 



168 


THEORY OE LIGHT 


CHAP. 


and : 

/ . cos cf> (334) 

Finally, for the reflected wave we obtain fully analogous 
expressions, which we distinguish, as before, by accenting 
them, noting that q' 0 = q 0 and d' 0 ~ tt — 6 0 . For the 
electric intensity of field we obtain in this way the 
components : 

— f • ^ cos </>' cos 9 0 . . . . (335) 

c o 

and : 

/' sin <f>' (336) 

C 0 

For the magnetic intensity of field we get the components : 

f sin </>' cos 9 0 (337) 

and : 

f cos cf)' (338) 

At the boundary surface the tangential components of 
all the field-strengths are continuous, that is, the sum of 
each pair of field-strengths, of the same kind, of the 
incident and reflected wave is equal to the corresponding 
field-strength of the refracted wave. This gives us four 
boundary conditions; namely, firstly, for the electric 
field-strength in the incident plane, by (331), (335), (326): 

fq 0 cos cf> cos 6 0 - fq 0 cos <j>' cos 6 0 = ~ g/i’? (339) 

secondly, for the electric field-strength perpendicular to 
the plane of incidence, by (332), (336), (327) : 

fq 0 sin <f> + f'q 0 sin f A£ . . (340) 

thirdly, for the magnetic intensity of field in the plane of 
incidence, by (333), 337), (329) : 

-/sin<£cos0 o + /'sinf cos 0 O = “A-“y • ( 341 ) 
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and fourthly, for the magnetic intensity of field perpen- 
dicular to the plane of incidence, by (324), (338), (330) : 

/cos ^ + f cos f ~ A-- C ° S Sin 9 - ■ (342) 

If, as in § 7, we set : 

/' = jxf and A = fJ-J . . . . (343) 

the four equations may be satisfied by definite values of 
( f i, (j>\ ji and ; that is, there is in fact an incident wave 
of definite form and polarization which, besides producing 
a reflected wave of a definite but different form and 
different polarization, only calls up in the crystal the 
refracted wave here introduced. In view of (304) the 
values in question come out as : 


, , gsin(fl o + 0) 

an< ^ 7](sin0 o cos# o cos0 + sin0)— f sin # 0 cos 0 O sin 0 ^ ' 



si ncf) sin (0 O — 0) 
sin sin (9 0 + 9) 


(346) 


sin <j> . cos 8 . sin 2 9 0 
£ . sin (0 O + 9) 


For the special case where the electrical intensity of 
field of the refracted wave is perpendicular to the wave- 
normal and lies in the plane of incidence (£ = — sin 8, 
7 ] = cos 9, £ = 0, 8 = 0) the values of /l and assume 
the similarly named values (23) and (24) (angle of inci- 
dence 9 0 , angle of refraction 9), whereas the angles $ and 
<j>' both vanish. 

§ 75. If we next enquire how the incident wave must 
be constituted, for the same angle of incidence 6 0 , in 
order that only the second of the two refracted waves 
may come about, we may of course find the answer in 
the same way by starting out from the wave -function of 

the second refracted wave : g x (t — ^- 2 ). The correspond- 
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ing calculations then again lead to a definite incident 
wave g which vibrates in the azimuth and a definite 
reflected wave g r which vibrates in the azimuth If 5 
as in § 7, we then set 

g' = erg and g 1 = a x g 

the last four equations of the preceding section give the 
corresponding values of ^ 2 , (f>* 2 , a and oq, if we substitute 
for 0, 8 , f, r), £ the quantities calculated for the second 
refracted wave. 

For the special case where the electric intensity of field 
of the refracted wave is perpendicular both to the wave- 
normal and to the plane of incidence (£ = 0 , 77 = 0 , 
£ = 1 , 8 = 0 ) the values of a and oq become equal to the 
similarly named values (23) and (24), while the angles 

cf> 2 and <j >' 2 both become equal to 

From this we may now obtain the solution of the 
general problem — namely that of finding, for an arbitrarily 
given incident wave, the form of the two refracted waves 
as well as the form and the polarization of the one 
reflected w T ave. We resolve the incident wave into the 
two wave-components / and g , each of which produces 
independently of the other a refracted and a reflected 
wave according to a definite law. Only we may not 
here, as was possible in the case of isotropic bodies, 

resolve the vibrations in the azimuths 0 and ~ — that is, 

in and perpendicular to the plane of incidence, but rather 
in the azimuths </q and </> 2 , as given by (344), if we sub- 
stitute for 6, i, rj, £ the values calculated for the two 
refracted waves from the angle of incidence 0 O according 
to § 69. Then . / and oq . g are the wave -functions of 
the two refracted waves, whereas /x . / and a . g, which 
have the same normal, combine to form the reflected wave 
according to their azimuths ef>\ and <j>' 2 . If the incident 
light is linearly polarized, then so is the reflected light ; 
if it is natural light, then the reflected light is partially 
polarized. 
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It is obvious that the results thus obtained now allow us 
to treat all the special applications previously made for 
isotropic bodies also for crystals. For example, we may 
calculate the coefficients of reflection and transmission 
(pand 1 - p), normal incidence (0 O = 0), grazing incidence 

9 0 = “-j, angle of polarization (p = 0), total reflection 

(0 imaginary). It would lead us too far to discuss all 
these questions in detail here. 




PART THREE 

DISPERSION OP ISOTROPIC BODIES 




CHAPTER I 

FUNDAMENTAL EQUATIONS 

§ 76. Reverting from now on to the consideration of 
isotropic bodies we shall next occupy ourselves with a 
phenomenon which intrudes itself to a greater or a lesser 
degree and which we have been compelled to ignore up 
to the present because it has no place at all in the original 
theory of Maxwell — namely dispersion, We encountered 
this peculiar circumstance quite early, in § 9, and we 
have already envisaged the only escape from the difficulty 
in which it involves the theory : this consists in letting 
fall the assumption that space is filled absolutely con- 
tinuously by matter and introducing instead an atomistic 
point of view. 

In taking this step we essentially discard the point of 
view to which we have hitherto clung without exception 
in all our discussions and enter a field in which, to be 
able to make progress at all, we are from the very outset 
at the mercy of more or less arbitrary hypotheses. This 
difficulty is inherent in the nature of the problem and 
may be mitigated to a certain extent by assuming the 
atoms or, expressed more generally, the smallest element- 
ary bricks ( Bausteine ) of matter to be extremely tiny 
and extremely numerous. While this assumption pro- 
vides us with the possibility, on the one hand, of attri- 
buting to the material body which is composed of atoms 
new kinds of properties which have not been expressible 
in the theory hitherto developed, the body yet behaves 
practically, provided that the phenomena in question do 
not occur in dimensions of too small an order of magni- 
tude, as if it were perfectly homogeneous. A good 
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example of this is given by the propagation of sound in 
a body constituted of atoms, for which, provided the 
wave-length is not too small, exactly the same laws hold 
as in an absolutely homogeneous body of corresponding 
density and elasticity. Regarded from this point of view 
the whole theory hitherto discussed appears as a limiting 
theory which is valid for comparatively slow and “ coarse- 
grained 55 events, and our present task is to extend the 
range of validity of the theory a ; step further towards 
the side of more rapid and “ finer-grained ” events. 

§ 77. To preserve connection with what has gone before 
we retain, as in crystal optics, the fundamental equations 
(1), (2) and (3) and only introduce a more general relation- 
ship in place of equation (4), which expresses the pro- 
portionality between the electric induction 2> and the 
electric intensity of field E. Such a relationship is sug- 
gested to a certain extent by a graphical physical inter- 
pretation which can be given, by III, § 26, to the difference 
between induction and field strength. For if we set, by 
111(141): 


I) - E 
4tt 


= M . 


(348) 


M may be regarded as the electric moment of the unit of 
volume, that is, as the sum of the moments of all the 
infinitely small electric dipoles which are to be supposed 
contained in the volume 1 (taken sufficiently small). 

The advance which we now intend to make beyond 
this view consists in ascribing individually to each of 
these dipoles a real physical existence, that is, we no 
longer assume an infinite number of infinitely small 
dipoles in immediate contact with one another, to enquire 
into the individual size of which has no meaning, but 
rather we assume that their number, size and distances 
from one another are real and finite. The most essential 
feature of this hypothesis is that now there is only a 
single medium which is continuous in the true sense, 
namely, the vacuum in which all electromagnetic field 
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effects occur, whereas the material bodies share only in a 
secondary way in these processes, in that their electrically 
charged components, the ions and electrons, fly about 
individually in the vacuum. This endows the critical 
velocity c, the velocity of propagation of light in a 
vacuum, with a much deeper meaning than it had in the 
original theory, where it was only one velocity of pro- 
pagation among many. 

§ 78. If we now enquire what general relationship may 
take the place of equation (4) if we adopt the new view 
suggested above, we must first reflect that by equation 
(4) the electric moment of all the dipoles contained in 
volume 1 is proportional to the simultaneous field-strength. 
As we now wish to ascribe a real existence to the dipoles, 
we shall also have to give them a certain independence ; 
that is, we shall not demand that their electrical moment 
owing to some intimate coupling shall instantaneously 
participate in all the fluctuations of the electric field that 
acts on them, but rather we shall generalize the law for 
the relationship between the electrical moment and the 
exciting electrical intensity of field by distinguishing 
between the fluctuations of the electrical moment of the 
dipoles and the fluctuations of the exciting electrical 
intensity of field. The former will, indeed, be determined 
by the latter, not by means of a simple relationship of 
proportionality, but by means of a special equation of 
vibration in which also the individual properties of the 
dipoles, in particular their inertia, play a part. There 
are various possibilities for choosing the form of the 
equation of vibration. For example, we may assume the 
dipoles to be rigid and capable of rotation and having 
constant electrical moments. Their moments of inertia 
then enter into the equation of vibration and the elec- 
trically neutral state is characterized by the completely 
chaotic arrangement of the dipoles. 

The phenomena of optical dispersion here under dis- 
cussion can be best accounted for, at least to a first 
order of approximation, if we represent the vibration of 
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a dipole by means of the equation for the forced oscil- 
lations of a point-mass in carrying an invariable electric 
charge e and oscillating about a definite position of 
equilibrium which has an equal and opposite charge. In 
vectorial form this equation runs (II, 203) : 

mr + mco 0 2 r = e . E r • . . . . (349) 

Here <o 0 denotes the proper frequency of the oscillating 
dipole, r its distance from the position of equilibrium, 
and hence, by III, § 26, e . r its electric moment, and, if 
there are N such dipoles in unit volume : 


N . e . r ~ M (350) 

is the electric moment per unit of volume. 

The exciting field-strength E' must be carefully dis- 
tinguished from the total field-strength E. For the former 
arises only from the electric charges which are situated 
outside the dipole in question, since a dipole does not 
excite itself ; the latter contains besides, however, the 
field due to the dipole in question itself. The difference 
between E' and E would be vanishingly small only if a 
sufficiently small volume v with the electric moment M . v 
were to produce no appreciable field in its interior. But 
this is in no way true. Rather, the intensity of field 
which an arbitrary small sphere, polarized uniformly with 
the moment M per unit of volume, produces in its interior 
is everywhere in that interior, according to III (189), 

equal to — ^ M, that is, independent of the size of the 

u 

sphere. It is only when we add this amount to the 
exciting field-strength E r that we obtain the total intensity 
of field E, thus : . 

jE'-~M = JS .... (351) 

If we combine the last three equations we get the 
relationship between the electric moment M and the 



i. ■ FUNDAMENTAL EQUATIONS 179 

electric intensity of field E, namely, the equation of 
vibration, in the form : 


mM + ( mco 0 3 - ~ Ne 2 )M = Ne 2 E . (352) 


and, if we combine this with (348), we arrive at the 
desired generalization of the fundamental equation (4). 

The expression on the left-hand side of the equation of 
vibration is often supplemented by a damping term in 
M and a positive coefficient, as in I (20) or in III (375), 
which allows us to take into consideration any damping 
that may be caused either by the loss of electromagnetic 
energy or by collisions with neighbouring oscillators ; we 
shall not, however, consider this extension here. 

If we use the abbreviations : 


477 Ac 2 
m 


(353) 

(354) 


the equation of vibration assumes the simpler form : 

M + o/ 0 2 M = — . E . . . . (355) 


The equations (1), (2), (3), (348) and (355) form the 
basis of the theory of dispersion which we have described. 
This theory is largely due to H. A. Lorentz. Our first 

deduction from it is that co 0 2 >|, which imposes a certain 

upper limit on the density of distribution of a definite 
kind of oscillator. For otherwise E and M would have 
different signs in the statical field. 

§ 79. In order first to link up with the original theory 
of Maxwell we consider processes which occur so slowly 
that in the equation of vibration (355) the term in M 
may be neglected in comparison with M ; this amounts 
to the frequency of the vibrations of the dipoles, caused 
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by the exciting wave, being small compared with their 
proper frequency. It then follows from (355) that : 


M = 


477 0>n' 2 


E 


and, if (348) is taken into account : 


D = (l+ , 2 >E 

<jj o 

Here, therefore, as in (4), we have strict proportionality 
between the electric induction and electric intensity of 
field, the factor of proportionality, namely, the dielectric 
constant, having the value : 


6 = 1 + (356) 

a>'(f 

This is the region to which the consequences deduced 
in § 9 refer. _ The dispersion vanishes entirely and, by 
(22), we get for the refractive index : 


n 2 = 1 + ■ . . (357) 

co' 


The formulae (356) and (357), which, as we see, contain 
a physical explanation of the nature of the dielectric 
constants, and hence also of the electric induction, may 
be tested experimentally by measuring the dependence 
of the dielectric constants or of the refractive index on 
the density of distribution N of the dipoles. We must 
note, however, that, by (353) and (354), a as well as co 0 ' 2 
depends on N. 

If we substitute the values in question in (357), we 
may write : 


n 2 - 1 _ „ 47re 2 

n 2 + 2 iY * 3mcV 


(358) 


Since the constants e , m and co 0 are independent of N, 
this relationship is very convenient for determining the 
connection between the quantities n and N. Indepen- 
dently of H. A. Lorentz it has also been derived by the 
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Danish physicist L. Lorenz and is satisfactorily confirmed 
by measurement. 

§ 80. For treating the general case it is advantageous 
to ehminate the magnetic intensity of field H and the 
electric induction D from the equations (1) and (348). 
In view of (2) this gives 

E + 4 ttM = c 2 AE (359) 

This equation together with (355) may serve for obt aining 
solutions of the dispersion problem. 



CHAPTER II 
PLANE WAVES 


§ 81 . Passing on now to develop the most important 
laws of the electromagnetic processes that occur in a 
dispersive body from the fundamental equations that 
have been established, we follow along the line of reasoning 
adopted in § 3 in treating isotropic bodies and first con- 
sider a plane wave; that is, we assume that, of the 
three space co-ordinates x, y, z, only the co-ordinate x 
comes into action. We also assume the wave to be 
linearly polarized, the vibrations being in the y-direction, 
as formerly for the /-wave. Then E x = 0, Ez — 0, 1I X =Q > 
H y — 0, and the equations (355) and (359) reduce to : 

M,J + a) 0 ' 2 My = ~E t . . (360) 

and: 

Ey + inMy = • • (361) 

Compared with the differential equation (6) for a body 
free from dispersion this equation exhibits the funda- 
mental difference that it has no longer to satisfy the 
differential equations without introducing a definite 
assumption regarding the form of the wave-function. 
Nothing remains but to set up a particular solution of 
the equations and then, by combining a sufficient number 
of particular solutions, also to write down the general 
solution, exactly as was done in dealing with the problems 
of total reflection in § 12. 

As there, we also assume : 



182 


( 362 ) 



chap. ii. PLANE WAVES 183 

and, correspondingly : 

M y = o ce v c/ (363) 


This denotes a singly periodic wave with the velocity 
of propagation : 


c 

n 


(364) 


Hence we also call the constant n the refractive index, 
as in the case of bodies free from dispersion. A pecu- 
liarity in this case, however, is, as we shall see, that n 
can also be imaginary. Then, as in the case of total 
reflection, the wave is spatially damped. 

The assumed values of E y and M y , when substituted in 
equations (360) and (361), satisfy them in actual fact 
provided that the following relationships hold between 
the constants : 

n 2 _ i _ 4 . 77 . 0 c — — 2 — 2 . . (365) 


Since the constants a and a / 0 are determined by the 
constitution of the body, there is accordingly for any 
frequency co of the wave a definite index of refraction n 
and a definite ratio a of the electric moment M to the 
field-strength E. It then follows from ( 1 ) that the 
magnetic field-strength is given by : 

Hz = nE y (366) 

§ 82. Let us next enquire into the dependence of the 
refractive index n on the frequency co of the wave. By 
(365) n is either real or purely imaginary. In the former 
case we choose n positive; that is, we make the wave 
proceed in the positive direction of x. For very small 
values of co we obtain the relationship (357), as is obvious. 
If co increases, n also increases, and, in fact, to an unlimited 
extent if co approaches indefinitely close to the value o / 0 
(cf. Pig. 24). After this value has been exceeded n 
becomes purely imaginary and equal to — i#c, where k 
is real and positive, since otherwise E y would become 
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infinite for x = oo . This lasts until o> reaches the value 
Voj'q 2 + a. In this interval k decreases continuously 
from oo to 0, as is indicated in Fig. 24 by the dotted 
line. From then onwards the refractive index again 
becomes real and increases from zero asymptotically to 
the value 1, which it attains for a> — oo . The result is 
that the whole spectral region from a> = 0 to o> • = oo is 
divided into three different parts (Fig. 24) by an 
£C absorption band 55 which stretches from o> = o/ 0 to 
o> = 's/cdq 2 + a. To the left and to the right of the 
absorption band there is no absorption, but the refractive 



index always increases as the frequency increases ; the 
dispersion is, as we say, “ normal.’ 5 But there is the 
essential difference that on the left, for smaller frequencies, 
the index of refraction is always greater than ^/7 (cf. 
(356)), whereas on the right, for greater frequencies, it is 
always less than 1. Within the absorption region there 
is no propagation of energy at all, the vibrations are 
stationary and are spatially damped, as in the case of 
total reflection in the optically less dense medium. 

Concerning the ratio a of the amplitudes of the dipole- 
moment M and the field-strength E, it is to be noted that 
by (365) it is positive on the left of the absorption band, 
but negative on the right and vanishes entirely when 
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^ = 00 ; that is, for infinitely rapid vibrations the dipoles 
play no part in the wave, the body behaves optically 
like a perfect vacuum. 

§ 83. The laws of dispersion and absorption which 
have here been derived become a little more generalized 
if, in accordance with the idea already mentioned in 
§ 78, a small damping term kM, where k>0, is added to 
the left-hand side of the equation of vibration (355). 
The term hM y then becomes added in (360), and the 
same assumption (362) and (363), instead of leading to 
(365), leads to the values : 


n 2 — 1 = 4ttoc = . * — 

oj J — 1 + l/CtO 


(367) 


whereas the relationship (366) remains unaltered. 

The only difference from the result previously obtained * 
is that n is no longer either real or purely imaginary, but 
is always complex. If in place of n we now unite gener- 
ally n — i/c, then by (362) n again represents the real 
refractive index, k the absorption index ; and Eig. 24 is 
to be changed in that the curve for n no longer extends 
to infinity, then falling to zero and remaining zero in the 
whole region of absorption, but rather n remains con- 
tinuous everywhere, rising to a steep maximum after 
entering the absorption region and then falling continu- 
ously to zero, corresponding to so-called “ anomalous 55 
dispersion, whereas conversely the curve for the absorp- 
tion index k keeps close to the axis of abscissse in the 
two regions of normal dispersion and has appreciable 
values only in the region of the absorption band.’ 

§ 84. A further generalization of practical importance 
enters into formula (367) if the dispersing body contains 
different kinds of simultaneously vibrating dipoles whose 
number per unit volume we shall denote by N v N 2 , 
N* In place of (348) we then have the relationship : 


2M 1 . 


D -E 
4-7T 


(368) 
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which must he summed up over all the kinds of dipoles ; 
for each individual kind there is a vibration equation of 
the form (349), with or without a damping term, with 
constant coefficients co v a l9 and possibly also k v while the 
equations (351) and (359) become generalized to : 


E' - y EM X = E 


(369) 


E + 4w2 , M 1 = c 2 AE . . . (370) 

The same assumption (362) for the electric intensity of 
field and (363) for the dipole-moments M v M 2 , . . . with 
the amplitudes oc 1? oc 2 , a 3 . . . then leads, as in (367), to 
the general relationships : 


and : 


n 2 — 1 1 s? __ 

n 2 + 2 3 ^ cox 2 — cd 2 + ilcito 

A ( n ‘ 2 + gK 

7r0Cl 3 ( a )} 2 a > 2 + ikjcx)) 


(371) 

(372) 


which determine the complex index of refraction n — i/c 
and the amplitudes of the vibrating dipoles of all kinds 
for every wave-frequency co. In geometrical language, 
as shown diagrammatically in Pig. 24, this relationship 
states that every kind of oscillator in the spectral region 
gives rise to a particular absorption band, caused by its 
proper constants, within which anomalous dispersion and 
appreciable absorption take place, whereas outside the 
absorption band the dispersion is always normal. But 
there is this difference, that in every intervening region 
between two absorption bands the real refractive index 
increases from very small to very great values, whereas 
in the outer region on the left (small values of co) it is 
always greater than V € and in the outer region on the 
right (great values of cu) it is always less than 1. In 
general we may say that every absorption-band whose 
frequency exceeds the wave -frequency a >, increases the 
real refractive index, whereas every absorption-band 
whose frequency is less than oj 3 reduces the refractive 
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index. Hence a refractive index which is greater than 1 
always points to the presence of higher proper frequencies. 

In all the following considerations we shall assume as 
known the complex refractive index n — ik , by which the 
whole optical behaviour of the body is characterized. 

§ 85. Concerning the laws of refraction and reflection 
of a plane periodic wave at the surface of a dispersive 
body or of two dispersive bodies in contact, these may be 
derived directly from the formulse which hold for non- 
dispersive bodies. For if we reflect that the expressions 
which we have set up in (362) and (366) for the electric 
and for the magnetic field-strength of a wave which is 
progressing in the dispersive body are contained as 
special cases in the equations (8), and further, that the 
boundary conditions at the surface of separation, which 
express the continuity of the tangential field-components, 
are the same for dispersive as for non-dispersive bodies, 
it is clear that all the formulae that were derived earlier 
for the reflected and the refracted wave also retain their 
validity here. The only difference is that here we must 
in general insert for the refractive index n and conse- 
quently, by (20), also for the angle of refraction 6 1 a 
complex value. This also makes the constants / a, a, ju, v 
a v defined by the formulse (23) and (24), of the reflected 
and refracted wave complex, as well as in the case of 
total reflection. The physical significance of this circum- 
stance is expressed, as we saw earlier, by the occurrence 
of a sudden phase-change on reflection and refraction. 
Since this sudden change of phase is different for the 
wave which vibrates in the angle of incidence, the /-wave, 
from what it is in the case of the wave which vibrates 
perpendicularly to the plane of incidence, the g- wave, 
linearly polarized light is not linearly polarized after 
reflection, but elliptically polarized, and the measure- 
ment of the phase-difference of the two-wave components 
by means of a compensator, say a Fresnel rhomb (§ 27) 
or a sheet of mica (§ 64), leads to a relationship which, 
combined with the measurement of the refractive index, 
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gives us the two equations which are necessary if we wish 
to calculate the two constants n and k , by which all the 
optical properties of the body are determined. 

According to the general discussion of § 18 the intensity 
of radiation of the reflected and the refracted light is 
defined by the squares of the absolute values of the coef- 
ficients fjL and a, the quantities : 

M 2 = Pil and M 2 = Px . . . (373) 

directly representing the reflection coefficients of the two 
wave-components, while in the case of the refracted wave 
the transmission coefficient is, by § 10, calculated more 
simply from the reflection coefficient than from the 
quantities fi x and a v For the special case of normal 
incidence (9 = 0) we get from (23) : 

u — a 1 

p\\ = p±- YTe'x 

( n — t'/ci) — (n - ik) (% - n)~ + (« x — k ) 2 

(n x — iK X ) *f ( n — ik) (% -f ^) 2 + (/c x + /c ) 2 ' 

From this we see among other things that the reflection 
coefficient vanishes only if both n 1 -=^ n and also k x = k, 
that is, if there is no optical difference at all between the 
two bodies in contact. Hence a “ black body/’ that is, 
one which absorbs all light that falls on it, cannot have 
a plane surface. On the other hand, the reflection coeffi- 
cient can very well approximate closely to 1, namely 
always when any one of the four quantities n v n , k v k 
is very great compared with the rest. In this case we 
speak of cc metallic 55 reflection as in the case of an absolute 
conductor of electricity (III, § 92). 

§ 86. The method which we have here adopted to 
derive the laws of dispersion and absorption from the 
hypotheses assumed concerning the nature and the dis- 
position of the vibrating dipoles may be replaced by a 
totally different method, whose results may therefore 
simultaneously be used to test those just obtained. 
For instead of first setting up differential equations for 
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the electromagnetic phenomena that occur in a body by 
passing on from the assumption of an atomic constitution 
to that of a continuous constitution of the dispersive 
body and then forming particular integrals of these 
equations, there is essentially nothing to prevent our 
performing this transition later (indeed, we shall see 
that this deepens our comprehension of the phenomena 
in question) and first considering directly the electro- 
magnetic waves that are emitted by the individual 
vibrating dipoles. On this view there is only one medium 
in which the waves propagate themselves, namely in the 
vacuum; in this vacuum there are embedded a great 
number of vibrating dipoles, each of which forms the 
centre of a spherical wave which propagates itself accord- 
ing to known laws independently of the remaining dipoles. 
For by III, 396 the external electromagnetic field pro- 
duced by a dipole vibrating with the moment f(t) is 
uniquely determined for all points of space and for all 
times. So we arrive at a particular solution of the 
problem, corresponding to the equations (362) and (363), 
if we assume the vacuum to contain firstly a periodic wave 
advancing with the velocity c in the direction of the 
a;-ax:is and vibrating in the direction of the 2 -axis with 
the frequency — this is the “ primary wave ” — and 
secondly a great number of spherical waves, te secondary ” 
waves, which are emitted by the individual dipoles in 
accordance with the vibrations which they emit with the 
same frequency ; these secondary waves superpose them- 
selves on each other and also on the primary wave. The 
result is the “ effective ” wave, whose field-strengths E 
and II are formed directly by vectorial addition of the 
corresponding field-strengths of the primary and the 
secondary waves. Since the dipoles lie very closely 
together, we may replace the addition by an integration. 

Account is taken of the interactions of the dipoles 
among themselves and with the primary wave in that the 
equation (349) now holds, just as before, between the 
moment of vibration of a dipole and of the field-strength 
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E' which excites it, the equation (351) or (369) again 
holding for the relationship between E and E ' . 

It is clear that these relationships give rise to a perfectly 
definite law for the propagation of an effective wave in 
the body in question ; and if there is no inherent contra- 
diction in any part of the theory, this law must be identical 
with that deduced earlier, although here we have made 
use neither of the differential equations of dispersion nor 
of the conception of electric induction, but only of the 
laws of waves in a pure vacuum. 

This postulate is fulfilled in actual fact, as we can 
easily show by direct calculation. It must be par- 
ticularly emphasized that compounding the secondary 
spherical waves due to the dipoles vectorially must give 
and does give two results — firstly a wave which advances 
with the velocity c and is equal and opposite to the 
primary wave, which also advances with the velocity 
c and hence just neutralizes it; secondly a wave which 

advances with the velocity -, where n is determined by 

(365) or (371) respectively. In all directions other than 
that of the a;-axis the secondary waves destroy each other 
by interference. 

We see that this new method of regarding the problem 
gives us a direct understanding of the nature of dispersion. 
But its physical significance goes still further. For it 
also leads to a new derivation of the laws of refraction 
and reflection. Let us imagine, for simplicity, a plane 
periodic wave to fall from a vacuum on to the plane 
surface of a dispersive body at the angle of incidence 6. 
This wave, which we may again call the primary wave, 
passes completely undisturbed through the body; for 
it propagates itself with the velocity c in the body just 
as in the vacuum. According to this view, then, there 
is no meaning in speaking of a surface of separation 
or contact. The influence of the body makes itself felt 
solely by the simultaneous vibration of its dipoles, and 
we may again regard the effective wave-process as the 
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superposition of the primary wave on the secondary 
waves emitted by the vibrating dipoles of the body, 
by making the appropriate assumptions for the moments 
of the dipoles, which must of course satisfy the equation 
of vibration (349). The result of the calculation again 
confirms in every detail the laws derived in the preceding 
section by a totally different method. For the result of 
superposing the secondary spherical waves is firstly a 
plane wave which completely destroys the primary wave 
by interference ; and secondly a plane wave which 
advances in the direction of the ray refracted according 
to Snell’s law, and thirdly a wave which is reflected 
back into the vacuum at the angle 9. The amplitudes 
of the waves also come out as above calculated. 

Although this method of deriving the laws of reflection 
and refraction is more complicated mathematically, it 
is more significant physically than the former method. 
For since it involves no boundary condition it gives a 
direct and complete explanation of the phenomena of 
reflection and refraction, which, according to it, are 
nothing but the combined effects of spherical waves that 
have been emitted and interfere in a particular way. 

§ 87. All the preceding discussion refers exclusively to 
such particular solutions of the field-equations as corre- 
spond to singly periodic waves. If we now turn more 
generally to the investigation of non-periodic plane waves 
which advance in the ^-direction within a dispersive 
body, it must be emphasized that the conception of 
velocity of propagation here loses completely the general 
meaning which it has for non-dispersive bodies. For in 
a dispersive body it is not always possible to represent 
the wave-function as a definite function of a single 
argument which depends only on x and t, as was done in 
equations (8) ; or, in other words, the wave does not 
advance unchanged, but becomes deformed in the course 
of time, and it is impossible to recognize with certainty 
any selected point of the wave again at a later time. 

From this it follows at once that we are at liberty to 
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define an infinite number of different kinds of velocities 
of propagation according to the characteristic feature 
that we fix on as the wave advances. For example, we 
can fix on the value of the field-strength E y at a point 
x at any time t and enquire what is the neighbouring 
point x + dx at which the field-strength has the same 
value at the time t + dt (E y = const.). This gives the 
following value for the velocity of propagation of the 
field-intensity E y : 


dEy 

dx __ dt 

dt dEy 

dx 


(375) 


a value which is always constant for non-dispersive 
mediums, but for dispersive mediums it is constant only 
for singly periodic waves, and hence in general has no 
deeper significance. The same applies if we consider the 
velocity of propagation of a maximum or a minimum 
of the field-strength. Its value is obtained from (375) 
dE 

if we substitute —■ in it for E y and is of course different 
in general from (375). 

It follows from this discussion that it is altogether 
impossible to define for a dispersive body a velocity of 
propagation which is constant in space and in time, 
so long as we make no assumption about the special 
form of the wave. It is easy to see, moreover, that for 
different forms of the wave the velocity of propagation 
may be defined in quite different ways. We shall see this 
confirmed in several particularly important cases. 

If we now enquire into the general integral of the wave 
equation (360) and the associated equation of vibration 
(361) we can find an essentially adequate solution by 
the same method as was successful in the case of the 
problem of total reflection, in § 14. For we may express 
any arbitrary function over any arbitrarily great region 
of its argument as a Fourier series, and then every indi- 
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vidual term of the Fourier series represents a s ingl y 
periodic occurrence which obeys the laws of propagation 
found above. Since the differential equations (360) 
and (361) are linear and homogeneous, the individual 
processes become superposed on each other and hence 
it follows that every wave which advances in a dispersive 
body may be regarded as composed of a number of singly 
periodic waves each of which advances unchanged with 

the velocity of propagation ~ which is characteristic of it 

and which is constant in space and in time. But since 
n is different for the individual waves, the phases of 
the individual partial wave become displaced relatively 
to one another and so cause the deformation of the total 
wave. 

§ 88. We shall now fix our attention on the special 
case where the dispersive body, which we shall suppose 
stretches from x = 0 to x — oo , is perfectly unexcited 
at first and that from a certain moment of time onwards, 
say t — 0, a singly periodic wave of frequency a> is con- 
tinuously incident on its surface in the normal direction. 
With what velocity will the disturbance propagate itself 
into the body ? We may be tempted to assume at first 
sight that the wave simply propagates itself into the body 

with the velocity as is the case with a permanently 

periodic wave. But this is not so. For here we are 
concerned with a non-periodic wave which becomes 
deformed as it advances. The complete solution of this 
problem has been given by Sommerfeld, 1 who resolved 
the waves into a Fourier integral. But we can find the 
answer to the above question without making special 
calculations if we apply the line of reasoning which was 
fully discussed in § 86. 

The primary wave also advances in the dispersive 
body with the velocity c. The secondary waves due 

1 A. Sommerfold, Weber-Festsehrift, B. G. Teubner, p. 338, 
1912. 


O 
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to the vibrating dipoles do the same. Rut the result 
of this joint action is not the same as in the ease of a 
permanently periodic wave. For the dipoles do not 
vibrate periodically; they were at rest initially and, on 
account of their inertia, they will only gradually be made 
to vibrate, the later the greater their abscissae x, and these 
vibrations, strictly speaking, will acquire their constant 
character only after an infinitely long period of time 
and so give rise to the singly periodic wave-motion whose 

c 

velocity of propagation is -. So long as this is not the 

case the primary wave will not be neutralized by the 
secondary waves and so it manifests itself as a disturbance 
which advances in the dispersive body with the velocity 
c. Since on the other hand the secondary waves, no 
matter what their form may be, also propagate themselves 
with the velocity c it follows that the desired velocity 
with which the head of the wave advances into the body, 
the “front velocity / 5 is always equal to c. In contrast 

with the front velocity is the phase- velocity / 5 -, which 

7b 

a permanently singly periodic wave has; in future we 
shall denote it by u : 

\ = « ( 376 ) 

We must not regard it as an inherent contradiction 
to the theory that the phase-velocity u may also be greater 
than the front velocity c (n<l). For the wave is not 
singly periodic and the head of the wave that penetrates 
into the body has no constant phase. Hence there can 
be no question of the head of the wave which advances 
with the front velocity being passed by the hinder part 
of the wave which advances with the phase-velocity. 

It further follows from this that the velocity of propaga- 
tion of a light-signal emitted into the body, the “signal 
velocity / 5 can at most be equal to c. This maximum 
value is attained in an ideal detector, that is, a receiver 
which reacts to even the smallest disturbance. Other- 
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wise the signal velocity must be made correspondingly 
smaller according to the sensitivity of the receiver. The 
value of the phase-velocity does not come into question 
at all in the matter ; for it is not possible to send signals 
by means of singly periodic waves. 

§ 89. We shall now fix our attention on a non-periodie 
wave of a form which is in principle still more important, 
namely a wave which is nearly singly periodic. For this 
case is realized in nature in all so-called monochromatic 
rays which are actually produced (§ 16). The general 
expression for a wave of this kind which is nearly periodic 
and which advances in a dispersive body is, by (362) : 


where : 


E„ 




OJJV — CO 1 ^ ^ ^ 
COi 


(377) 

(378) 


Here m is the order number with which the frequency 
aim increases, the total number N is arbitrarily great and 
Cm is any complex constant. Let us first consider the 
form of such a wave, that is, E,, as a function of a; at a 
definite time t, say represented by a curve with the 
abscissa x. This curve is formed by the superposition 
of a number of simple sine curves of very nearly the 
same wave-length. The result is, as emerged from our 
exhaustive discussion in § 16, a nearly periodic curve of 
the same wave-length, but with a slowly and in general 
irregularly changing amplitude and phase. 

Let us consider the conditions a little more closely. 
The decisive features for the amplitude and phase of 
the resultant wave are the phase-differences with which 
the individual singly periodic partial waves meet. If 
we fix on any one of them of frequency a> and phase - 

velocity then the phase of any other partial wave of 

')b 

Q 

frequency a/ and phase-velocity can, by (377), for all 
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points x and all times t 3 be set equal to the phase of the 
first wave plus the phase-difference : 

(a / — a>)t — - (c o'n' — con) -f - 6' ~ d . . (379) 

where 6' and 6 denote the angular arguments of the 
complex constants O' and C. So long as this phase- 
difference remains appreciably unchanged for all com- 
binations of any two partial vibrations, the partial waves 
interfere in the same way and the resultant wave has a 
constant amplitude and phase. Now if the time t is 
assumed constant the phase-difference changes only with 
the abscissa x. If we pass from the point x to a neigh- 
bouring point x x >x the phase-difference (379) changes 

by : 

- Xl ~ ~ - ( an ’ - wri) . . . (380) 


Hence so long as the two points lie so close together 
that for all pairs of partial vibrations : 



the resultant wave forms between the two points a simple 
sine curve of wave-length : 


or, respectively. 


which, by (378), are appreciably the same. But if we 
pass so far along the axis of abscissae with the point x t 
that the relationship (381) ceases to be valid, the inter- 
ference effect of the partial waves changes and the result- 
ing wave acquires a new amplitude and a new phase- 
constant. Accordingly we may say : the form of the 
resultant wave E u at any time t may be regarded in 
general as compounded of a number of successive groups 
of simple sine -curves. Every single group stretches over 
a region of abscissae for which the relationship (381) is 
fulfilled ; a characteristic of the group is a definite phase- 
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difference (379) for any two partial vibrations, and con- 
sequently also a definite amplitude and a definite phase- 
constant. All these quantities vary from group to group. 
The number of wave-periods that belong to a group is, 
by (332) : 


x i — x (aq — x ) . nco 
A 2ttc 


(383) 


where x and 
the group, 
relationship : 


x 1 denote the initial and the final points of 
This number may be very large; for the 


(aq — x ) . nco _ 

— — L ~> > 1 

27TC 


(384) 


is, on account of (378), certainly compatible with the 
relationship (381). 

Let us enquire further into the changes that occur in 
the course of the time t firstly when we keep a definite 
point x in space fixed. Then, by reasoning which is fully 
analogous to that given above, we find that within an 
interval of time t 1 — t which is so small that for every 
pair of partial vibrations : 


U -t<< 


co 


(385) 


the resultant field-strength E y executes pure sine vibra- 
tions of frequency a>, but that at later times t a change in 
the amplitude and the phase-constant occurs. 

The process is a little easier to picture if we fix our 
attention not on a definite point x of space but on a 
definite phase of the resultant wave, say the peak of the 
crest of a wave. This crest advances with the phase- 

velocity ^ and retains its value and its velocity of propaga- 
tion along a certain distance which may include many 
wave-lengths. But after that it will slightly change its 
value and also its velocity of propagation. The latter 
will always be nearly equal to the phase-velocity, but the 
value of the crest of the wave can increase or decrease 
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by an arbitrary amount. In other words, any definite 
wave-crest on which we fix our attention does not remain 
in the same group in the course of time, but passes over 
into a neighbouring group. 

This circumstance leads us to enquire with what 
velocity a definite selected group advances. The answer 
is contained in the above given property of a group that 
for it the phase-difference (379) for every combination 
of any two partial vibrations retains a definite invariable 
value. Hence if we set the value (379) for %, t equal to 
that for x v t v we get : 

x i ~ x __ ( a > / ~~ C0 ) c .... (386) 

q — t co'n' — con 


which is independent of x and t. This velocity, which is 
constant for all points and all time and which may, on 
account of (378), also be written in the form : 


c 


clou 

d(ajn) 


(387) 


represents the required group -velocity, 
with the phase- velocity u, according 
relationship : 



v dco 


It is connected 
to (376), by the 


(388) 


The difference in u and v of course depends on the 
circumstance already pointed out above, that a definite 
selected wave-crest advancing with its phase-velocity 
changes from group to group in the course of time. 

It is only in the case of dispersionless bodies =* Q 

that v = u, and every wave-crest remains permanently 

in its group. In the case of normal dispersion 


dco 


> 0 


we have v< u, that is, in the course of time a definite 
w T ave-crest moves out of its group forwards into the next 
group and so catches up one group after another. In 
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the case of anomalous dispersion the reverse occurs. But 
we must bear in mind that a group does not form a 
configuration which is permanent for all time. Eor since 
the differential-quotient in (386) is not, strictly speaking, 
equal to the differential quotient in (387) but has a 
different value for every combination of two partial 
vibrations, in a fairly great interval of time t 1 — t the 
distance x x — x ' traversed will be appreciably different 
for different combinations, that is, the group will fall 
apart. Hence in rather great intervals of time a gradual 
change also occurs between the groups in that every 
group eventually splits up and fuses partly with the 
preceding group and partly with the following group. 
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GEOMETRICAL OPTICS OE NON-HOMOGENEOUS 
BODIES. RELATIONSHIPS TO QUANTUM 
MECHANICS 

§ 90. An optically non-homogeneous body is charac- 
terized in having its refractive index n dependent in 
some definite way on the space-co-ordinates x, y, z. We 
shall also assume that the body is dispersive, that is, we 
also make n depend on the frequency oj. But we shall 
assume that it is isotropic, so that at a definite point 
and for a definite frequency n is the same in all directions. 
The laws of propagation of light-waves for such a body 
are in general very complicated, but an essential simplifi- 
cation may be achieved by restricting ourselves to the 
region of geometrical optics or ray-optics which, by § 28, 
embraces those phenomena for whose representation the 
laws which hold for plane waves suffice. 

Let us again begin with singly periodic vibrations. 
Then every wave-function has the phase : 

cot - (f>(x, y,z) -© . . . . (389) 

where w denotes the frequency, and </> a certain function 
of position, the “ eikonal ” ; the conception of <f> is 
generalized from (362). 

If we define a wave-surface as a surface whose points 
at a definite time t are all in the same phase 0, then a 
wave-surface is represented by the equation : 

<f>(x,y,z)= const (390) 

The family of wave-surfaces which results when we allow 
the constant to run through all its values is the same for 
all times and hence is fixed in the body. If we fix on a 
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definite phase 0, then, by (389), it will be found at the 
time t in the surface <f> = cot — 0, and so as the time 
increases it passes uniformly from smaller to greater values 
of (/>. If ? on the other hand, we fix on a definite point 
x, y , z in space, the phase at it increases uniformly with 
the time : 



(391) 


The characteristic theorem of geometrical optics is that 
every sufficiently small portion of a wave-surface can be 
regarded as a plane wave which propagates itself inde- 
pendently in the direction of its normal with the phase- 

velocity u = - which it has at the point in question. 

Since the body is isotropic, the direction of the normal is 
simultaneously that of the ray; and, if any particular 
wave-surface is known, we obtain an infinitely near -wave- 
surface by marking off from every point of the given 
wave-surface the distance u . (It in the direction of the 
normal and then connect the end-points of all these 
distances together. Thus if ds is the normal distance 
between the two surfaces, we have : 

ds = udt. 

Since, on the other hand, the phase 0 remains the 
same during the propagation of the wave — for u is the 
phase-velocity — we have, by (389) : 

codt — • d<f> — 0 

and this, combined with the preceding equation, gives 
the following value for the differential coefficient of <j> in 
the direction of the normal : 

d<f> ___ co __ no) _ d& (392) 

ds u c ds 

By means of this differential equation, the so-called 
eikonal equation, it is possible to calculate from a single 
wave-surface all the remaining wave-surfaces, since n is 
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given as a function of x, y, z. The whole of geometric 
optics is in essence contained in this equation. For if 
the system of wave-surfaces is known, its orthogonal curves 
represent the rays which propagate themselves through 
the body. 

§ 91. The laws of geometric optics, besides allowing 
themselves to be expressed by equation (392), may also 
be represented by Fermat’s principle, the derivation of 
which may be made clear in a w r ay similar to that used 
in § 37 for a special case. 

According to Fermat’s principle, the time which a 
definite phase requires to arrive with its characteristic 
velocity u from a definite point P to a definite other point 
Q of the body along the path of the ray which goes from 
P to Q is less than that taken along any other line 
connecting P and Q, or, symbolically : 

8 J ™ = 0 = S J nils . . . . (393) 

If we substitute for the variations in the equation : 


8 n 


Ms 


dn- , dn~ , 
fo 8x + ty 8y + 


S z 




and then, as in I, § 108, integrate by parts so as to 
reduce the variations S dx s 8dy, Sdz to 8x, 8s, after- 

wards setting the coefficients individually equal to zero, 
the differential equations for the path of the ray from P 
to Q is obtained in the form : 





3 n m d n . an 
~dx * dy ’ dz 


(394) 


These equations, by which the path of the rays, can be 
calculated if the starting-point and the direction of the 
ray are given, also contain all the laws of geometric 
optics. In content they are completely identical with 
equation (392) for the wave-surfaces, as can also be seen 
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directly if we imagine the construction described in § 90 
to be carried out for a wave-surface adjacent to a given 
wave-surface and if we calculate the resultant change in 
the direction of the normal, that is, the difference between 
the normals of the new and the original surface. We 
then arrive at exactly the relationship (394). 

§ 92. In order to be sure of our subsequent procedure 
we shall first formulate a little more accurately the con- 
ditions under which geometric optics may be applied. 
In general the differential equations of dispersion (359) 
and (355) hold. The former contains only the universal 
constant c, in the latter we assume for a non-homogeneous 
body that the constants co 0 ' and a are given functions of 
the space-co-ordinates x, y, z. 

The two differential equations are then satisfied by a 
generalization of the assumptions (362) and (363) made 
for plane waves in a homogeneous body : 

E = if, M = cLifs . . . . (395) 

where a is now to be a certain function of position, while 
i//, the wave-function, depends in some way on phase and 
time. 

By substituting in (355) and (359) we obtain : 

g /.if + ci)' 0 2 oLif = -^-if (396) 

and : 

f + 4:rrOiif — C 2 A*/r 
or 

^ = |.A . . . (397) 

where n is again related to oc by (365). 

If wc now introduce the assumption that the process is 
singly periodic in time, that is : 

f - e*(«*-*> • (398) 

where the function of position <f> , being a generalization 
from (389), may also be complex, then; 

if = — co^ip 


. (398a) 
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and the relationship (396) leads to the same value of n 
as the relationship (365) deduced for plane waves in 
homogeneous bodies. So the whole problem is in general 
reduced to a differential equation (397) in the wave- 
function j ft. 

Substituting the equation (398) we get : 


c 2 



\dzJ 


or, if ds again denotes the element of length of the normal 
to the surface ijj = const. : 


This relationship becomes the same as the fundamental 
equation (392) of geometric optics if and only if the 
term in A <j> can be neglected. Since the equation (399) 
is complex, two conditions follow from this. Firstly the 
differential coefficients of the imaginary part of <f> must 
be very small or, what comes to the same thing, according 
to (398), the amplitude of the wave-function must be 
only slowly variable in space ; secondly, the term in the 
second differential coefficient of the real part of <f> must 
be small compared with each of the other two terms of 
the equation (399). In view of (382) this means that 
the radii of curvature of the wave-surfaces must be great 
compared with the wave-length — the same postulate as 
we set up earlier in § 28. Or, the radii of curvature of 
the rays must be great compared with the wave-length. 
So long as these conditions are fulfilled, the laws of 
geometric optics may be applied as a first approximation. 
But as soon as they are transgressed at any point the 
exact wave-equation (397) enters into force in place of 
them. 

§ 93. We now turn from singly periodic rays to the 
natural monochromatic rays, that is, nearly singly periodic 
rays of the same kind as those we investigated in § 80. 
We saw there that such a ray may be regarded as a 
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progressive series of consecutive groups of singly periodic 
rays whose velocity of propagation v is characterized by 
the condition that the difference 0' — © of the phases 
of any two partial vibrations a>' and o> remains constant. 
Thus by (389) : 

(a/ — a>) . t — (<£' — <f>) = const, 
or: 

W -g)*-o. 

From this we get, taking (392) into account : 

(Is to' — CO 

at n a) 7io) 

which agrees perfectly with (386), (387) and (388). 

Since the group velocity v is different from the phase- 
velocity u, a point which advances with the group- 
velocity will in the course of time change its phase, and 
by (389) we have : 

(IS 00 00 

tit dt ds V 

and by (391) and (392) : 

d ® = (400) 

(it U v lU 

Combined with (376) and (387) we get from this : 

dn 

d® — — (401) 

dt dn n 
060 60 

The differential coefficients have here been written with 
round 3’s because n also depends on position. 

For a non-dispersive body j— = 0, so that in this case 

the phase 0 remains constant, that is, the group-velocity 
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coincides with the phase-velocity. This holds for noil- 
homogeneous as well as for homogeneous bodies. 

§ 94. The fundamental equations of geometric optics 
above derived for a non-homogeneous dispersive isotropic 
body for the case of a nearly singly periodic ray exhibit 
in their formal structure a remarkable analogy with the 
fundamental equations of classical mechanics for a free 
point-mass which moves in a given statical conservative 
field of force (I, § 49) — an analogy which expresses itself 
in the circumstance that a wave-group (restricted to a 
small portion of space) moves in the direction of its ray 
according to the same laws as the point-mass moves in 
the direction of its orbital curve. It is, indeed, possible, 
by appropriately allocating the optical constants in the 
first case to the mechanical constants in the second case, 
to make the optical and the mechanical equations 
identical. This will now be shown. 

For this purpose we first consider our two equations 
(391) and (401), which specify the time-change of phase 
& for a stationary point and for a point moving along 
the ray with the group -velocity v. We compare these 
equations with the two corresponding equations of 
mechanics for the integral of action IF, considered as a 
function of x, y , z , t, namely (I, 420) : 


dW 

dt 


= -E 


(402) 


for a stationary point, and (I, 419) : 


dW 

dt 


= E — 2U 


(403) 


for a point which moves on its orbital curve according to 
the laws of motion. Here E is the energy and E — 2U is 
the kinetic potential or Lagrange’s function, namely, the 
difference of the kinetic energy : 

E — U = ~mv 2 .... (404) 


and the potential energy U, which we consider as a given 
function of the co-ordinates x , y, z. 
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In addition to the formulae for the time derivatives of 
© on the one hand and W on the other we have those 
for the space derivatives, namely, equation (392) for © 
and the equations I (417) for W, which express that 
the space gradient of W is equal to the momentum of the 
moving point-mass, that is, that the direction of the 
velocity at every moment is perpendicular to the surface 


W = const, and that the value of the gradient ~ 

OS 

represents the product of the mass and the velocity : 


dW 

— == mv 

ds 


(405) 


If wo now identify the velocity v of the moving point- 
mass in size and direction with the velocity of propaga- 
tion of the wave-group designated by v earlier, and if we 
reflect that the ray-direction is perpendicular to the wave- 
surface 0 = const., we shall be led to set the functions 
IF and 0 proportional to one another. It then follows 
by comparison with (391) and (402) that the optical 
constant to is proportional to the mechanical constant E. 
Since these two quantities have different dimensions, the 
factor of proportionality is not a pure number. Therefore 
we set : 

E = A • (406) 

Ztt 


and interpret h as a certain constant whose dimensions 
are the product of an energy and a time. Then (391) 
and (402) become identical, if we further set : 

W = . . . . (407) 

Ztt 

This also establishes the relationship between the equa- 
tions (405) and (392) : 

3 W __ h 3© _ nojh __ nE 
mV ds 2 tt ds 2 itc c 


( 408 ) 
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and from this we get, if v is calculated from (404), the 
following value of the refractive index : 


n = yj • V 2m{E 

Mi 


V) 


(409) 


The dependence of the refractive index on the function 
of position U is caused by the non-homogeneity of the 
body, the dependence on the energy constant E or, 
respectively, to, is due to the dispersion of the body in 
which the wave-group in question is propagated. 

It is now easy to convince ourselves of the fact that 
through the relationships which have been introduced 
between the mechanical and the optical quantities the 
two systems of equations become fully identical. In 
particular the equation (387) for the group-velocity : 


1 d{a>n) 

C doj 


1 d(En) 
c 3 E 


. (410) 


leads identically to the value of v in (404), if we substitute 
from (409), and in the same way (401) transforms directly 
into (403). 

So we may perfectly generally enunciate the theorem : 
according to the laws of classical mechanics a point-mass 
m moves in a field of force of given potential energy U 
exactly like a singly periodic ray in a non-homogeneous 
dispersive isotropic body whose index of refraction is : 


2ttC 
oo h 


Mt 


u 


. (411) 


Just as the motion of the point-mass is determined by 
the initial position, the initial direction and the energy E, 
so the propagation of the ray corresponding to it is 
determined by the initial point, the initial direction and 
2ttE 


the frequency to = 


h 


Although a far-reaching agreement is established in 
this way, between the laws of classical mechanics and 
those of geometrical optics, the agreement is not perfect. 
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For there is a fundamental difference in the following 
respect. Classical mechanics demands absolute validity ; 
it is complete in itself and is not bound by assumptions 
concerning orders of magnitude. Of geometrical optics 
we know, however, that it does not hold generally, but 
by § 92 only if the wave-length (382) corresponding to 
the frequency to of the ray is small compared with the 
radius of curvature of the path of the ray. As soon as 
this condition is infringed the laws of geometrical optics 
become inadequate and we have to fall back on the 
general wave-equation (397) or (399). 

§ 95. We would be able to rest satisfied with these 
results if classical mechanics accounted for all the dynami- 
cal processes that occur in nature. But this is not the 
case. 

Rather the facts of experience compel us to assume 
that in the case of very small quick movements certain 
of the circumstances arise which are quite foreign to the 
basic assumptions of classical mechanics; for example, 
in the motion of an electron about an atomic nucleus 
certain perfectly definite values of the energy E appear 
to be favoured — a theorem which is totally foreign to 
classical theory. 

Hence if it has been shown to be necessary to change 
classical mechanics, it can only be a question of generaliz- 
ing the laws hitherto established for it; for these laws 
would then retain their significance under all circum- 
stances. Now we must bear in mind that classical 
mechanics can never be generalized from within itself ; 
for, like every system that is perfect in itself, it bears the 
stamp of inner completeness and conceals no germs which 
allow of further development. Rather, the impulse must 
come from without. A beautiful example of this is given 
by the theory of relativity, which, starting from optics, 
has penetrated into mechanics and has added a totally 
new side to it. 

In the present case, too, help has come from the direc- 
tion of optics. L. de Broglie first adduced the long- 
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familiar analogy, discussed in the previous section, 
between classical mechanics and geometric optics to 
account for quantum phenomena, and E. Schrodinger 
then developed the idea further, showing that the exten- 
sion which the laws of classical mechanics must incur in 
order to adapt themselves to reality is of exactly the 
same type as that which must be applied to geometric 
optics when the limits within which it holds are over- 
stepped. In other words, the laws of the new mechanics 
are found simply by retaining throughout the analogy of 
mechanics with optics and accordingly modifying the 
laws of the old mechanics so that the basic difference 
mentioned at the end of the previous section is eliminated. 
The equations (406) and (411), which are characteristic 
for the analogy between mechanics and optics, then remain 
valid beyond the range of validity of geometric optics, 
and out of classical mechanics, ray mechanics, we thus 
get wave mechanics, which contains the former as a 
special case. Its characteristic feature is that it no longer 
allows the motion of a point-mass to be represented by 
the motion of a geometrical point, but that the point-mass 
is in a certain sense resolved into a number of waves of 
a definite kind, just as a light-ray, strictly speaking, does 
not allow itself to be represented by a single curve, but 
arises only from the combined action of optical waves. 

The first great success of this idea is shown by intro- 
ducing the wave-equation (397) into the new mechanics. 
This gives with the value of (n) in (411) : 

8rr 2 7n f a>Jl tt\ y A f , 

= • • • ( 412 ) 

or, if we substitute the expression (39 Sa) for Yfs : 

A# + 8 A(g _ J7 ) # _ 0 (413) 

for which, by (406), we may also write : 

+ = 0 . . (414) 
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The peculiarity of this differential equation of Schro- 
dinger consists in the fact that the coefficient (E — U) 
can also become negative in some circumstances. If the 
constant E is taken so great that for all points of the 
space — and the wave has no spatial limits — (E — U) > 0, 
there are finite and continuous solutions of the equation 
for all values of E. But if E is not sufficiently great to 
make the coefficient positive for even the greatest values 
of U, there are finite and continuous functions iff every- 
where, which satisfy the equation (414) only if E has 
certain definite values E v jE 2 , E s , — , the so-called proper 
values ( Eigenwerte ). Corresponding to these values we 
have the proper functions (. Eigenfunktionen ) if/ v 0 a# 
ifj$ . . . Perhaps this startling result appears less striking 
if we consider that for negative values of (E — U) the 
value (409) of the refractive index becomes imaginary 
and so geometric optics loses its meaning entirely. 

But the fact that gives the differential equation (414) 
its great importance is that the proper values for E 
calculated from it agree exactly in all cases with those 
proper values which had already been calculated earlier 
by Heisenberg, Born and Jordan from the equations of 
quantum mechanics in matrix form; these equations 
had been based directly on experimental facts and had 
been developed independently of any particular physical 
assumptions. This subsequent agreement between results 
which had been obtained by two totally independent 
methods is a definite indication of their physical signifi- 
cance, and hence there can be no doubt that Schrodinger’s 
differential equation, in virtue of its close relationship to 
classical mechanics, deepens our insight into the nature of 
quantum phenomena. 

It is true that much still remains to be done before we 
shall be clear about the physical nature of the wave- 
function if; and before we can satisfactorily answer all the 
questions that thrust themselves upon us. 

§ 96. We shall conclude by showing how it is possible 
to calculate the numerical order of magnitude of the 
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universal constant h in (406) by using the hypothesis 
that was introduced in the preceding paragraph. In tin 
first case it is clear that if we assume h to be infinite! 
small, wave-optics becomes merged into geometric optic, 
and hence quantum mechanics becomes identical with 
classical mechanics. For with this assumption it follows 
from (406) that for every value of E there is an infinitely 
great value of co, and hence by (382) an infinitely small 
value for the wave-length A ; so the condition that the 
wave-length must be small compared with the radius of 
curvature of the orbital curve is always fulfilled. » 

From this we see that the value of h will come out th 
greater the greater the radii of curvature of the orbit* 
curve at which the deviations from classical mechanic 
become appreciable. We shall now take as the basis ( 
our estimate the assumption that for an electron whic 
revolves around a stationary atomic nucleus which bear 
a single positive charge the deviation from classica 
mechanics becomes appreciable when its orbital radio 
shrinks to atomic size. 

We have by classical mechanics in general that if 
denotes the orbital radius, v the orbital velocity, e th 
charge, m the mass of the electron, then in the case 
uniform circular motion the centripetal force is equal 
the attractive force of the nucleus, thus : 


and : 


mv z 

f? 



On the other hand, we get for the wave-length A 
which is to be allocated to the motion here considered, 
by (382), (411), (406) and (404) : 


h h 

V%m(E ~U) mv 
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and by (415) : 

= (416) 

e v m 

t ms the ratio of the wave-length A to the radius of 
curvature r is : 

* = h (417) 

r z^/mr 

If this ratio is to assume an appreciable value the 
umerator and the denominator of the fraction must be 
the same order of magnitude, that is : 

h^e\/ mr (418) 

If, in accordance with our above assumption, we set r 
ual to the value 10~ 7 cm., then in electrostatic units : 

e = 4*77 . 10~ 10 (erg . cm)-, m — 9*02 . 10“ 28 grm. 

7t~4*5 . 10~ 27 erg. sec. . . . (419) 

The more exact value of the quantum of action is : 
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